Theoretical description of halogen bonding - an insight based on the natural orbitals for chemical valence combined with the extended-transition-state method (ETS-NOCV)

An innovative deep eutectic solvent: chalcogen bonding as the primary driving force

Chalcogen bonding (ChB) interactions have drawn intensive attention in the last few decades as interesting alternatives to hydrogen bonding. The applications of ChB were mostly centered on the solid state and have rarely been explored in solution. In this work, a novel strategy for forming ChB-based deep eutectic solvents (DESs) was exploited. We set forth the formation, physicochemical properties, and interaction sites in detail. This work not only provides a new idea to design DES systems but also to exploit the potential application of ChB complexes.

Occupied-Virtual Orbitals for Chemical Valence with Applications to Charge Transfer in Energy Decomposition Analysis

In this article, we introduce the occupied-virtual orbitals for chemical valence (OVOCV). The OVOCVs can replace or complement the closely related idea of the natural orbitals for chemical valence (NOCV). The input is a difference density matrix connecting any initial single determinant to any final determinant, at a given molecular geometry, and a given one-particle basis. This arises in problems such as orbital rearrangement or charge transfer (CT) in energy decomposition analysis (EDA). The OVOCVs block-diagonalize the density difference operator into 2 × 2 blocks, which are spanned by one level that is filled in the initial state (the occupied OVOCV) and one that is empty (the virtual OVOCV). By contrast, the NOCVs fully diagonalize the density difference matrix and therefore are orbitals with mixed occupied-virtual character. Use of the OVOCVs makes it much easier to identify the donor and acceptor orbitals. We also introduce two different types of EDA methods with the OVOCVs and, most importantly, a charge decomposition analysis method that fixes the unreasonably large CT amount obtained directly from NOCV analysis. The square of the CT amount associated with each NOCV pair emerges as the appropriate value from the OVOCV analysis. When connecting the same initial and final states, this value is identical to the CT amount obtained from the independent absolutely localized molecular orbital (ALMO) complementary occupied-virtual orbital pair (COVP) analysis. The total, summed over all pairs, is also exactly the same as the independently suggested excitation number, as proved herein. Several examples are presented to compare NOCVs and OVOCVs: stretched H2+, a strong halogen bond between tetramethylthiourea and iodine, coordination of ethene in Zeise's salt, and binding in the Cp3La···C≡NCy complex.

Activation of perfluoroalkyl iodides by anions: extending the scope of halogen bond activation to C(sp3)-H amidation, C(sp2)-H iodination, and perfluoroalkylation reactions

A simple, efficient, and convenient activation of perfluoroalkyl iodides by tBuONa or KOH, without expensive photo- or transition metal catalysts, allows the promotion of versatile α-sp³ C-H amidation reactions of alkyl ethers and benzylic hydrocarbons, C-H iodination of heteroaryl compounds, and perfluoroalkylations of electron-rich π bonds. Mechanistic studies show that these novel protocols are based on the halogen bond interaction between perfluoroalkyl iodides and tBuONa or KOH, which promote homolysis of perfluoroalkyl iodides under mild conditions.

Explaining the High Catalytic Activity in Bis(indenyl)methyl Zirconium Cation Using Combined EDA-NOCV/QTAIM Approach

The main purpose of this study is to elucidate some discrepancies already observed in the catalytic activity values of some zirconocene methyl cations. The EDA-NOCV scheme was employed for a theoretical description of the interactions between an ethylene molecule and five catalysts of zirconocene methyl cation. The nature of the chemical interactions has been elucidated through the QTAIM topological analysis. The steric hindrance due to the ligands was evaluated qualitatively by means of an IRI-based analysis and quantitively through Fisher information. The findings prove that the indenyl ligand seems to favor the orbital interaction between the ethylene molecule and the metal centre of zirconocene methyl cation. Both electrostatic and orbital contributions play a crucial role in stabilising the studied complexes. Based on the NOCV deformation density contributions, the strongest orbital interaction is reached with the bis(indenyl)methyl zirconium cation, which is the only one exhibiting covalent interactions. Especially, the strong contribution of π-back donation (occurring from the occupied orbitals of the zirconium atom to the π* anti-bonding orbital of ethylene) may be a key to understand why this catalyst has a higher polymerisation yield than the other studied catalysts. This work suggests a perspective for predicting values of catalytic activity when theoretically designing novel catalysts of zirconocene type.

The Importance of Electrostatics and Polarization for Noncovalent Interactions: Ionic Hydrogen Bonds vs Ionic Halogen Bonds

A series of 26 hydrogen-bonded complexes between Br⁻ and halogen, oxygen and sulfur hydrogen-bond (HB) donors is investigated at the M06-2X/6–311 + G(2df,2p) level of theory. Analysis using a model in which Br⁻ is replaced by a point charge shows that the interaction energy (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta E}_{Int}$$\end{document}ΔEInt) of the complexes is accurately reproduced by the scaled interaction energy with the point charge (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta E}_{Int}^{PC}$$\end{document}ΔEIntPC).This is demonstrated by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta E}_{Int}=0.86{\Delta E}_{Int}^{PC}$$\end{document}ΔEInt=0.86ΔEIntPC with a correlation coefficient, R² =0.999. The only outlier is (Br-H-Br)⁻, which generally is classified as a strong charge-transfer complex with covalent character rather than a HB complex. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta E}_{Int}^{PC}$$\end{document}ΔEIntPC can be divided rigorously into an electrostatic contribution (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta E}_{ES}^{PC}$$\end{document}ΔEESPC) and a polarization contribution (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta E}_{Pol}^{PC}$$\end{document}ΔEPolPC).Within the set of HB complexes investigated, the former varies between -7.2 and -32.7 kcal mol⁻¹, whereas the latter varies between -1.6 and -11.5 kcal mol⁻¹. Compared to our previous study of halogen-bonded (XB) complexes between Br⁻ and C–Br XB donors, the electrostatic contribution is generally stronger and the polarization contribution is generally weaker in the HB complexes. However, for both types of bonding, the variation in interaction strength can be reproduced accurately without invoking a charge-transfer term. For the Br⁻···HF complex, the importance of charge penetration on the variation of the interaction energy with intermolecular distance is investigated. It is shown that the repulsive character of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta E}_{Int}$$\end{document}ΔEInt at short distances in this complex to a large extent can be attributed to charge penetration.

Assessing the Orbital Contribution in the "Spodium Bond" by Natural Orbital for Chemical Valence-Charge Displacement Analysis

The term "spodium bond" (SpB) has been recently proposed to describe the noncoordinative interaction that can be established between a polarized group 12 metal and a mild Lewis base (LB). Most of the systems showing short metal-donor distances compatible with SpB are characterized by the coexistence of multiple weak interactions, including hydrogen and halogen bonding, making the assessment of real importance of SpB difficult. Here, we show that the relative importance of each contribution can be probed by dissecting the orbital component of the interaction through the extended transition state-natural orbital for chemical valence-charge displacement analysis (ETS-NOCV-CD). The latter gives useful information about relative energies and electrons involved, for model systems ([(thiourea)₂MX₂]···LB; M = Zn, Cd, and Hg; X = Cl and I; and LB = CH₂S, CH₂O, CH₃CN, and CO) and a variety of structures extracted from experimentally characterized adducts, allowing us to demonstrate the lack of a direct correlation between a favorable metal-base distance and the presence of an orbital contribution for the SpB.

New Sensitive and Selective Chemical Sensors for Ni2+ and Cu2+ Ions: Insights into the Sensing Mechanism through DFT Methods

We report the synthesis and theoretical study of two new colorimetric chemosensors with special selectivity and sensitivity to Ni²⁺ and Cu²⁺ ions over other metal cations in the CH₃CN/H₂O solution. Compounds (E)-4-((2-nitrophenyl)diazenyl)-N,N-bis(pyridin-2-ylmethyl)aniline (A) and (E)-4-((3-nitrophenyl)diazenyl)-N,N-bis(pyridin-2-ylmethyl)aniline (B) exhibited a drastic color change from yellow to colorless, which allows the detection of the mentioned metal cations through different techniques. The interaction of sensors with these metal ions induced a new absorption band with a hypsochromic shift to the characteristic signal of the free sensors. A theoretical study via time-dependent density functional theory (TD-DFT) was performed. This method has enabled us to reproduce the hypsochromic shift in the maximum UV-vis absorption band and explain the selective sensing of the ions. For all of the systems studied, the absorption band is characterized by a π → π* transition centered in the ligand. Instead of Ni²⁺ and Cu²⁺ ions, the transition is set toward the σ* molecular orbital with a strong contribution of the 3d_x²-_y² transition (π → 3d_x²-_y²). These absorptions imply a ligand-to-metal charge transfer (LMCT) mechanism that results in the hypsochromic shift in the absorption band of these systems.

Disentanglement of orthogonal hydrogen and halogen bonds via natural orbital for chemical valence: A charge displacement analysis

As known, the electron density of covalently bound halogen atoms is anisotropically distributed, making them potentially able to establish many weak interactions, acting at the same time as halogen bond donors and hydrogen bond acceptors. Indeed, there are many examples in which the halogen and hydrogen bond coexist in the same structure and, if a correct bond analysis is required, their separation is mandatory. Here, the advantages and limitations of coupling the charge displacement analysis with natural orbital for chemical valence method (NOCV-CD) to separately analyze orthogonal weak interactions are shown, for both symmetric and asymmetric adducts. The methodology gives optimal results with intermolecular adducts but, in the presence of an organometallic complex, also intramolecular interactions can be correctly analyzed. Beyond the methodological aspects, it is shown that correctly separate and quantify the interactions can give interesting chemical insights about the systems.

Halogen Bonds in Ligand-Protein Systems: Molecular Orbital Theory for Drug Design

Halogen bonds are highly important in medicinal chemistry as halogenation of drugs, generally, improves both selectivity and efficacy toward protein active sites. However, accurate modeling of halogen bond interactions remains a challenge, since a thorough theoretical investigation of the bonding mechanism, focusing on the realistic complexity of drug–receptor systems, is lacking. Our systematic quantum-chemical study on ligand/peptide-like systems reveals that halogen bonding is driven by the same bonding interactions as hydrogen bonding. Besides the electrostatic and the dispersion interactions, our bonding analyses, based on quantitative Kohn–Sham molecular orbital theory together with energy decomposition analysis, reveal that donor–acceptor interactions and steric repulsion between the occupied orbitals of the halogenated ligand and the protein need to be considered more carefully within the drug design process.

Charge Displacement Analysis-A Tool to Theoretically Characterize the Charge Transfer Contribution of Halogen Bonds

Theoretical bonding analysis is of prime importance for the deep understanding of the various chemical interactions, covalent or not. Among the various methods that have been developed in the last decades, the analysis of the Charge Displacement function (CD) demonstrated to be useful to reveal the charge transfer effects in many contexts, from weak hydrogen bonds, to the characterization of σ hole interactions, as halogen, chalcogen and pnictogen bonding or even in the decomposition of the metal-ligand bond. Quite often, the CD analysis has also been coupled with experimental techniques, in order to give a complete description of the system under study. In this review, we focus on the use of CD analysis on halogen bonded systems, describing the most relevant literature examples about gas phase and condensed phase systems. Chemical insights will be drawn about the nature of halogen bond, its cooperativity and its influence on metal-ligand bond components.

Electrically Tunable Electron Transfer and Binding Interaction between Hydrated Ions and Graphene Oxide

Density functional theory simulations were carried out to study the binding interaction between hydrated Na⁺/Cl^- and graphene oxide (GO) under electric fields. External electric fields can modify the binding interactions of the hydrated ions with GO. The field-dependent binding energy is mainly controlled by the orbital interaction driven by the field-dependent electron transfer, in which miscellaneous electron-transfer routes in the interfaces between hydrated ions and GO surface were disclosed. The electric field is able to influence the electron-transfer degree for each route, thereby creating various electron acceptor-donor coupling interactions. Furthermore, we preliminarily explored the effect of the electric field on the interlayer structure of bilayer GO with NaCl and water confined inside. Electric fields can enlarge the interlayer spacing through tuning of the hydrated ion-GO interactions. Our simulations present a new understanding of hydrated ion-GO interactions in the presence of an electric field, which is expected to be valuable in the electrical modulation of GO nanomaterials.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

The intrinsic strength of the halogen bond: electrostatic and covalent contributions described by coupled cluster theory

36 halogen-bonded complexes YXAR_m (X: F, Cl, Br; Y: donor group; AR_m acceptor group) have been investigated at the CCSD(T)/aug-cc-pVTZ level of theory. Binding energies, geometries, NBO charges, charge transfer, dipole moments, electrostatic potential, electron and energy density distributions, difference density distributions, vibrational frequencies, local stretching and bending force constants, and relative bond strength orders n have been calculated and used to order the halogen bonds according to their intrinsic strength. Halogen bonding is found to arise from electrostatic and strong covalent contributions. It can be strengthened by H-bonding or lone pair delocalization. The covalent character of a halogen bond increases in the way 3c-4e (three-center-four-electron) bonding becomes possible. One can characterize halogen bonds by their percentage of 3c-4e bonding. FCl-phosphine complexes can form relatively strong halogen bonds provided electronegative substituents increase the covalent contributions in form of 3c-4e halogen bonding. Binding energies between 1 and 45 kcal mol^-1 are calculated, which reflects the large variety in halogen bonding.

Electronic structure and bonding of the dinuclear metal M2(CO)10 decacarbonyls: applications of natural orbitals for chemical valence

The nature of the chemical metal-metal bond in M₂(CO)₁₀ (M = Mn, Re, Tc) dinuclear decacarbonyls complexes was investigated for the first time using the natural orbital chemical valence (NOCV) approach combined with the extended transition state (ETS) for energy decomposition analysis (EDA). The optimized geometries carried out at different levels of theory BP86, BLYP, BLYPD and BP86D, showed that the latter method, i.e., BP86D, led to the best agreement with X-ray experimental measurements. The BP86D/TZP results revealed that the computed covalent contribution to the metal-metal bond are 60.5%, 54.1% and 52.0% for Mn-Mn, Re-Re and Tc-Tc, respectively. The computed total interaction energies resulting from attractive terms (ΔE _orb and ΔE _eles), correspond well to experimental predictions, based on bond lengths and energy interaction analysis for the studied complexes.

Molecular properties of metal difluorides and their interactions with CO2 and H2O molecules: a DFT investigation

A computational study of metal difluorides (MF₂; M = Ca to Zn) and their interactions with carbon dioxide and water molecules was performed. The structural parameter values obtained and the results of AIM analysis and energy decomposition analysis indicated that the Ca-F bond is weaker and less ionic than the bonds in the transition metal difluorides. A deformation density plot revealed the stablizing influence of the Jahn-Teller effect in nonlinear MF₂ molecules (e.g., where M= Sc, Ti, Cr). An anaysis of the metal K-edge peaks of the difluorides showed that shifts in the edge energy were due to the combined effects of the ionicity, effective nuclear charge, and the spin state of the metal. The interactions of CO₂ with ScF₂ (Scc3 geometry) and TiF₂ (Tic2 geometry) caused CO₂ to shift from its usual linear geometry to a bent geometry (η²(C=O) binding mode), while it retained its linear geometry (η¹(O) binding mode) when it interacted with the other metal difluorides. Energy decomposition analysis showed that, among the various geometries considered, the Scc3 and Tic2 geometries possessed the highest interaction energies and orbital interaction energies. Heavier transition metal difluorides showed stronger affinities for H₂O, whereas the lighter transition metal (Sc and Ti) difluorides preferred CO₂. Overall, the results of this study suggest that fluorides of lighter transition metals with partially filled d orbitals (e.g., Sc and Ti) could be used for CO₂ capture under moist conditions. Graphical abstract Interaction of metal difluorides with carbon dioxide and water.

Quantitative Assessment of Halogen Bonding Utilizing Vibrational Spectroscopy

A total of 202 halogen-bonded complexes have been studied using a dual-level approach: ωB97XD/aug-cc-pVTZ was used to determine geometries, natural bond order charges, charge transfer, dipole moments, electron and energy density distributions, vibrational frequencies, local stretching force constants, and relative bond strength orders n. The accuracy of these calculations was checked for a subset of complexes at the CCSD(T)/aug-cc-pVTZ level of theory. Apart from this, all binding energies were verified at the CCSD(T) level. A total of 10 different electronic effects have been identified that contribute to halogen bonding and explain the variation in its intrinsic strength. Strong halogen bonds are found for systems with three-center-four-electron (3c-4e) bonding such as chlorine donors in interaction with substituted phosphines. If halogen bonding is supported by hydrogen bonding, genuine 3c-4e bonding can be realized. Perfluorinated diiodobenzenes form relatively strong halogen bonds with alkylamines as they gain stability due to increased electrostatic interactions.

On the Stability of Cis- and Trans-2-Butene Isomers. An Insight Based on the FAMSEC, IQA, and ETS-NOCV Schemes

In the present account, the real space fragment attributed molecular system energy change (FAMSEC) approach, interacting quantum atoms energy decomposition scheme as well as molecular orbitals based the extended transition state scheme coupled with natural orbitals for chemical valence (ETS-NOCV) have been, for the first time, successfully used to delineate factors of importance for stability of the 2-butene conformers (cis-eq, cis-TS, trans-eq, trans-TS). Our results demonstrate that atoms of the controversial H-H contact in cis-eq (i) are involved in attractive interaction dominated by the exchange-correlation term, (ii) are weekly stabilized, (iii) show trends in several descriptors found in other typical H-bonds, and (iv) are part of most stabilized CH-HC fragment (loc-FAMSEC = -3.6 kcal/mol) with most favourably changed intrafragment interactions on trans-eq→cis-eq. Moreover, lower stability of cis-eq vs. trans-eq is linked with the entire HCCH (ethylenic) fragment which destabilized cis-eq (mol-FAMSEC, +3.9 kcal/mol) the most. Although the H-H contact can be linked with smaller, relative to trans-, rotational energy barrier in cis-2-butene, we have proven that to rationalize this phenomenon one must account for changes in interactions between various fragments that constitute the entire molecule. Importantly, we discovered a number of comparable trends in fundamental properties of equivalent molecular fragments on a methyl group rotation; for example, interaction between BP-free H-atoms in trans-eq (involving CH bonds of the methyl and ethylenic units) and BP-linked H-atoms in cis-eq. Clearly, rotational energy barrier cannot be entirely (i) rationalized by the properties of or (ii) attributed to the H-H contact in cis-eq. © 2016 Wiley Periodicals, Inc.

The Halogen Bond

The halogen bond occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a halogen atom in a molecular entity and a nucleophilic region in another, or the same, molecular entity. In this fairly extensive review, after a brief history of the interaction, we will provide the reader with a snapshot of where the research on the halogen bond is now, and, perhaps, where it is going. The specific advantages brought up by a design based on the use of the halogen bond will be demonstrated in quite different fields spanning from material sciences to biomolecular recognition and drug design.

Halogen bonds with coordinative nature: halogen bonding in a S–I+–S iodonium complex

A thorough characterization of unusually strong I⁺–S halogen bonds indicates that halogen bonding is not merely an electrostatic contact but also has a strong coordinative nature, thus expanding the scope of halogen bonding.

Modification of the supramolecular structure of [(thione)IY] (Y = Cl, Br) systems by cooperation of strong halogen bonds and hydrogen bonds

Strong halogen bonds and hydrogen bonds cooperated to control the supramolecular structures of [(thione)IY] (Y = Cl, Br) compounds.

ETS-NOCV description of σ-hole bonding

The ETS-NOCV analysis was applied to describe the σ-hole in a systematic way in a series of halogen compounds, CF3-X (X = I, Br, Cl, F), CH3I, and C(CH3)nH3-n-I (n = 1,2,3), as well as for the example germanium-based systems. GeXH3, X = F, Cl, H. Further, the ETS-NOCV analysis was used to characterize bonding with ammonia for these systems. The results show that the dominating contribution to the deformation density, Δρ 1 , exhibits the negative-value area with a minimum, corresponding to σ-hole. The "size" (spatial extension of negative value) and "depth" (minium value) of the σ-hole varies for different X in CF3-X, and is influenced by the carbon substituents (fluorine atoms, hydrogen atoms, methyl groups). The size and depth of σ-hole decreases in the order: I, Br, Cl, F in CF3-X. In CH3-I and C(CH3)nH3-n-I, compared to CF3-I, introduction of hydrogen atoms and their subsequent replacements by methyl groups lead to the systematic decrease in the σ-hole size and depth. The ETS-NOCV σ-hole picture is consistent with the existence the positive MEP area at the extension of σ-hole generating bond. Finally, the NOCV deformation density contours as well as by the ETS orbital-interaction energy indicate that the σ-hole-based bond with ammonia contains a degree of covalent contribution. In all analyzed systems, it was found that the electrostatic energy is approximately two times larger than the orbital-interaction term, confirming the indisputable role of the electrostatic stabilization in halogen bonding and σ-hole bonding.