Electronic criterion for categorizing the chalcogen and halogen bonds: sulfur–iodine interactions in crystals

Intermolecular Metal-Involving Pnictogen Bonding: The Case of σ-(SbIII)-Hole···dz2[PtII] Interaction

Cocrystallizations of trans-[PtX'2(NCNR2)2] (R2 = Me2, X' = Cl 1a, Br 1b, I 1c; R2 = (CH2)5, X' = I 2c) with SbX3 (X = Cl, Br, I) gave 1:2 cocrystals 1a·2SbCl3, 1b·2SbBr3, 1c·2SbCl3, 1c·2SbBr3, 1c·2SbI3, and 2c·2SbI3. In all six X-ray structures, the association of the molecular coformers is achieved mainly by SbIII···dz2[PtII] metal-involving intermolecular pnictogen bonding. Density functional theory (DFT) calculations (based on experimentally determined geometries) using both gas-phase and solid-state approximations revealed that a σ-(Sb)-hole interacts with an area of negative potential associated with the dz2-orbital of the positively charged platinum(II) sites, thus forming a pnictogen bond whose energy falls in the range between -7.3 and -16.9 kcal/mol.

Dynamics of non-covalent interactions during the P-O bond cleavage reaction by ribonuclease A

In this work, an atomistic-scale investigation of the phosphodiester P-O bond cleavage reaction by the enzyme ribonuclease A was carried out using computer simulation techniques. It is shown that during the reaction the network of non-covalent interactions in the active center of the ribonuclease changes significantly, while the role of these non-covalent interactions is different: coordination of the corresponding groups, electron density transfer, and ligand holding in the active center. It is shown that the process of proton transfer from Asp121 to His119 is the first stage of this reaction; at the same time, the hydrogen bond between the phosphate ligand and the imino group of Arg39 is broken, which, although keeping the ligand in the active center, does not allow the ligand to orient itself more conveniently for subsequent proton transfers. Furthermore, the key step of this reaction occurs: proton transfer with the participation of imidazole rings His12 and His119, in which the guiding role is played by several hydrogen bonds with the participation of Phe120, and the role of an electron density carrier is played by the pnictogen bond between the oxygen of the phosphate ligand and the pyridine-like nitrogen of the imidazole ring His119, which was detected for the first time.

Noncovalent Interactions of Silatranes and Germatranes with the Surface of Stoichiometric and Hydroxylated 2D Silica

Analysis of noncovalent interactions formed by the surface of a 2D silica bilayer and atrane molecules, as precursors of functional layers immobilized on a surface of silicatene, was performed. For this purpose, the systems of substituted silatranes and germatranes adsorbed on silicatene surfaces with different amounts of hydroxyl groups (0, 2, 4, and 30 per cell) were simulated by using quantum chemical modeling with periodic boundary conditions and full-electron basis sets. The observation results for interaction energy changes in the systems "atrane molecule-silicatene surface" with increasing silanol number of the silicatene surface can be used to predict the optimal degree of silicatene hydroxylation in order to control the effective progress of atrane deprotection on activated 2D silica materials. In addition to the typical hydrogen bonds, the ability of atranes to form noncovalent O···O and O···Hal interactions was discovered. In these bonds, the oxygen or halogen atoms of atranes act as electron-donor centers in relation to the silicatene oxygen atoms. The observed weakening of the Ge-O covalent bonds in germatranes, on which further deprotection reaction depends, is more manifested than for the Si-O bonds in adsorbed silatranes.

Electron delocalization in defect-containing graphene and its influence on tetrel bond formation

Using the advanced analyses of electron density and fermionic potential, we show how electron delocalization influences the ability of defect-containing graphene to form tetrel bonds. The Cg atoms of a vacancy defect can produce one nonpolar interaction, alongside a peculiar polar Cg⋯Cg bond. The latter stems from the presence of a localized electron pair on a vacancy defect Cg atom and the local depletion of electron localization on another Cg atom. This interaction is an example of intralayer tetrel bond. In the presence of an absorbed molecule of bisphenol A diglycidyl ether (DGEBA), graphene is able to form incipient tetrel Cg⋯O bonds with an ether group oxygen. In contrast to an epoxy group oxygen, the disposition of the ether oxygen often causes the orientation of electron-rich π-domains of graphene carbon on the weakly expressed electrophilic region of the oxygen. In the case of graphene with a point Si defect, the Si atom can form quite strong Si⋯C interactions with the DGEBA aryl carbons. In contrast to other noncovalent bonds, this interaction significantly alters the electron (de)localization on the Si atom and in the aryl ring. The reliability of the obtained results is enhanced by the use of multiple 2D periodic models with defects located at different positions along the DGEBA skeleton.

Applicability of transferable multipole pseudo-atoms for restoring inner-crystal electronic force density fields. Chemical bonding and binding features in the crystal and dimer of 1,3-bis(2-hydroxyethyl)-6-methyluracil

We considered it timely to test the applicability of transferable multipole pseudo-atoms for restoring inner-crystal electronic force density fields. The procedure was carried out on the crystal of 1,3-bis(2-hydroxyethyl)-6-methyluracil, and some derived properties of the scalar potential and vector force fields were compared with those obtained from the experimental multipole model and from the aspherical pseudo-atom model with parameters fitted to the calculated structure factors. The procedure was shown to accurately replicate the general vector-field behavior, the peculiarities of the quantum potentials and the characteristics of the force-field pseudoatoms, such as charge, shape and volume, as well as to reproduce the relative arrangement of atomic and pseudoatomic zero-flux surfaces along internuclear regions. It was found that, in addition to the quantum-topological atoms, the force-field pseudoatoms are spatially reproduced within a single structural fragment and similar environment. In addition, the classical and nonclassical hydrogen bonds in the uracil derivative crystal, as well as the H...O, N...O and N...C interactions in the free π-stacked dimer of the uracil derivative molecules, were studied using the potential and force fields within the concepts of interatomic charge transfer and electron lone pair donation-acceptance. Remarkably, the nitrogen atoms in the N...O and N...C interactions behave rather like a Lewis base and an electron contributor. At the same time, the hydrogen atom in the H...O interaction, being a Lewis acid, also participates in the interatomic electron transfer by acting as a contributor. Thus, it has been argued that, when describing polar interatomic interactions within orbital-free considerations, it makes more physical sense to identify electronegative (electron occupier) and electropositive (electron contributor) atoms or subatomic fragments rather than nucleophilic and electrophilic sites.

Anion···Anion [AuI4]-···[AuI2]- Complex Trapped in the Solid State by Tetramethylammonium Cations

A discrete π-hole···σ-hole dimer is synthesized and X-ray characterized. It presents a perfect thumbtack geometry where the σ-hole of the linear [AuI₂]^- anion points to the π-hole located above the central Au-atom of the [AuI₄]^- anion. Such discrete π-hole···σ-hole dimers are unprecedented in literature, since all mixed-valence gold(I/III) iodide compounds reported to date form infinite ···([AuI₄]^-···[AuI₂]^-) _n ·· chains in the solid state. If an excess of iodine is used for the synthesis, triiodide [I₃]^- ions are partially incorporated into the [AuI₂]^- sites, forming infinite chains. The nature of the anion···anion interaction has been studied considering two possibilities: (i) a π-hole coinage bond or (ii) σ-hole halogen bond using high-level density functional theory calculations, the quantum theory of atoms in molecules, and the noncovalent interaction plot index.

Halogen Bonding Involving Gold Nucleophiles in Different Oxidation States

A single-crystal X-ray diffraction (XRD) study of diaryliodonium tetrachloroaurates (or, in the recent terminology, tetrachloridoaurates), [(p-XC₆H₄)₂I][AuCl₄] (X = Cl, 1; Br, 2), was performed for 1 (the structure is denoted as 1a to show similarity with the isomorphic structure 2a) and two polymorphs─2a (obtained from MeOH) and 2b (from 1,2-C₂H₄Cl₂). Examination of the XRD data for these three structures revealed 2-center C-X···Au^III (X = Cl and Br) and 3-center bifurcated C-Br···(Cl-Au) halogen bonding (abbreviated as XB) between the p-Cl or p-Br atoms of the diaryliodonium cations and the gold(III) atom of [AuCl₄]^-. The noncovalent nature of Au^III-involving interactions, the nucleophilicity of the gold(III) atoms, and the electrophilic role of p-X atoms of the diaryliodonium cations in the XBs were studied by a set of complementary computational methods. Combined experimental and theoretical studies allowed the recognition of the d-nucleophilicity of the [d⁸Au^III] atom which, regardless of its rather substantial formal 3+ charge, can function as a d-nucleophilic partner of XB. This conclusion was also supported by theoretical calculations performed for the structures' refcodes BINXOM and ICSD 62511; the obtained data verified the nucleophilicity of Au^III toward a K⁺ ions or a σ-(Cl)-hole, respectively. All our results, together with consideration of relevant literature, indicate that gold atoms in the three oxidation states (0, I, and even III) exhibit nucleophilicity in XBs.

Crystallographic and Theoretical Study of Osme Bonds in Nitrido-Osmium(VI) Complexes

Osme bonds have been recently defined as the attractive interaction between an element of group 8 acting as an electrophile and any atom or group of atoms acting as a nucleophile. To date, the known examples of osme bonds in X-ray structures involve mostly the highly reactive OsO4 and amines and amine oxides. In this work, evidence supporting the existence of osme bonds in osmium(VI) derivatives is reported. In particular, nitrido-osmium(VI) complexes that present square-pyramidal geometries are well disposed to participate in osme bonds opposite to the Os≡N bond. By using a combination of experimental and theoretical results, the existence and importance of this new class of σ-hole interactions is demonstrated in the solid state of several nitrido-osmium(VI) derivatives.

Can We Merge the Weak and Strong Tetrel Bonds? Electronic Features of Tetrahedral Molecules Interacted with Halide Anions

Using the orbital-free quantum crystallography approach, we have disclosed the quantitative trends in electronic features for bonds of different strengths formed by tetrel (Tt) atoms in stable molecular complexes consisting of electrically neutral tetrahedral molecules and halide anions. We have revealed the role of the electrostatic and exchange-correlation components of the total one-electron static potential that are determined by the equilibrium atomic structure and by kinetic Pauli potential, which reflects the spin-dependent electron motion features of the weak and strong bonds. The gap between the extreme positions in the electrostatic and total static potentials along the line linking the Tt atom and halide anion is wide for weak bonds and narrow for strong ones. It is in very good agreement with the number of minima in the Pauli potential between the bounded atoms. This gap exponentially correlates with the exchange-correlation potential in various series with a fixed nucleophilic fragment. A criterion for categorizing the noncovalent tetrel bonds (TtB) based on the potential features is suggested.

The Distance between Minima of Electron Density and Electrostatic Potential as a Measure of Halogen Bond Strength

In this study, we present results of a detailed topological analysis of electron density (ED) of 145 halogen-bonded complexes formed by various fluorine-, chlorine-, bromine-, and iodine-containing compounds with trimethylphosphine oxide, Me3PO. To characterize the halogen bond (XB) strength, we used the complexation enthalpy, the interatomic distance between oxygen and halogen, as well as the typical set of electron density properties at the bond critical points calculated at B3LYP/jorge-ATZP level of theory. We show for the first time that it is possible to predict the XB strength based on the distance between the minima of ED and molecular electrostatic potential (ESP) along the XB path. The gap between ED and ESP minima exponentially depends on local electronic kinetic energy density at the bond critical point and tends to be a common limiting value for the strongest halogen bond.

Non-covalent interactions in the glutathione peroxidase active center and their influence on the enzyme activity

In this computational work, the structure of the active center of the enzyme glutathione peroxidase (in three forms -SeH, -SeOH and -Se(O)OH) and the non-covalent interactions in it were investigated using modern quantum chemistry methods. The non-covalent interactions are described in detail. The presence of σ-hole interactions (chalcogen, tetrel and pnictogen bonds) formed mostly by a selenium atom as an electrophile in the glutathione peroxidase active center is confirmed for the first time. It is shown that a number of non-covalent interactions stabilize intermediates along the catalytic cycle and that modelling of the whole enzyme active center is necessary for accurate predictions of thermodynamic parameters, in particular, activation barriers.

Synthesis of Polysubstituted Ferrocenesulfoxides

The purpose of the study is to design synthetic methodologies, especially directed deprotometalation using polar organometallic reagents, to access polysubstituted ferrocenesulfoxides. From enantiopure 2-substituted (SiMe₃, PPh₂) S-tert-butylferrocenesulfoxides, a third substituent was first introduced at the 5 position (SiMe₃, I, D, C(OH)Ph₂, Me, PPh₂, CH₂NMe₂, F) and removal of the trimethylsilyl group then afforded 2-substituted ferrocenesulfoxides unreachable otherwise. Attempts to apply the “halogen dance” reaction to the ferrocenesulfoxide series led to unexpected results although rationalized in light of calculated pK_a values. Further functionalizations were also possible. Thus, new enantiopure, planar chiral di- and trisubstituted ferrocenes have been obtained, in addition to several original 2-substituted, 2,3- and 2,5-disubstituted, 2,3,5-trisubstituted and even 2,3,4,5-tetrasubstituted ferrocenesulfoxides, also enantiopure.

Inorganic–Organic {dz2-MIIS4}···π-Hole Stacking in Reverse Sandwich Structures. The Case of Cocrystals of Group 10 Metal Dithiocarbamates with Electron-deficient Arenes

Cocrystallization of the dithiocarbamate complexes [M(S2CNEt2)2] (M = Ni 1, Pd 2, Pt 3) and X-substituted perfluoroarenes (X = I, Br; 1,2-dibromoperfluorobenzene FBrB and 1,2-diiodoperfluorobenzene FIB) gives isomorphous cocrystals of...

Chameleonic Metal-bound Isocyanides: π-Donating CuI-center Imparts a Nucleophilicity to the Isocyanide Carbon toward Halogen Bonding

In the structures of the isostructural cocrystals [CuI3(CNXyl)3]·CHX3 (X = Br, I), two adjacent CuI-bound isocyanide groups, whose carbon lone pairs are blocked by the ligation, exhibit nucleophilic properties induced...

Disulfide metathesis via sulfur⋯iodine interaction and photoswitchability

The idea of constitutional dynamic chemistry (CDC) and dynamic combinatorial chemistry (DCC) is widespread in the literature using the chemistry of disulfides. The synthesis of unsymmetrical diaryl disulfides is challenging due to the presence of a weak S-S bond. We report herein the synthesis of unsymmetrical diaryl disulfides from two symmetrical disulfides via a cross-metathesis reaction which was controlled by a weak sulfur⋯iodine (S⋯I) interaction. The unsymmetrical disulfides were stable in acetonitrile solution in the presence of N-iodosuccinimide (NIS), and found to be reversibly photoswitchable to the symmetrical disulfides under visible light irradiation.

Studies of Nature of Uncommon Bifurcated I-I···(I-M) Metal-Involving Noncovalent Interaction in Palladium(II) and Platinum(II) Isocyanide Cocrystals

Two isostructural trans-[MI2(CNXyl)2]·I2 (M = Pd or Pt; CNXyl = 2,6-dimethylphenyl isocyanide) metallopolymeric cocrystals containing uncommon bifurcated iodine···(metal-iodide) contact were obtained. In addition to classical halogen bonding, single-crystal X-ray diffraction analysis revealed a rare type of metal-involved stabilizing contact in both cocrystals. The nature of the noncovalent contact was studied computationally (via DFT, electrostatic surface potential, electron localization function, quantum theory of atoms in molecules, and noncovalent interactions plot methods). Studies confirmed that the I···I halogen bond is the strongest noncovalent interaction in the systems, followed by weaker I···M interaction. The electrophilic and nucleophilic nature of atoms participating in I···M interaction was studied with ED/ESP minima analysis. In trans-[PtI2(CNXyl)2]·I2 cocrystal, Pt atoms act as weak nucleophiles in I···Pt interaction. In the case of trans-[PdI2(CNXyl)2]·I2 cocrystal, electrophilic/nucleophilic roles of Pd and I are not clear, and thus the quasimetallophilic nature of the I···Pd interaction was suggested.

Competition between Inter and Intramolecular Tetrel Bonds: Theoretical Studies Complemented by CSD Survey

Crystal structures document the ability of a TF₃ group (T=Si, Ge, Sn, Pb) situated on a naphthalene system to engage in an intramolecular tetrel bond (TB) with an amino group on the adjoining ring. Ab initio calculations evaluate the strength of this bond and evaluate whether it can influence the ability of the T atom to engage in a second, intermolecular TB with another nucleophile. A very strong CN^- anionic base can approach the T either along the extension of a T-C or T-F bond and form a strong TB with an interaction energy approaching 100 kcal/mol, although this bond is weakened a bit by the presence of the internal T⋅⋅⋅N bond. The much less potent NCH base engages in a correspondingly longer and weaker TB, less than 10 kcal/mol. Such an intermolecular TB is weakened by the presence of the internal TB, to the point that it only occurs for the two heavier tetrel atoms Sn and Pb.

The explicit role of electron exchange in the hydrogen bonded molecular complexes

We applied a set of advanced bonding descriptors to establish the hidden electron density features and binding energy characteristics of intermolecular DH∙∙∙A hydrogen bonds (OH∙∙∙O, NH∙∙∙O and SH∙∙∙O) in 150 isolated and solvated molecular complexes. The exchange-correlation and Pauli potentials as well as corresponding local one-electron forces allowed us to explicitly ascertain how electron exchange defines the bonding picture in the proximity of the H-bond critical point. The electron density features of DH∙∙∙A interaction are governed by alterations in the electron localization in the H-bond region displaying itself in the exchange hole. At that, they do not depend on the variations in the exchange hole mobility. The electrostatic interaction mainly defines the energy of H-bonds of different types, whereas the strengthening/weakening of H-bonds in complexes with varying substituents depends on the barrier height of the exchange potential near the bond critical point. Energy variations between H-bonds in isolated and solvated systems are also caused the electron exchange peculiarities as follows from the corresponding potential and the interacting quantum atom analyses complemented by electron delocalization index calculations. Our approach is based on the bonding descriptors associated with the characteristics of the observable electron density and can be recommended for in-depth studies of non-covalent bonding.

Orbital-Free Quantum Crystallographic View on Noncovalent Bonding: Insights into Hydrogen Bonds, π⋅⋅⋅π and Reverse Electron Lone Pair⋅⋅⋅π Interactions

A detailed analysis of a complete set of the local potentials that appear in the Euler equation for electron density is carried out for noncovalent interactions in the crystal of a uracil derivative using experimental X-ray charge density. The interplay between the quantum theory of atoms in molecules and crystals and the local potentials and corresponding inner-crystal electronic forces of electrostatic and kinetic origin is explored. Partitioning of crystal space into atomic basins and atomic-like potential basins led us to the definite description of interatomic interaction and charge transfer. Novel physically grounded bonding descriptors derived within the orbital-free quantum crystallography provided the detailed examination of π-stacking and intricate C=O⋅⋅⋅π interactions and nonclassical hydrogen bonds present in the crystal. The donor-acceptor character of these interactions is revealed by analysis of Pauli and von Weizsäcker potentials together with well-known functions, e. g., deformation electron density and electron localization function. In this way, our analysis throws light on aspects of these closed-shell interactions hitherto hidden from the description.

Bifurcated Halogen Bonding Involving Two Rhodium(I) Centers as an Integrated σ-Hole Acceptor

The complexes [RhX(COD)]2 (X = Cl, Br; COD = 1,5-cyclooctadiene) form cocrystals with σ-hole iodine donors. X-ray diffraction studies and extensive theoretical considerations indicate that the d z 2-orbitals of two positively charged rhodium(I) centers provide sufficient nucleophilicity to form a three-center halogen bond (XB) with the σ-hole donors. The two metal centers function as an integrated XB acceptor, providing assembly via a metal-involving XB.

Synthesis and supramolecular organization of the iodide and triiodides of a polycyclic adamantane-based diammonium cation: the effects of hydrogen bonds and weak I⋯I interactions

Adamantane-like divalent building blocks and iodide or polyiodide anions combine into supramolecular architectures with the help of various noncovalent forces ranging from strong hydrogen bonds to secondary and weak I⋯I interactions.

The X-ray constrained wavefunction of the [Mn(CO)4{(C6H5)2P-S-C(Br2)-P(C6H5)2}]Br complex: a theoretical and experimental study of dihalogen bonds and other noncovalent interactions

The synthesis and X-ray structure determination of the [Mn(CO)₄{(C₆H₅)₂P-S-C(Br₂)-P(C₆H₅)₂}]Br complex (1) are described. The C-Br...Br dihalogen bond present in 1 has been characterized by means of topological studies of the electron density. Both the quantum theory of atoms in molecules and the electron localization function approaches have been applied to several theoretically calculated wavefunctions as well as to an X-ray constrained wavefunction. In addition, a number of theoretical techniques, such as the source function, the reduced density gradient method and the interacting quantum atoms approach, among others, have been used to analyse the dihalogen bond as well as several intramolecular interactions of the type C-H...Br-C which have also been detected in 1. The results show clearly that while bonding in the latter interactions are dominated by electrostatic components, the former has a high degree of covalency.

The shortest chalcogen...halogen contacts in molecular crystals

The survey of the shortest contacts in structures deposited in the Cambridge Structural Database shows that chalcogen...halogen, halogen...halogen and chalcogen...chalcogen interactions can compete as cohesion forces in molecular crystals. The smallest parameter δ (defined as the interatomic distance minus the sum of relevant van der Waals radii) for Ch...X contacts between chalcogens (Ch: S, Se) and halogens (X: F, Cl, Br, I) is present only in 0.86% out of 30 766 deposited structures containing these atoms. Thus, in less than 1% of these structures can the Ch...X forces be considered as the main type of cohesion forces responsible for the molecular arrangement. Among the 263 structures with the shortest Ch...X contact, there are four crystals where no contacts shorter than the sums of van der Waals radii are present (so-called loose crystals). The smallest δ criterion has been used for distinguishing between the bonding (covalent bond) and non-bonding contacts and for validating the structural models of crystals.

Noncovalent Bonds, Spectral and Thermal Properties of Substituted Thiazolo[2,3-b][1,3]thiazinium Triiodides

The interrelation between noncovalent bonds and physicochemical properties is in the spotlight due to the practical aspects in the field of crystalline material design. Such study requires a number of similar substances in order to reveal the effect of structural features on observed properties. For this reason, we analyzed a series of three substituted thiazolo[2,3-b][1,3]thiazinium triiodides synthesized by an iodocyclization reaction. They have been characterized with the use of X-ray diffraction, Raman spectroscopy, and thermal analysis. Various types of noncovalent interactions have been considered, and an S…I chalcogen bond type has been confirmed using the electronic criterion based on the calculated electron density and electrostatic potential. The involvement of triiodide anions in the I…I halogen and S…I chalcogen bonding is reflected in the Raman spectroscopic properties of the I–I bonds: identical bond lengths demonstrate different wave numbers of symmetric triiodide vibration and different values of electron density at bond critical points. Chalcogen and halogen bonds formed by the terminal iodine atom of triiodide anion and numerous cation…cation pairwise interactions can serve as one of the reasons for increased thermal stability and retention of iodine in the melt under heating.

Chalcogen bonding of two ligands to hypervalent YF4 (Y = S, Se, Te, Po)

The ability of two NH₃ ligands to engage in simultaneous chalcogen bonds to a hypervalent YF₄ molecule, with Y = S, Se, Te, Po, is assessed via quantum calculations. The complex can take on one of two different geometries. The cis structure places the two ligands adjacent to one another in a pseudo-octahedral geometry, held there by a pair of σ-hole chalcogen bonds. The bases can also lie nearly opposite one another, in a distorted octahedron containing one π-hole and one strained σ-hole bond. The cis geometry is favored for Y = S, while Te, and Po tend toward the trans structure; they are nearly equally stable for Se. In either case, the binding energy rises rapidly with the size of the Y atom, exceeding 30 kcal mol^-1 for PoF₄.

Identification of the Tetrel Bonds between Halide Anions and Carbon Atom of Methyl Groups Using Electronic Criterion

The consideration of the disposition of minima of electron density and electrostatic potential along the line between non-covalently bound atoms in systems with Hal-···CH₃⁻Y (Hal- = Cl, Br; Y = N, O) fragments allowed to prove that the carbon atom in methyl group serves as an electrophilic site provider. These interactions between halide anion and carbon in methyl group can be categorized as the typical tetrel bonds. Statistics of geometrical parameters for such tetrel bonds in CSD is analyzed. It is established that the binding energy in molecular complexes with tetrel bonds correlate with the potential acting on an electron in molecule (PAEM). The PAEM barriers for tetrel bonds show a similar behavior for both sets of complexes with Br- and Cl- electron donors.