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.

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.

Understanding noncovalent bonds and their controlling forces

The fundamental underpinnings of noncovalent bonds are presented, focusing on the σ-hole interactions that are closely related to the H-bond. Different means of assessing their strength and the factors that control it are discussed. The establishment of a noncovalent bond is monitored as the two subunits are brought together, allowing the electrostatic, charge redistribution, and other effects to slowly take hold. Methods are discussed that permit prediction as to which site an approaching nucleophile will be drawn, and the maximum number of bonds around a central atom in its normal or hypervalent states is assessed. The manner in which a pair of anions can be held together despite an overall Coulombic repulsion is explained. The possibility that first-row atoms can participate in such bonds is discussed, along with the introduction of a tetrel analog of the dihydrogen bond.

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.

Weak Interactions Get Strong: Synergy between Tetrel and Alkaline-Earth Bonds

Weak and strong noncovalent interactions such as tetrel bonds and alkaline-earth bonds, respectively, cooperate and get reinforced when acting together in ternary complexes of general formula RN··· SiH₃F···MY, where MY is a Be or Mg derivative and RN is a N-containing Lewis base with different hybridization patterns. Cooperativity has been studied in the optimized MP2/aug'-cc-pVTZ ternary complexes by looking at changes on geometries, binding energies, ²⁹Si NMR chemical shifts, and topological features according to the atoms in molecules theoretical framework. Our study shows that cooperativity in terms of energy is in general significant: more than 40 kJ/mol, and up to 83.6 kJ/mol in the most favorable case. The weakest the isolated interaction, the strongest the reinforcement in the ternary complex; in this sense, the tetrel bond is shortened enormously, between 0.3 and 0.6 Å. This dramatic reinforcement of the tetrel bond is also nicely reflected in the positive variations of the ²⁹Si chemical shifts in all the ternary complexes. At the same time the ternary complexes are characterized by the presence of totally planar silyl group, due to the pentacoordination of the Si atom. Both the hybridization of the N base and the geometry imposed by the alkaline-earth ligands have a strong influence on the binding energies, as they modify the donor ability of N and the Lewis acid character of the alkaline-earth metal.

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.

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.

Alkaline-earth (Be, Mg and Ca) bonds at the origin of huge acidity enhancements

The interaction between alkaline-earth derivatives with the general formula X₂M (X = H, F and Cl; M = Be, Mg and Ca) and a set of Lewis bases, including first and second-row hydrides, namely YH_n (Y = O, N, F, S, P and Cl) hydrides, as well as other typical cyclic organic bases, such as aniline, 1H-1,2,3-triazole, 1H-tetrazole and phenylphosphine, was investigated using the G4 ab initio composite method. Contrary to what was expected, it was found that the interactions involving Mg and Ca derivatives were not necessarily weaker than those between beryllium bonds. The origin is two-fold: larger deformation of the interacting systems when Be-derivatives are involved and appearance of secondary non-covalent interactions in the formation of some of the Mg- and Ca-containing complexes. Hence, the dissociation of the latter complexes may require higher enthalpies than that of the Be complexes. These deformations are triggered by a significant redistribution of electron density of the two interacting moieties, which also result in dramatic changes in the reactivity of the interacting compounds and in particular in the intrinsic basicity of the Lewis bases investigated, to the point that conventional bases, such as ammonia or aniline, upon complexation with MCl₂ (M = Be, Mg and Ca), become stronger Brønsted acids than phosphoric acid, whereas other bases, such as 1H-tetrazole, become stronger acids than perchloric acid.

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.

Nature and strength of chalcogen–π bonds

We have analyzed the chalcogen–π bonding mechanism in a systematic series of model systems through Kohn–Sham molecular orbital theory and a quantitative energy decomposition scheme.

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.

The mutual influence of Y···N and H···H interactions in XHY···NCH···HM complexes (X = F, Cl, Br; Y = S, Se; M = Li, Na, BeH, MgH): tuning of the chalcogen bond by dihydrogen bond interaction

The equilibrium structures, interaction energies, and bonding properties of ternary XHY···NCH···HM complexes are studied by ab initio calculations, where X = F, Cl, Br, Y = S, Se, and M = Li, Na, BeH, MgH. The ab initio calculations are carried out at the MP2/aug-cc-pVTZ level. The results indicate that all optimized Y···N and H···H binding distances in the ternary complexes are smaller than the corresponding values in the binary systems. The calculated cooperative energies (Ecoop) are between −0.20 kcal/mol in BrHS···NCH···HBeH and −3.29 kcal/mol in FHSe···NCH···HNa. For a given Y and M, the estimated Ecoop values increase as X = F > Cl > Br. In addition, the selenium-bonded complexes exibit larger Ecoop values than those of the sulfur-bonded counterparts. The cooperativity between Y···N and H···H interactions is further analyzed by quantum theory of atoms in molecules and natural bond orbital methods. Cooperative effects make an increase in the J(Y–N) and J(H–H) spin–spin coupling constants of the ternary complexes with respect to the binary systems.

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.

An intensified π-hole in beryllium-doped boron nitride meshes: its determinant role in CO2 conversion into hydrocarbon fuels

DFT investigations on beryllium-doped boron nitride meshes or sheets (BNs) predict the existence of a very reactive kind of novel material capable of spontaneously reducing the first hydrogenation step in the CO₂ conversion mechanism.

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.

Chalcogen bonds in complexes of SOXY (X, Y = F, Cl) with nitrogen bases

SOF2, SOFCl, and SOCl2 were each paired with a series of N bases. The potential energy surface of the binary complexes were characterized by MP2 calculations with double and triple-ξ basis sets, extrapolated to complete sets. The most stable configurations contained a S···N chalcogen bond with interaction energies as high as 6.8 kcal/mol. These structures are stabilized by a Nlp → σ*(S-Z) electron transfer (Z = O, F, Cl), complemented by Coulombic attraction of N to the σ-hole opposite the Z atom. N···S-F and N···S-Cl chalcogen bonds are stronger than N···S═O interactions. Formation of each chalcogen bond elongates all of the internal covalent bonds within SOXY, especially the S-Cl bond. Halogen-bonded (N···Cl-S) complexes were also observed, but these are more weakly bound, by less than 3 kcal/mol.