Aerogen Bonding Interaction: A New Supramolecular Force?

Intramolecular Noncovalent Carbon Bonding Interaction Stabilizes the cis Conformation in Acylhydrazones

Noncovalent carbon bonding, a recently explored σ-hole interaction, was hitherto supposed to be a weak and structure-guided interaction. Here, its role in the intramolecular stabilization of the cis conformation of the amide moiety in acylhydrazones is described. The calculations reveal an electron donation from the lone pair of the nitrogen atom to the empty antibonding C-N orbital [LP(N)→BD*(C-N)] with a concomitant stabilization energy of E⁽²⁾ =1.2 kcal mol^-1 .

Regium bonds between Mn clusters (M = Cu, Ag, Au and n = 2-6) and nucleophiles NH3 and HCN

The most stable geometries of the coinage metal (or regium) atom (Cu, Ag, Au) clusters Mn for n up to 6 are all planar, and adopt the lowest possible spin multiplicity. Clusters with even numbers of M atoms are thus singlets, while those with odd n are open-shell doublets. Examination of the molecular electrostatic potential (MEP) of each cluster provides strong indications of the most likely site of attack by an approaching nucleophile, generally one of two positions. A nucleophile (NH3 or HCN) most favorably approaches one particular M atom of each cluster, rather than a bond midpoint or face. In the closed-shell clusters, the interaction energies are highly dependent upon the intensity of the MEP, but this correlation fades for the open-shell systems studied in this work. The strength of the interaction is also closely related to the basicity of the nucleophile. Regium bond energies can be more than 30 kcal mol-1 and tend to follow the Au > Cu > Ag order. These interaction energies are in large part derived from Coulombic attraction, with a smaller orbital interaction contribution.

Boron triel bonding: a weak electrostatic interaction lacking electron-density descriptors

In an effort to describe π-hole interactions, we undertook accurate high-resolution X-ray diffraction analyses of single crystals of 1,4-dinitrobenzene, a co-crystal of cis-tartaric acid and bis-pyridine N-oxide and the hydrochloride of B-4-pyridinylboronic acid. We selected these three compounds owing to the π-hole accessibility features that the sp² hybridized B, C and N atoms provide, thus allowing us to compare the fundamental characteristics of π-hole interactions using Bader's Atom in Molecules (AIM) theory. This particular study required extremely accurate experimental diffraction data, because the interaction of interest is weak. As shown by the experimental charge density maps of the -YO₂ (Y = B, C, N) units, we assign the depletion of electron-density present in the central boron, carbon and nitrogen atoms (electrophilic π-holes) as the main origin for the establishment of intermolecular Lewis acid-Lewis base attractive interaction with complementary electron-rich regions. Unexpectedly, the Bader's analyses of both experimentally and theoretically calculated charge distribution maps for the solid involving the - BO₂H₂ group do not show the presence of bond paths, neither of the bond critical points, between the interacting electron rich sites and the boron or carbon atoms featuring the electron hole. In contrast, these topological descriptors of chemical interactions for the AIM theory were easily located in the solid-state structures of the compounds involving the carboxylic and the nitro groups.

Thermodynamics of Anion Binding by Chalcogen Bonding Receptors

The application of chalcogen bonding (ChB) to anion recognition is an underdeveloped area of host-guest supramolecular chemistry. The chemical instability of heavier chalcogen derivatives may in part be responsible for the lack of progress. Herein, the synthesis of a new structurally simple, tellurium-based ChB binding motif is reported, the robust stability of which has enabled the thermodynamic properties for ChB halide anion binding in polar aprotic and wet protic organic solvent media to be elucidated. The thermodynamic data reveals how the subtle interplay between ChB host, anion guest and solvent dictates halide binding selectivity and affinity trends. These findings help to provide a deeper insight into the nature of the ChB-anion interaction.

Tetrel Bonding Interactions in Perchlorinated Cyclopenta- and Cyclohexatetrelanes: A Combined DFT and CSD Study

In this manuscript, we combined DFT calculations (PBE0-D3/def2-TZVP level of theory) and a Cambridge Structural Database (CSD) survey to evaluate the ability of perchlorinated cyclopenta- and cyclohexatetrelanes in establishing tetrel bonding interactions. For this purpose, we used Tr₅Cl₁₀ and Tr₆Cl₁₂ (Tr = Si and Ge) and HCN, HF, OH⁻ and Cl⁻ as electron donor entities. Furthermore, we performed an Atoms in Molecules (AIM) analysis to further describe and characterize the interactions studied herein. A survey of crystal structures in the CSD reveals that close contacts between Si and lone-pair-possessing atoms are quite common and oriented along the extension of the covalent bond formed by the silicon with the halogen atom.

Chalcogen Bonding in Protein-Ligand Complexes: PDB Survey and Quantum Mechanical Calculations

A chalcogen bond is a nonclassical noncovalent interaction which can stabilise small-molecule crystals as well as protein structures. Here, we systematically explore the stabilising potential of chalcogen bonding in protein-ligand complexes in the Protein Data Bank (PDB). We have found that a large fraction (23 %) of complexes with a S/Se-containing ligand feature close S/Se⋅⋅⋅O/N/S contacts. Eleven non-redundant representative potential S/Se⋅⋅⋅O chalcogen-bond motifs were selected and truncated to model systems and seven more model systems were prepared by S-to-Se substitution. These systems were then subjected to analysis by quantum chemical (QM) methods-electrostatic potential, geometry optimisation or interaction energy calculations, including solvent effects. The QM calculations indicate that chalcogen bonding does indeed play a dominant role in stabilising some of the interaction motifs studied. We thus advocate further exploration of chalcogen bonding with the aim of potential future use in structure-based drug design.

Comparison between Tetrel Bonded Complexes Stabilized by σ and π Hole Interactions

The σ-hole tetrel bonds formed by a tetravalent molecule are compared with those involving a π-hole above the tetrel atom in a trivalent bonding situation. The former are modeled by TH₄, TH₃F, and TH₂F₂ (T = Si, Ge, Sn) and the latter by TH₂=CH₂, THF=CH₂, and TF₂=CH₂, all paired with NH₃ as Lewis base. The latter π-bonded complexes are considerably more strongly bound, despite the near equivalence of the σ and π-hole intensities. The larger binding energies of the π-dimers are attributed to greater electrostatic attraction and orbital interaction. Each progressive replacement of H by F increases the strength of the tetrel bond, whether σ or π. The magnitudes of the maxima of the molecular electrostatic potential in the two types of systems are not good indicators of either the interaction energy or even the full Coulombic energy. The geometry of the Lewis acid is significantly distorted by the formation of the dimer, more so in the case of the σ-bonded complexes, and this deformation intensifies the σ and π holes.

Coordination of Be and Mg Centres by HCN Ligands - Be…N and Mg…N Interactions

ωB97XD/aug-cc-pVTZ calculations were performed for clusters of Z2+ cations (Z=Be and Mg) and HCN molecules (up to six molecules). The clusters of Be(CH3 )2 and Mg(CH3 )2 with HCN species were also calculated to analyse the influence of the Be/Mg-C formally covalent bonds on interactions of Be or Mg centre with ligands. The beryllium and magnesium centres possess different areas of a positive electrostatic potential that depend on a number of HCN ligands in the cluster considered. Numerous correlations between geometrical, energetic and topological parameters of the clusters considered are discussed since various theoretical approaches are applied; Quantum Theory of 'Atoms in Molecules', Natural Bond Orbital method and decomposition of the energy of interaction. The Be/Mg…N interactions classified as beryllium and magnesium bonds possess numerous characteristics which are known for the hydrogen bonds. Different types of coordination of Be and Mg centres analysed here exist also in crystal structures.

Regium-π bonds: An Unexplored Link between Noble Metal Nanoparticles and Aromatic Surfaces

The ability of metal clusters involving elements from group 11 (Ag, Cu, Au) to favorably interact with π systems of different size and electronic nature was evaluated at the PBE0-D3/def2-TZVPP//PBE0-D3/def2-TZVP level of theory. The M₉ clusters (M=Cu, Ag, Au) were used as σ-hole and σ-lump donors, and benzene, trifluorobenzene, and hexafluorobenzene as aromatic rings. In addition, the study was expanded to the analysis of extended π systems by using naphthalene and anthracene as well as their corresponding perfluorinated derivatives. Furthermore, Bader's theory of Atoms in Molecules as well as natural bonding orbital and spin-density calculations were used to further investigate and characterize the regium-π and σ-lump complexes described herein. Apparently, regium-π bonds have not previously been described in the literature and may be of great importance in the understanding of organocatalytic processes involving aromatic substrates as well as in the design of new materials based on this novel subclass of σ-hole bonding.

Chalcogen 'like-like' Interactions Involving Trisulphide and Triselenide Compounds: A Combined CSD and Ab Initio Study

In this manuscript, we combined a search in the Cambridge Structural Database (CSD) and ab initio calculations (RI-MP2/def2-TZVPD level of theory) to analyze the ability of trisulphide and triselenide moieties to establish chalcogen ‘like-like’ interactions. A preliminary CSD inspection revealed two predominant structural patterns, depending on the anti or syn conformation adopted by the substituents of the S₃/Se₃ bridge, leading to bifurcated or double chalcogen bonding interactions, respectively. In order to analyze these two relevant structural motifs we have used a series of S and Se derivatives Ch₃X₂ (Ch = S and Se and X = H, F, CN, and CF₃) which act as both electron donor (using the lone pairs) and acceptor (using the σ-holes) entities. Besides, we have carried out “atoms in molecules” (AIM) and natural bonding orbital (NBO) analyses to further describe and characterize the chalcogen bonding interactions described herein. As far as we know, chalcogen···chalcogen interactions involving trichalconides (S₃/Se₃) have not been previously described in literature a may be of great importance in the preparation and characterization of new solids based on this subclass of σ-hole bonding.

Pseudohalide Tectons within the Coordination Sphere of the Uranyl Ion: Experimental and Theoretical Study of C-H···O, C-H···S, and Chalcogenide Noncovalent Interactions

A series of uranyl thiocyanate and selenocyanate of the type [R₄N]₃[UO₂(NCS)₅] (R₄ = ⁿBu₄, Me₃Bz, Et₃Bz), [Ph₄P][UO₂(NCS)₃(NO₃)] and [R₄N]₃[UO₂(NCSe)₅] (R₄ = Me₄, ⁿPr₄, Et₃Bz) have been prepared and structurally characterized. The resulting noncovalent interactions have been examined and compared to other examples in the literature. The nature of these interactions is determined by the cation so that when the alkyl groups are small, chalcogenide···chalcogenide interactions are present, but this "switches off" when R = ⁿPr and charge assisted U═O···H-C and S(e)···H-C hydrogen bonding remain the dominant interaction. Increasing the size of the chain to ⁿBu results in only S···H-C interactions. The spectroscopic implications of these chalcogenide interactions have been explored in the vibrational and photophysical properties of the series [R₄N]₃[UO₂(NCS)₅] (R₄ = Me₄, Et₄, ⁿPr₄, ⁿBu₄, Me₃Bz, Et₃Bz), [R₄N]₃[UO₂(NCSe)₅] (R₄ = Me₄, ⁿPr₄, Et₃Bz) and [Et₄N]₄[UO₂(NCSe)₅][NCSe]. The data suggest that U═O···H-C interactions are weak and do not perturb the uranyl moiety. While the chalcogenide interactions do not influence the photophysical properties, a coupling of the U═O and δ(NCS) or δ(NCSe) vibrational modes is observed in the 77 K solid state emission spectra. A theoretical examination of representative examples of Se···Se, C-H···Se, and C-H···O═U by molecular electrostatic potentials and NBO and AIM methodologies gives a deeper understanding of these weak interactions. C-H···Se are individually weak but C-H···O═U interactions are even weaker, supporting the idea that the -yl oxo's are weak Lewis bases. An Atoms in Molecules study suggests that the chalcogenide interaction is similar to lone pair···π or fluorine···fluorine interactions. An oxidation of the NCS ligands to form [(UO₂)(SO₄)₂(H₂O)₄]·3H₂O was also noted.

Aerogen bonds formed between AeOF2 (Ae = Kr, Xe) and diazines: comparisons between σ-hole and π-hole complexes

The interaction between KrOF₂ or XeOF₂ and the 1,2, 1,3, and 1,4 diazines is characterized chiefly by a Kr/XeN aerogen bond, as deduced from ab initio calculations. The most stable dimers take advantage of the σ-hole on the aerogen atom, wherein the two molecules lie in the same plane. The interaction is quite strong, as much as 18 kcal mol^-1. A second class of dimer geometry utilizes the π-hole above the aerogen atom in an approximate perpendicular arrangement of the two monomers; these structures are not as strongly bound: 6-8 kcal mol^-1. Both sorts of dimers contain auxiliary CHF H-bonds which contribute to their stability, but even with their removal, the aerogen bond energy remains as high as 14 kcal mol^-1. The nature and strength of each specific interaction is confirmed and quantified by AIM, NCI, NBO, and electron density shift patterns. There is not a great deal of sensitivity to the identity of either the aerogen atom or the position of the two N atoms in the diazine.

Concurrent aerogen bonding and lone pair/anion-π interactions in the stability of organoxenon derivatives: a combined CSD and ab initio study

In this manuscript the ability of organoxenon fluorides to establish concurrent π-hole aerogen bonding and lone pair/anion-π interactions has been studied at the RI-MP2/def2-TZVP level of theory. The presence of both an aromatic system (benzene, trifluorobenzene and pentafluorobenzene) and a xenon atom makes these molecules suitable for simultaneously establishing both interactions. In this regard, we have used CH₃CN, NH₃, O(CH₃)₂, Cl^-, CN^- and BF₄^- as neutral and charged electron donors, respectively. Moreover, the NBO analysis showed that orbital effects contribute to the global stabilization of the complexes studied. Furthermore, we have used Bader's theory of "atoms in molecules" to analyse and characterize the noncovalent interactions described herein from a charge-density perspective. Finally, several examples retrieved from the CSD are also included, highlighting the impact of these interactions in the solid state chemistry of Xe.

Sigma-hole carbon-bonding interactions in carbon-carbon double bonds: an unnoticed contact

In this manuscript, we combine high-level ab initio calculations on some small complexes and a CSD survey to analyze the existence of unprecedented noncovalent carbon bonds in X₂C[double bond, length as m-dash]CH₂Y systems (Y = electron-rich atom or group). The methylene group is usually seen as a weak hydrogen bond donor when interacting with an electron-rich atom. However, we demonstrate that when the electron-rich atom is located equidistant from the two H atoms and along the C[double bond, length as m-dash]C bond a σ-hole noncovalent carbon-bonding interaction is established, instead of a bifurcated hydrogen bond, as derived from Atoms-in-Molecules (AIM) and Natural bond orbital (NBO) analyses. The physical nature of the interaction has been analyzed using the Symmetry Adapted Perturbation Theory (SAPT) method. The results indicate that electrostatics is very important followed by either the induction or dispersion terms in anionic and neutral complexes, respectively. In addition the CSD analysis reveals the existence of such interactions, giving reliability to our calculations, which are much more numerous for neutral than for anionic Y systems.

Hydrogen bonds and other interactions as a response to protect doublet/octet electron structure

MP2/aug-cc-pVTZ calculations were performed for complexes linked by hydrogen bonds. Three types of proton donating species were taken into account: H₂O, CCl₃H, and H₃O⁺. These calculations are supported by the natural bond orbital (NBO) method and the quantum theory of atoms in molecules (QTAIM) approach. Numerous correlations between parameters of H-bonded systems were found. The most important are those which show the response of the system on the H-bond formation; for example, the increase of polarization of the A-H bond correlates with the strength of the hydrogen bond. Similar relationships were found for the σ-hole bonds while the π-hole bonds do not follow the trends known for the hydrogen bonds. Graphical abstract Hydrogen bonds and other interactions as a response to protect doublet/octet electron structureᅟ.

Implications of monomer deformation for tetrel and pnicogen bonds

Monomer rearrangement raises the interaction energy by up to 20 kcal mol⁻¹and intensifies its σ-hole by a factor of 1.5–2.9.

Intermolecular and very strong intramolecular C–Se⋯O/N chalcogen bonds in nitrophenyl selenocyanate crystals

Different bonding strengths of C–Se⋯O/N chalcogen bonds involved in polymorphic o-NSC (1a/1b) and monomorphic p-NSC (2) result in different thermal properties.

Structure guided or structure guiding? Mixed carbon/hydrogen bonding in a bis-Schiff base of N-allyl isatin

A supramolecular motif listed as ‘carbon bonded’ or ‘hydrogen bonded’ may have the character of both. We highlight the hybrid character of the non-covalent interaction in a bis-Schiff base of N-allyl isatin by combining theory and experiment.

Non-covalent tetrel bonding interactions in hemidirectional lead(ii) complexes with nickel(ii)-salen type metalloligands

Tetrel bonding interactions have been investigated in hetero-dinuclear nickel(ii)/lead(ii) complexes using MEP and DFT calculations.

On the σ, π and δ hole interactions: a molecular orbital overview

With molecular orbital theory it is possible to distinguish and design σ, π and the elusive δ electrostatic holes.

On the importance of Pb⋯X (X = O, N, S, Br) tetrel bonding interactions in a series of tetra- and hexa-coordinated Pb(ii) compounds

We study Pb⋯X (X = N, S, Br) tetrel bonding interactions in a series of tetra- and hexa-coordinated Pb(ii) compounds that exhibit hemidirectional coordination.

Pb⋯X (X = N, S, I) tetrel bonding interactions in Pb(ii) complexes: X-ray characterization, Hirshfeld surfaces and DFT calculations

We report the synthesis and X-ray characterization of four new Pb(ii) complexes of nicotinoylhydrazone and picolinoylhydrazone-based ligands and three different anionic co-ligands (acetate, thiocyanate and iodide) exhibiting relevant tetrel bonding interactions.

New Type of Halogen Bond: Multivalent Halogen Interacting with π- and σ-Electrons

MP2/aug-cc-pVTZ calculations were performed for complexes of BrF₃ and BrF₅ acting as Lewis acids through the bromine centre, with species playing a role of Lewis base: dihydrogen, acetylene, ethylene, and benzene. The molecular hydrogen donates electrons by its σ-bond, while in remaining moieties—in complexes of hydrocarbons; such an electron transfer follows from π-electrons. The complexes are linked by a kind of the halogen bond that is analyzed for the first time in this study, i.e., it is the link between the multivalent halogen and π or σ-electrons. The nature of such a halogen bond is discussed, as well as various dependencies and correlations are presented. Different approaches are applied here, the Quantum Theory of Atoms in Molecules, Natural Bond Orbital method, the decomposition of the energy of interaction, the analysis of electrostatic potentials, etc.

Systematic Coupled Cluster Study of Noncovalent Interactions Involving Halogens, Chalcogens, and Pnicogens

The noncovalent interactions of 32 complexes involving pnicogens, chalcogens, and halogens atoms were investigated at the CCSD(T)/aug-cc-pVTZ level of theory. Two different types of complexes could be distinguished on the basis of geometric parameters, electron difference densities, and the charge transfer mechanisms associated with each type. In the type I conformation, the monomers adopt a skewed orientation allowing charge to be transfer between both monomers, whereas in the type II conformation the monomers adopt a linear arrangement, maximizing charge transfer in only one direction. Type I complexes involving the interaction between pnicogens and chalcogens cannot be unambiguously defined as chalcogen or pnicogen bonds, they are an admixture of both. The charge transfer dependence on the conformation adopted by the complexes described in this work can serve as a novel conformationally driven design concept for materials.

The Origin of Chalcogen-Bonding Interactions

Favorable molecular interactions between group 16 elements have been implicated in catalysis, biological processes, and materials and medicinal chemistry. Such interactions have since become known as chalcogen bonds by analogy to hydrogen and halogen bonds. Although the prevalence and applications of chalcogen-bonding interactions continues to develop, debate still surrounds the energetic significance and physicochemical origins of this class of σ-hole interaction. Here, synthetic molecular balances were used to perform a quantitative experimental investigation of chalcogen-bonding interactions. Over 160 experimental conformational free energies were measured in 13 different solvents to examine the energetics of O···S, O···Se, S···S, O···HC, and S···HC contacts and the associated substituent and solvent effects. The strongest chalcogen-bonding interactions were found to be at least as strong as conventional H-bonds, but unlike H-bonds, surprisingly independent of the solvent. The independence of the conformational free energies on solvent polarity, polarizability, and H-bonding characteristics showed that electrostatic, solvophobic, and van der Waals dispersion forces did not account for the observed experimental trends. Instead, a quantitative relationship between the experimental conformational free energies and computed molecular orbital energies was consistent with the chalcogen-bonding interactions being dominated by n → σ* orbital delocalization between a lone pair (n) of a (thio)amide donor and the antibonding σ* orbital of an acceptor thiophene or selenophene. Interestingly, stabilization was manifested through the same acceptor molecular orbital irrespective of whether a direct chalcogen···chalcogen or chalcogen···H-C contact was made. Our results underline the importance of often-overlooked orbital delocalization effects in conformational control and molecular recognition phenomena.

Two faces of triel bonds in boron trihalide complexes

The N⋅⋅⋅B triel bonds in complexes of boron trihalides, BX₃ (X = F, Cl, Br, and I), with species acting as Lewis bases through the nitrogen center, NH₃ , N₂ , and HCN, are analyzed theoretically (MP2/aug-cc-pVTZ calculations). It is confirmed that stronger Lewis acid properties of the boron center are observed for the BCl₃ moiety than for the BF₃ one in complexes with the strong Lewis base (NH₃ ); while the opposite order is observed for complexes with the weak Lewis base (N₂ ). The BX₃ NCH complexes (for X = Cl, Br, and I) are characterized by two tautomeric forms and by two corresponding N⋅⋅⋅B distances, the shorter one possesses characteristics of the covalent bond. In a case of the BF₃ NCH complex one energetic minimum is observed. Ab initio calculations are supported by an analysis of molecular electrostatic potentials (EPs) and electron density distributions. The quantum theory of 'atoms in molecules' and the decomposition of the energy of interaction are applied. The aforementioned acidity orders as well as the existence of two tautomers for some of complexes result partly from the electrostatic interactions' balance; the EP distribution is different for the BF₃ species than for the other BX₃ species where X = Cl, Br, and I. © 2017 Wiley Periodicals, Inc.

Cooperativity between the hydrogen bonding and σ-hole interaction in linear NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 complexes (X = Cl, Br; Z = Ar, Kr): a comparative study

Quantum chemical calculations are performed to investigate the cooperativity of hydrogen bonding with halogen or aerogen bonding interactions in linear NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 clusters, where X = Cl, Br and Z = Ar, Kr. To understand the cooperativity mechanism in these systems, the corresponding binary NCX···NCH and O3Z···NCH complexes are also considered. The binding distances, interaction energies, and bonding properties of the NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 clusters are analyzed in detail. It is found that the cooperative effects in the hydrogen bonding tend to strengthen X···N and Z···N interactions. For both NCX···(NCH)n and O3Z···(NCH)n clusters, a small bond shrinkage is observed from n = 4 to n = 5, which suggests that the cooperativity effects are almost saturated in the larger clusters (n > 5). As the size of the X or Z atom is increased, the magnitude of the cooperative energy in these systems is also increased, which is mainly ascribed to changes in electrostatic potentials and orbital interactions. Our results indicate that the cooperative effects lead to a substantial change in the 14N nuclear quadrupole coupling constants of the NCH molecule.

Halogen and Hydrogen Bonding between (N-Halogeno)-succinimides and Pyridine Derivatives in Solution, the Solid State and In Silico

A study of strong halogen bonding within three series of halogen-bonded complexes, derived from seven para-substituted pyridine derivatives and three N-halosuccinimides (iodo, bromo and chloro), has been undertaken with the aid of single-crystal diffraction, solution complexation and computational methods. The halogen bond was compared with the hydrogen bond in an equivalent series based on succinimide. The halogen-bond energies are in the range -60 to -20 kJ mol^-1 and change regularly with pyridine basicity and the Lewis acidity of the halogen. The halogen-bond energies correlate linearly with the product of charges on the contact atoms, which indicates a predominantly electrostatic interaction. The binding enthalpies in solution are around 19 kJ mol^-1 less negative due to solvation effects. The optimised geometries of the complexes in the gas phase are comparable to those of the solid-state structures, and the effects of the supramolecular surroundings in the latter are discussed. The bond energies for the hydrogen-bonded series are intermediate between the halogen-bond energies of iodine and bromine, although there are specific differences in the geometries of the halogen- and hydrogen-bonded complexes.