ETS-NOCV description of σ-hole bonding

Unveiling the Inverse Sandwich Complexes of XeO3: A Computational Exploration

Our study introduces the design of inverse sandwich (iSw) complexes incorporating a noble gas compound: xenon trioxide (XeO3). Through comprehensive computational analyses, we have investigated the critical factors influencing their stability by employing a variety of state-of-the-art computational tools. We demonstrated that the coordination number of xenon in the iSw complex of XeO3 with 18-crown-6 is influenced by the presence of a rare, weakly stabilizing Xe···Xe interaction between the XeO3 molecules. Furthermore, we observed that the stability of iSw complexes of 1,3,5-triphenylbenzene (TPB) and its derivatives is not solely attributed to aerogen bonding, but also involves contributions from C-H···O interactions and back-donation from the lone pair of Xe to the antibonding C-C orbitals of TPB. Additionally, the significant contributions from orbital interactions and dispersion interactions in the TPB derivatives highlight the multifaceted amphoteric properties of XeO3 and reveal that the iSw complexes of TPB and derivatives are not predominantly governed by electrostatic interactions, contrary to conventional belief.

Static and Dynamical Quantum Studies of CX3-AlX2 and CSiX3-BX2 (X = F, Cl, Br) Complexes with Hydrocyanic Acid: Unusual Behavior of Strong π-Hole at Triel Center

The set of TX₃-TrX₂ (T = C, Si, Ge; Tr = B, Al, Ga; X = F, Cl, Br) molecules offers a rather unique opportunity to study both σ-hole and π-hole dimerization on the tetrel and triel ends, respectively. According to the molecular electrostatic potential (MEP) distribution, the π-hole extrema (acidic sites) were more intense than their σ-hole counterparts. The molecules owning the most (CX₃-AlX₂) and least (SiX₃-BX₂) intense π-holes were chosen to evaluate their capacities to attract one and two HCN molecules (Lewis bases). We discovered that the energetic characteristics of π-hole dimers severely conflict with the monomers MEP pattern since the weakest π-hole monomer forms a dimer characterized by interaction energy compared to those created by the monomers with noticeably greater power in the π-hole region. This outcome is due to the deformation of the weakest π-hole donor. Furthermore, the MEP analysis for monomers in the geometry of respective dimers revealed a "residual π-hole" site that was able to drive second ligand attachment, giving rise to the two "unusual trimers" examined further by the NCI and QTAIM analyses. Apart from them, the π-hole/π-hole and σ-hole/π-hole trimers have also been obtained throughout this study and described using energetic and geometric parameters. The SAPT approach revealed details of the bonding in one of the "unusual trimers". Finally, Born-Oppenheimer Molecular Dynamics (BOMD) simulations were carried out to investigate the time evolution of the interatomic distances of the studied complexes as well as their stability.

Probing the Directionality of S···O/N Chalcogen Bond and Its Interplay with Weak C-H···O/N/S Hydrogen Bond Using Molecular Electrostatic Potential

The directionality of the chalcogen bond (Ch-bond) formed by S and its interplay with other weak interactions have important chemical and biological implications. Here, dimers made of CH₃-S-X and O/N containing nucleophiles are studied and found to be stabilized by coexisting S···O/N and C-H···O/N interactions. Based on experimentally accessible electron density and molecular electrostatic potentials (MESPs), we showed that reciprocity between S···O/N and C-H···O/N interactions in the stability of cumulative molecular interaction (ΔE) was dependent on the strength of the σ-hole on S (V_s,max). Direct correlation between ΔE of dimers with V_s,max of S supports the electrostatic nature of the Ch-bond. Such interplay of the Ch-bond is necessary for its directionality in complex nucleophiles (carbonyl groups) with multiple electron-rich centers, which is explained using MESP. A correlation between the MESP minima in the π-region and the strength of the S-π interaction explains the directional selectivity of the Ch-bond.

Bonding analysis of the C2 precursor Me3E–C2–I(Ph)FBF3 (E = C, Si, Ge)

A series of possible precursors for generating C2 with the general formula Me3E–C2–I(Ph)FBF3 [E = C (1), Si (2), and Ge (3)] has been theoretically investigated using quantum chemical calculations. The equilibrium geometries of all species show a linear E–C2–I+ backbone. The inspection of the electronic structure of the Me3E–C2 bond by energy decomposition analysis coupled with the natural orbital for chemical valence (EDA-NOCV) method suggests a combination of electron sharing C–C σ-bond and v weak π-dative bond between Me3C and C2 fragments in the doublet state for species 1 (E = C). For species 2 (Si) and 3 (Ge), the analysis reveals σ-dative Me3E–C2 bonds (E = Si, Ge; Me3E←C2) resulting from the interaction of singly charged (Me3E)+ and (C2–IPh(BF4))− fragments in their singlet states. The C2–I bond is diagnosed as an electron sharing σ-bond in all three species, 1, 2 and 3.

π-Hole Tetrel Bonds—Lewis Acid Properties of Metallylenes

The MP2/aug-cc-pVTZ calculations were performed on the dihalometallylenes to indicate their Lewis acid and Lewis base sites. The results of the Cambridge Structural Database search show corresponding and related crystal structures where the tetrel center often possesses the configuration of a trigonal bipyramid or octahedron. The calculations were also carried out on dimers of dichlorogermylene and dibromogermylene and on complexes of these germylenes with one and two 1,4-dioxide molecules. The Ge⋯Cl, Ge⋯Br, and Ge⋯O interactions are analyzed. The Ge⋯O interactions in the above mentioned germylene complexes may be classified as the π-hole tetrel bonds. The MP2 calculations are supported by the results of the Quantum Theory of Atoms in Molecules (QTAIM) and the Natural Bond Orbital (NBO) approaches.

A Combined Experimental/Quantum-Chemical Study of Tetrel, Pnictogen, and Chalcogen Bonds of Linear Triatomic Molecules

Linear triatomic molecules (CO₂, N₂O, and OCS) are scrutinized for their propensity to form perpendicular tetrel (CO₂ and OCS) or pnictogen (N₂O) bonds with Lewis bases (dimethyl ether and trimethyl amine) as compared with their tendency to form end-on chalcogen bonds. Comparison of the IR spectra of the complexes with the corresponding monomers in cryogenic solutions in liquid argon enables to determine the stoichiometry and the nature of the complexes. In the present cases, perpendicular tetrel and pnictogen 1:1 complexes are identified mainly on the basis of the lifting of the degenerate ν 2 bending mode with the appearance of both a blue and a red shift. Van ′t Hoff plots of equilibrium constants as a function of temperature lead to complexation enthalpies that, when converted to complexation energies, form the first series of experimental complexation energies on sp¹ tetrel bonds in the literature, directly comparable to quantum-chemically obtained values. Their order of magnitude corresponds with what can be expected on the basis of experimental work on halogen and chalcogen bonds and previous computational work on tetrel bonds. Both the order of magnitude and sequence are in fair agreement with both CCSD(T) and DFA calculations, certainly when taking into account the small differences in complexation energies of the different complexes (often not more than a few kJ mol⁻¹) and the experimental error. It should, however, be noted that the OCS chalcogen complexes are not identified experimentally, most probably owing to entropic effects. For a given Lewis base, the stability sequence of the complexes is first successfully interpreted via a classical electrostatic quadrupole–dipole moment model, highlighting the importance of the magnitude and sign of the quadrupole moment of the Lewis acid. This approach is validated by a subsequent analysis of the molecular electrostatic potential, scrutinizing the σ and π holes, as well as the evolution in preference for chalcogen versus tetrel bonds when passing to “higher” chalcogens in agreement with the evolution of the quadrupole moment. The energy decomposition analysis gives further support to the importance/dominance of electrostatic effects, as it turns out to be the largest attractive term in all cases considered, followed by the orbital interaction and the dispersion term. The natural orbitals for chemical valence highlight the sequence of charge transfer in the orbital interaction term, which is dominated by an electron-donating effect of the N or O lone-pair(s) of the base to the central atom of the triatomics, with its value being lower than in the case of comparable halogen bonding situations. The effect is appreciably larger for TMA, in line with its much higher basicity than DME, explaining the comparable complexation energies for DME and TMA despite the much larger dipole moment for DME.

Yet another perspective on hole interactions

Hole interactions are known by different names depending on the key atom of the bond (halogen bond, chalcogen bond, hydrogen bond, etc.), and the geometry of the interaction (σ if in line, π if perpendicular to the Lewis acid plane). However, its origin starts with the creation of a Lewis acid by an underlying covalent bond, which forms an electrostatic depletion and a virtual antibonding orbital, which can create non-covalent interactions with Lewis bases. In this (maybe subjective) perspective, we will claim that hole interactions must be defined via the molecular orbital origin of the molecule. Under this premise we can better explore the richness of such bonding patterns. For that, we will study old, recent and new systems, trying to pinpoint some misinterpretations that are often associated with them. We will use as exemplars the triel bonds, a couple of metal complexes, a discussion on convergent σ-holes, and many cases of anti-electrostatic hole interactions.

Tuning the activity of known drugs via the introduction of halogen atoms, a case study of SERT ligands - Fluoxetine and fluvoxamine

The selective serotonin reuptake inhibitors (SSRIs), acting at the serotonin transporter (SERT), are one of the most widely prescribed antidepressant medications. All five approved SSRIs possess either fluorine or chlorine atoms, and there is a limited number of reports describing their analogs with heavier halogens, i.e., bromine and iodine. To elucidate the role of halogen atoms in the binding of SSRIs to SERT, we designed a series of 22 fluoxetine and fluvoxamine analogs substituted with fluorine, chlorine, bromine, and iodine atoms, differently arranged on the phenyl ring. The obtained biological activity data, supported by a thorough in silico binding mode analysis, allowed the identification of two partners for halogen bond interactions: the backbone carbonyl oxygen atoms of E493 and T497. Additionally, compounds with heavier halogen atoms were found to bind with the SERT via a distinctly different binding mode, a result not presented elsewhere. The subsequent analysis of the prepared XSAR sets showed that E493 and T497 participated in the largest number of formed halogen bonds. The XSAR library analysis led to the synthesis of two of the most active compounds (3,4-diCl-fluoxetine 42, SERT K_i = 5 nM and 3,4-diCl-fluvoxamine 46, SERT K_i = 9 nM, fluoxetine SERT K_i = 31 nM, fluvoxamine SERT K_i = 458 nM). We present an example of the successful use of a rational methodology to analyze binding and design more active compounds by halogen atom introduction. 'XSAR library analysis', a new tool in medicinal chemistry, was instrumental in identifying optimal halogen atom substitution.

Catalytic activity of a new Ru(ii) complex for the hydrogen transfer reaction of acetophenone and N-benzylideneaniline: synthesis, characterization and relativistic DFT approaches

The synthesis and characterization of a new ruthenium(ii) complex containing a hemilabile P^N-ligand are reported.

Halogen bonding in differently charged complexes: basic profile, essential interaction terms and intrinsic σ-hole

Studies on halogen bonds (XB) between organohalogens and their acceptors in crystal structures revealed that the XB donor and acceptor could be differently charged, making it difficult to understand the nature of the interaction, especially the negatively charged donor's electrophilicity and positively charged acceptor's nucleophilicity. In this paper, 9 XB systems mimicking all possibly charged halogen bonding interactions were designed and explored computationally. The results revealed that all XBs could be stable, with binding energies after removing background interaction as strong as -1.2, -3.4, and -8.3 kcal mol-1 for Cl, Br, and I involved XBs respectively. Orbital and dispersion interactions are found to be always attractive while unidirectional intermolecular electron transfer from a XB acceptor to a XB donor occurs in all XB complexes. These observations could be attributed to the intrinsic σ-hole of the XB donor and the intrinsic electronic properties of the XB acceptor regardless of their charge states. Intramolecular charge redistribution inside both the donor and the acceptor is found to be system-dependent but always leads to a more stable XB. Accordingly, this study demonstrates that the orbital-based origin of halogen bonds could successfully interpret the complicated behaviour of differently charged XB complexes, while electrostatic interaction may dramatically change the overall bonding strength. The results should further promote the application of halogens in all related areas.

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.

Tuning the molecular antenna effect using donor and acceptor substituents on the optical properties of the [(C5F5)2ThMCp2]2+ and [(C5F5)2ThMCpL2]+ complexes, where M = Fe, Ru and Os and L = CO and C5H5N

In the present work, a theoretical methodology based on DFT is used to establish the effects of the electron rearrangements on the optical properties in a series of Th-Transtion metal complexes.

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.

Exploiting directional long range secondary forces for regulating electrostatics-dominated noncovalent interactions

It has been well established that long range secondary electrostatic interactions (SEIs) have a significant effect on the stability of supramolecular complexes. However, general rules for exploiting SEIs in the rational design of diverse supramolecular complexes have been difficult to obtain. In this work, we outline a quantum chemical approach for understanding the strength of electrostatic interactions. This approach is seen to provide excellent correlation between the electrostatic force and the binding energy between two partners in hydrogen-bonded complexes, as well as that between two ions in ion-pair complexes. Furthermore, we illustrate how the understanding of the binding allows for the rational design of new complexes where the association constant between the two partners can be increased or decreased, as desired, by several orders of magnitude. Hence, the current work showcases a general, simple and powerful method of understanding and exploiting long range secondary electrostatic interactions.

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.

Asymmetric bifurcated halogen bonds

Halogen bonding (XB) is being extensively explored for its potential use in advanced materials and drug design. Despite significant progress in describing this interaction by theoretical and experimental methods, the chemical nature remains somewhat elusive, and it seems to vary with the selected system. In this work we present a detailed DFT analysis of three-center asymmetric halogen bond (XB) formed between dihalogen molecules and variously 4-substituted 1,2-dimethoxybenzene. The energy decomposition, orbital, and electron density analyses suggest that the contribution of electrostatic stabilization is comparable with that of non-electrostatic factors. Both terms increase parallel with increasing negative charge of the electron donor molecule in our model systems. Depending on the orientation of the dihalogen molecules, this bifurcated interaction may be classified as 'σ-hole - lone pair' or 'σ-hole - π' halogen bonds. Arrangement of the XB investigated here deviates significantly from a recent IUPAC definition of XB and, in analogy to the hydrogen bonding, the term bifurcated halogen bond (BXB) seems to be appropriate for this type of interaction.

Metal-interacted histidine dimer: an ETS-NOCV and XANES study

We have analyzed the metal coordination in a histidine dimer, hydrated with a water molecule, based on the extended transition state scheme with the theory of natural orbitals for chemical valence (ETS-NOCV).

Synergistic effects of halogen bond and π–π interactions in thiophene-based building blocks

The use of a thiophene-based building block (Pyr-T) to achieve a supramolecular assembly driven by halogen bonding (XB) is described. X-ray crystallography, NMR, Raman spectroscopy, and computational chemistry afford evidence of strong XB and π–π stacking.

On the origin of internal rotation in ammonia borane

The internal rotation in ammonia borane (AB) was studied on the basis of natural orbitals for chemical valence (NOCV) and eigenvectors for Pauli repulsion (NOPR). We found that the total hyperconjugation stabilization (ca. 5 kcal mol⁻¹), based on the charge transfer from the occupied σ (B–H) orbitals into the empty σ*(N–H), slightly favors the staggered conformation over the eclipsed one; however, the barrier to internal rotation in ammonia borane can be understood predominantly in a ‘classical’ way, as originating from the steric (Pauli) repulsion contributions (of the kinetic origin) that act solely between N–H and B–H bonds. Repulsion between the lone pair of ammonia and the adjacent B–H bonds was found to be dominant in absolute terms; however, it does not determine the rotational barrier. Similar conclusions on the role of CH↔HC repulsion appeared to be valid for isoelectronic ethane.

Figure

Halogen bonding interaction of chloromethane with several nitrogen donating molecules: addressing the nature of the chlorine surface σ-hole

The small and positive σ-hole on the outermost extension of the C–Cl covalent bond in CH₃Cl can adequately form halogen bonded interactions with the RN series of thirteen other molecules, R = FC, ClC, BrC, CH₃C, HOC, HSC, PCC, P, CCl₃C, SiH₃C, NCC, CNC, and NaC.