Noncovalent interactions in dimers and trimers of SO3 and CO

Chalcogen Noncovalent Interactions between Diazines and Sulfur Oxides in Supramolecular Circular Chains

The noncovalent chalcogen interaction between SO2/SO3 and diazines was studied through a dispersion-corrected DFT Kohn-Sham molecular orbital together with quantitative energy decomposition analyses. For this, supramolecular circular chains of up to 12 molecules were built with the aim of checking the capability of diazine molecules to detect SO2/SO3 compounds within the atmosphere. Trends in the interaction energies with the increasing number of molecules are mainly determined by the Pauli steric repulsion involved in these σ-hole/π-hole interactions. But more importantly, despite the assumed electrostatic nature of the involved interactions, the covalent component also plays a determinant role in its strength in the involved chalcogen bonds. Noticeably, π-hole interactions are supported by the charge transfer from diazines to SO2/SO3 molecules. Interaction energies in these supramolecular complexes are not only determined by the S···N bond lengths but attractive electrostatic and orbital interactions also determine the trends. These results should allow us to establish the fundamental characteristics of chalcogen bonding based on its strength and nature, which is of relevance for the capture of sulfur oxides.

Ab initio study for superior sensitivity of graphyne nanoflake towards nitrogen halides over ammonia

Graphyne (GYN) has received immense attention in gas adsorption applications due to its large surface area. The adsorption of toxic ammonia and nitrogen halides gaseous molecules on graphyne has been theoretically studied at ωB97XD/6-31 + G(d, p) level of DFT. The counterpoise corrected interaction energies of NH₃, NF₃, NCl₃, and NBr₃ molecules with GYN are - 4.73, - 2.27, - 5.22, and - 7.19 kcal mol^-1, respectively. Symmetry-adapted perturbation theory (SAPT0) and noncovalent interaction index (NCI) reveal that the noncovalent interaction between analytes and GYN is dominated by dispersion forces. The significant change in electronic behavior, i.e., energies of HOMO and LUMO orbitals and NBO charge transfer correspond to the pronounced sensitivity of GYN towards considered analytes, especially NBr₃. Finally, TD-DFT calculation reveals a decrease in electronic transition energies and shifting of adsorption to a longer wavelength. The recovery time for NX₃@GYN is observed in nanoseconds, which is many orders of magnitude smaller than the reported systems. The recovery time is further decreased with increasing temperature, indicating that the GYN benefits from a short recovery time as a chemical sensor.

CHAL336 Benchmark Set: How Well Do Quantum-Chemical Methods Describe Chalcogen-Bonding Interactions?

We present the CHAL336 benchmark set-the most comprehensive database for the assessment of chalcogen-bonding (CB) interactions. After careful selection of suitable systems and identification of three high-level reference methods, the set comprises 336 dimers each consisting of up to 49 atoms and covers both σ- and π-hole interactions across four categories: chalcogen-chalcogen, chalcogen-π, chalcogen-halogen, and chalcogen-nitrogen interactions. In a subsequent study of DFT methods, we re-emphasize the need for using proper London dispersion corrections when treating noncovalent interactions. We also point out that the deterioration of results and systematic overestimation of interaction energies for some dispersion-corrected DFT methods does not hint at problems with the chosen dispersion correction but is a consequence of large density-driven errors. We conclude this work by performing the most detailed DFT benchmark study for CB interactions to date. We assess 109 variations of dispersion-corrected and dispersion-uncorrected DFT methods and carry out a detailed analysis of 80 of them. Double-hybrid functionals are the most reliable approaches for CB interactions, and they should be used whenever computationally feasible. The best three double hybrids are SOS0-PBE0-2-D3(BJ), revDSD-PBEP86-D3(BJ), and B2NCPLYP-D3(BJ). The best hybrids in this study are ωB97M-V, PW6B95-D3(0), and PW6B95-D3(BJ). We do not recommend using the popular B3LYP functional nor the MP2 approach, which have both been frequently used to describe CB interactions in the past. We hope to inspire a change in computational protocols surrounding CB interactions that leads away from the commonly used, popular methods to the more robust and accurate ones recommended herein. We would also like to encourage method developers to use our set for the investigation and reduction of density-driven errors in new density functional approximations.

Effective adsorption of A-series chemical warfare agents on graphdiyne nanoflake: a DFT study

Chemical warfare agents (CWAs) are highly poisonous and their presence may cause diverse effects not only on living organisms but also on environment. Therefore, their detection and removal in a short time span is very important. In this regard, here the utility of graphdiyne (GDY) nanoflake is studied theoretically as an electrochemical sensor material for the hazardous CWAs including A-230, A-232, and A-234. Herein, we explain the phenomenon of adsorption of A-series CWAs on GDY nanoflake within the density functional theory (DFT) framework. The characterisation of adsorption is based on optimised geometries, BSSE-corrected energies, SAPT0, RDG, FMO, CHELPG charge transfer, QTAIM and UV-Vis analyses. The calculated counterpoise adsorption energies for reported complexes range from - 13.70 to - 17.19 kcal mol^-1. These adsorption energies show that analytes are physiosorbed onto GDY which usually takes place through noncovalent interactions. The noncovalent adsorption of CWAs on GDY is also attributed by the SAPT0, RDG and QTAIM analyses. These properties also reveal that dispersion factors dominate in the complexes among many noncovalent components (exchange, induction, electrostatic, steric and repulsion). In order to estimate the sensitivity of GDY, the %sensitivity and average energy gap variations are quantitatively measured by energies of HOMO and LUMO orbitals. In terms of adsorption affinity of GDY, UV-Vis analysis, CHELPG charge transfer and DOS analyses depict an appreciable response towards these toxic CWAs. Graphical abstract.

Participation of S and Se in hydrogen and chalcogen bonds

The heavier chalcogen atoms S, Se, and Te can each participate in a range of different noncovalent interactions. They can serve as both proton donor and acceptor in H-bonds. Each atom can also act as electron acceptor in a chalcogen bond.

Local Vibrational Mode Analysis of π–Hole Interactions between Aryl Donors and Small Molecule Acceptors

11 aryl–lone pair and three aryl–anion π –hole interactions are investigated, along with the argon–benzene dimer and water dimer as reference compounds, utilizing the local vibrational mode theory, originally introduced by Konkoli and Cremer, to quantify the strength of the π –hole interaction in terms of a new local vibrational mode stretching force constant between the two engaged monomers, which can be conveniently used to compare different π –hole systems. Several factors have emerged which influence strength of the π –hole interactions, including aryl substituent effects, the chemical nature of atoms composing the aryl rings/ π –hole acceptors, and secondary bonding interactions between donors/acceptors. Substituent effects indirectly affect the π –hole interaction strength, where electronegative aryl-substituents moderately increase π –hole interaction strength. N-aryl members significantly increase π –hole interaction strength, and anion acceptors bind more strongly with the π –hole compared to charge neutral acceptors (lone–pair donors). Secondary bonding interactions between the acceptor and the atoms in the aryl ring can increase π –hole interaction strength, while hydrogen bonding between the π –hole acceptor/donor can significantly increase or decrease strength of the π –hole interaction depending on the directionality of hydrogen bond donation. Work is in progress expanding this research on aryl π –hole interactions to a large number of systems, including halides, CO, and OCH3− as acceptors, in order to derive a general design protocol for new members of this interesting class of compounds.

Regium Bonds between Silver(I) Pyrazolates Dinuclear Complexes and Lewis Bases (N2, OH2, NCH, SH2, NH3, PH3, CO and CNH)

A theoretical study and Cambridge Structural Database (CSD) search of dinuclear Ag(I) pyrazolates interactions with Lewis bases were carried out and the effect of the substituents and ligands on the structure and on the aromaticity were analyzed. A relationship between the intramolecular Ag–Ag distance and stability was found in the unsubstituted system, which indicates a destabilization at longer distances compensated by ligands upon complexation. It was also observed that the asymmetrical interaction with phosphines as ligands increases the Ag–Ag distance. This increase is dramatically higher when two simultaneous PH3 ligands are taken into account. The calculated 109Ag chemical shielding shows variation up to 1200 ppm due to the complexation. Calculations showed that six-membered rings possessed non-aromatic character while pyrazole rings do not change their aromatic character significantly upon complexation.

(Isocyano group)⋯lone pair interactions involving coordinated isocyanides: experimental, theoretical and CSD studies

Both carbon and nitrogen centers in the coordinated isocyano group are capable of acting as a π-hole donor toward lone pairs.

Theoretical Investigation of Cyano-Chalcogen Dimers and Their Importance in Molecular Recognition

In this manuscript the different noncovalent interactions established between (HYCN)₂ dimers (Y=S, Se and Te) have been studied at the MP2 and CCSD(T) level of theory. Several homodimers have been taken into account, highlighting the capacity of these compounds to act both as electron donor and acceptor. The main properties studied were geometries, binding energy (E_b ), and molecular electrostatic potential (MEP). Given the wide application of chalcogen bonds, and more specifically of cyano-chalcogen moieties in molecular recognition, natural bond orbital (NBO), "atoms-in-molecules" (AIM), and electron density shift (EDS) analysis were also used to analyse the different noncovalent interactions upon complexation. The presence of hydrogen, chalcogen and dipole-dipole interactions was confirmed and their implications on molecular recognition were analysed.

Cooperative Effects in Weak Interactions: Enhancement of Tetrel Bonds by Intramolecular Hydrogen Bonds

A series of silyl and germanium complexes containing halogen atoms (fluorine and chlorine atoms) and exhibiting tetrel bonds with Lewis bases were analyzed by means of Møller-Plesset computational theory. Binding energies of germanium derivatives were more negative than silicon ones. Amongst the different Lewis bases utilized, ammonia produced the strongest tetrel bonded complexes in both Ge and Si cases, and substitution of the F atom by Cl led to stronger complexes with an ethylene backbone. However, with phenyl backbones, the fluorosilyl complexes were shown to be less stable than the chlorosilyl ones, but the opposite occurred for halogermanium complexes. In all the cases studied, the presence of a hydroxyl group enhanced the tetrel bond. That effect becomes more remarkable when an intramolecular hydrogen bond between the halogen and the hydrogen atom of the hydroxyl group takes places.

Chalcogen Bonding "2S-2N Squares" versus Competing Interactions: Exploring the Recognition Properties of Sulfur

Chalcogen bonding (CB) is the focus of increased attention for its applications in medicinal chemistry, materials science, and crystal engineering. However, the origin of sulfur's recognition properties remains controversial, and experimental evidence for supporting theories is still emerging. Here, a comprehensive evaluation of sulfur CB interactions is presented by investigating 2,1,3-benzothiadiazole X-ray crystallographic structures gathered from the Cambridge Structure Database (CSD), Protein Data Bank (PDB), and own laboratory findings. Through the systematic analysis of substituent effects on a subset library of over thirty benzothiadiazole derivatives, the competing interactions have been categorized into four main classes, namely 2S-2N CB square, halogen bonding (XB), S⋅⋅⋅S, and hydrogen-bonding (HB). A geometric model is employed to characterize the 2S-2N CB square motifs and discuss the role of electrostatic, dipole, and orbital contributions toward the interaction.

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.

Ligation-Enhanced π-Hole···π Interactions Involving Isocyanides: Effect of π-Hole···π Noncovalent Bonding on Conformational Stabilization of Acyclic Diaminocarbene Ligands

The reaction of cis-[PdCl₂(CNXyl)₂] (Xyl = 2,6-Me₂C₆H₃) with the aminoazoles [1 H-imidazol-2-amine (1), 4 H-1,2,4-triazol-3-amine (2), 1 H-tetrazol-5-amine (3), 1 H-benzimidazol-2-amine (4), 1-alkyl-1 H-benzimidazol-2-amines, where alkyl = Me (5), Et (6)] in a 2:1 ratio in the presence of a base in CHCl₃ at RT proceeds regioselectively and leads to the binuclear diaminocarbene complexes [(ClPdCNXyl)₂{μ-C(N-azolyl)N(Xyl)C═NXyl}] (7-12; 73-91%). Compounds 7-12 were characterized by C, H, N elemental analyses, high-resolution ESI⁺-MS, Fourier transform infrared spectroscopy, 1D (¹H, ¹³C) and 2D (¹H,¹H-COSY, ¹H,¹H-NOESY, ¹H,¹³C-HSQC, ¹H,¹³C-HMBC) NMR spectroscopies, and X-ray diffraction (XRDn). Inspection of the XRDn data and results of the Hirshfeld surface analysis suggest the presence in all six structures of intramolecular π-hole_isocyanide···π_arene interactions between the electrophilic C atom of the isocyanide moiety and the neighboring arene ring. These interactions also result in distortion of the Pd-C≡N-Xyl fragment from the linearity. Results of density functional theory calculations [M06/MWB28 (Pd) and 6-31G* (other atoms) level of theory] for model structures of 7-9 followed by the topological analysis of the electron density distribution within the framework of Bader's theory (QTAIM method) reveal the presence of these weak interactions also in a CHCl₃ solution, and their calculated strength is 1.9-2.2 kcal/mol. The natural bond orbital analysis of 7-9 revealed that π(C-C)Xyl → π*(C-N)isocyanide charge transfer (CT) takes place along with the intramolecular π-hole_isocyanide···π_arene interactions. The observed π(C-C)_Xyl → π*(C-N)_isocyanide CT is due to ligation of the isocyanide to the metal center, whereas in the cases of the uncomplexed p-CNC₆H₄NC and CNXyl species, the effects of CT are negligible. Available CCDC data were processed from the perspective of isocyanide-involving π-hole···π interactions, disclosed the role of metal coordination in the π-hole donor ability of isocyanides, and verified the π-hole_isocyanide···π_arene interaction effect on the stabilization of the in-conformation in metal-bound acyclic diaminocarbenes.

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.

Nature of bonding and cooperativity in linear DMSO clusters: A DFT, AIM and NCI analysis

This study aims to cast light on the nature of interactions and cooperativity that exists in linear dimethyl sulfoxide (DMSO) clusters using dispersion corrected density functional theory. In the linear DMSO, DMSO molecules in the middle of the clusters are bound strongly than at the terminal. The plot of the total binding energy of the clusters vs the cluster size and mean polarizabilities vs cluster size shows an excellent linearity demonstrating the presence of cooperativity effect. The computed incremental binding energy of the clusters remains nearly constant, implying that DMSO addition at the terminal site can happen to form an infinite chain. In the linear clusters, two σ-hole at the terminal DMSO molecules were found and the value on it was found to increase with the increase in cluster size. The quantum theory of atoms in molecules topography shows the existence of hydrogen and SO⋯S type in linear tetramer and larger clusters. In the dimer and trimer SO⋯OS type of interaction exists. In 2D non-covalent interactions plot, additional peaks in the regions which contribute to the stabilization of the clusters were observed and it splits in the trimer and intensifies in the larger clusters. In the trimer and larger clusters in addition to the blue patches due to hydrogen bonds, additional, light blue patches were seen between the hydrogen atom of the methyl groups and the sulphur atom of the nearby DMSO molecule. Thus, in addition to the strong H-bonds, strong electrostatic interactions between the sulphur atom and methyl hydrogens exists in the linear clusters.

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.

Theoretical insights into the π-hole interactions in the complexes containing triphosphorus hydride (P3H3) and its derivatives

The π-hole of triphosphorus hydride (P₃H₃) and its derivatives Z₃X₃ (Z = P, As; X = H, F, Cl, Br) was discovered and analyzed. MP2/aug-cc-pVDZ calculations were performed on the π-hole interactions in the HCN...Z₃X₃ complexes and the mutual influence between π-hole interactions and the hydrogen bond in the HCN...HCN...Z₃X₃ and HCN...Z₃X₃...HCN complexes studied. The π-hole interaction belongs to the typical closed-shell noncovalent interaction. The linear relationship was found between the most positive electrostatic potential of the π-hole (V_S,max) and the interaction energy. Moreover, the V_S,max of the π-hole was also found to be linearly correlated to the electrostatic energy term, indicating the important contribution of the electrostatic energy term to the π-hole interaction. There is positive cooperativity between the π-hole interaction and the hydrogen bond in the termolecular complexes. The π-hole interaction has a greater influence on the hydrogen bond than vice versa. The mutual enhancing effect between the π-hole interaction and the hydrogen bond in the HCN...HCN...Z₃X₃ complexes is greater than that in the HCN...Z₃X₃...HCN complexes.

Difference in Energy between Two Distinct Types of Chalcogen Bonds Drives Regioisomerization of Binuclear (Diaminocarbene)PdII Complexes

The reaction of cis-[PdCl₂(CNXyl)₂] (Xyl = 2,6-Me₂C₆H₃) with various 1,3-thiazol- and 1,3,4-thiadiazol-2-amines in chloroform gives a mixture of two regioisomeric binuclear diaminocarbene complexes. For 1,3-thiazol-2-amines the isomeric ratio depends on the reaction conditions and kinetically (KRs) or thermodynamically (TRs) controlled regioisomers were obtained at room temperature and on heating, respectively. In CHCl₃ solutions, the isomers are subject to reversible isomerization accompanied by the cleavage of Pd-N and C-N bonds in the carbene fragment XylNCN(R)Xyl. Results of DFT calculations followed by the topological analysis of the electron density distribution within the formalism of Bader's theory (AIM method) reveal that in CHCl₃ solution the relative stability of the regioisomers (ΔG_exp = 1.2 kcal/mol; ΔG_calcd = 3.2 kcal/mol) is determined by the energy difference between two types of the intramolecular chalcogen bonds, viz. S···Cl in KRs (2.8-3.0 kcal/mol) and S···N in TRs (4.6-5.3 kcal/mol). In the case of the 1,3,4-thiadiazol-2-amines, the regioisomers are formed in approximately equal amounts and, accordingly, the energy difference between these species is only 0.1 kcal/mol in terms of ΔG_exp (ΔG_calcd = 2.1 kcal/mol). The regioisomers were characterized by elemental analyses (C, H, N), HRESI⁺-MS and FTIR, 1D (¹H, ¹³C{¹H}) and 2D (¹H,¹H-COSY, ¹H,¹H-NOESY, ¹H,¹³C-HSQC, ¹H,¹³C-HMBC) NMR spectroscopies, and structures of six complexes (three KRs and three TRs) were elucidated by single-crystal X-ray diffraction.

On the Versatility of BH2 X (X=F, Cl, Br, and I) Compounds as Halogen-, Hydrogen-, and Triel-Bond Donors: An Ab Initio Study

In this manuscript, the ability of BH₂ X compounds (X=F, Cl, Br, and I) to establish halogen-, hydrogen-, and triel-bonding interactions was studied at the RI-MP2/aug-cc-pVQZ level of theory. Several homodimers were taken into account, highlighting the versatility of these compounds to act as both electron donors and electron acceptors. Natural bond orbital analysis showed that orbital effects were important contributors to the global stabilization of the σ- and π-hole bonded complexes studied. Finally, some X-ray solid-state structures retrieved from the Cambridge structural database were described to demonstrate the importance of these interactions involving boron derivatives in the solid state.

π-Hole aerogen bonding interactions

In this manuscript we combine high level ab initio calculations (RI-MP2/aug-cc-pVTZ) and the analysis of several crystal structures to demonstrate the existence of π-hole aerogen bonding interactions in Xe(iv) compounds. The ability of XeF4 and Xe(OMe)4 to interact with electron rich molecules is rationalized using several computational tools, including molecular electrostatic potential surfaces, energetic and geometric features of the complexes and "atoms in molecules" (AIM) and Natural Bond Orbital (NBO) analyses. We have found support for the π-hole interaction involving the xenon atom from the solid state architecture of several X-ray structures retrieved from the crystal structural depot. Particularly, π-hole aerogen bonding interactions are quite common in the solid state of Xe(IV) compounds.

Expanding the applicability of electrostatic potentials to the realm of transition states

Central to any reaction mechanism study, and sometimes a challenging job, is tracing a transition state in a reaction path. For the first time, electrostatic potentials (ESP) of the reactants were used as guiding tactics to predict whether there is a possibility of any transition state in a reaction surface. The main motive behind this strategy is to see whether the directionality nature of the transition state has something to do with the anisotropic natures of the ESP with their embedded directionalities. Strategically, some atmospherically important, but simple, reactions have been chosen for this study, which heretofore were believed to be barrierless. By carefully analysing the ESP maps of the reactants, regions of possible interactions were located. Using the bilinear interpolation of the 2D grids of the ESP surfaces, search co-ordinates were fine-tuned for a local gradient based approach for the search of a transition state. Out of the three reactions studied in this work, we were able to successfully locate transition states, for the first time, in two cases and the third one still proved to be barrierless. This gives a clear indication that though ESP maps can qualitatively predict the possibility of a transition state; it is not always true that there should definitely be a transition state, as some of the reaction surfaces may genuinely be barrierless. But, nevertheless this strategy definitely has credential to be tested for many more reactions, either new or already established, and may be applied to create the initial search co-ordinates for any well-established transition state search method. Moreover, we have observed that the analysis of the ESP maps of the reactants were very much useful in explaining the nature of interactions existing in those observed transition states and we hope the same can also be extended to any transition state in an electrostatically driven reaction potential energy surface.

Theoretical study on σ- and π-hole carbon⋯carbon bonding interactions: implications in CFC chemistry

In this manuscript the ability of CO₂ and several CFCs to establish noncovalent carbon⋯carbon interactions with atmospheric gases CO, ethene and ethyne has been studied at the RI-MP2/def2-TZVPD level of theory.

Weak interactions within nitryl halide heterodimers

Nitryl halides (XNO₂, X = F, Cl, Br and I) are versatile molecules that exhibit several types of interactions within XNO₂:YNO₂ heterodimers mainly governed by dispersion.

Exploration of unconventional π–hole and C–H⋯H–C types of supramolecular interactions in a trinuclear Cd(ii) and a heteronuclear Cd(ii)–Ni(ii) complex and experimental evidence for preferential site selection of the ligand by 3d and 4d metal ions

We report the synthesis of a trinuclear cadmium(ii) (1) and a di(phenoxido)-bridged dinuclear cadmium(ii)–nickel(ii) (2) complex derived from a bicompartmental (N₂O₄) Schiff base ligand.

Comparison of π-hole tetrel bonding with σ-hole halogen bonds in complexes of XCN (X = F, Cl, Br, I) and NH3

In addition to the standard halogen bond formed when NH₃approaches XCN (X = F, Cl, Br, I) along its molecular axis, a perpendicular approach is also possible, toward a π-hole that is present above the X–C bond.

A Study of π-π Stacking Interactions and Aromaticity in Polycyclic Aromatic Hydrocarbon/Nucleobase Complexes

We analysed the interactions and aromaticity electron-density delocalisation observed in π-π complexes between the phenalenyl radical and acenaphthylene, and the DNA and RNA nucleobases (adenine, guanine, cytosine, thymine and uracil). Interaction energies are obtained at the M06-2X/6-311++G(2df,p) computational level for gas phase and PCM-water conditions. For both the phenalenyl radical and acenaphthylene, the complexes formed with guanine are the most stable ones. Atoms in molecules and natural bond orbital results reveal weak π-π interactions between both interacting moieties, characterized by bond critical points between C⋅⋅⋅C and C⋅⋅⋅N atoms. Nucleus independent chemical shifts (NICS) indicate the retention of the aromatic character of the monomers in the outer region of the complex. The fluctuation indexes reveal a loss of electron delocalisation upon complexation for all cases except guanine complexes. Nevertheless, the interface region shows large negative NICS values, which is not associated with an increase of the aromaticity or electron-density delocalisation, but with magnetic couplings of both molecules, leading to an unrealistic description of the aromatic behaviour in that region.

1,1,2,2-Tetracyanocyclopropane (TCCP) as supramolecular synthon

The 1,1,2,2-tetracyanocyclopropane (TCCP) unit presents a synthetically accessible and versatile synthon that can interact with lone-pair or π-electrons by 'non-covalent carbon bonding'. Complexes of TCCP with common small molecules, anions, aromatics like fullerenes, amino acids and nucleobases were computed at the DFT BP86-D3/def2-TZVP level of theory. Binding energies vary between about -10 kcal mol(-1) for neutral guests and -15 to -50 kcal mol(-1) for anionic species. This is comparable to strong and very strong hydrogen bonding respectively. Thus, in addition to synthons that contain polarized hydrogen or halogen atoms, TCCP presents a new supramolecular synthon that awaits experimental exploitation.