Chalcogen bonds: Hierarchical ab initio benchmark and density functional theory performance study

Dynamic Combinatorial Chemistry of Ditellurides

Dynamic combinatorial chemistry enables the exploration of reversible chemical exchange reactions. Here, we report studies on the dynamic Te─Te bond exchange of chemically and structurally diverse ditellurides using 125Te NMR spectroscopy. The results revealed rapid and unbiased equilibration, forming a consistent 1:1 distribution of reactants and products across all chemical systems. Extending to multi-component mixtures, we observed robust equilibrium formation even in complex systems involving multiple ditelluride species that can be easily distinguished by 125Te NMR spectroscopy. In line with the NMR analyses showing that the Te─Te bond exchange efficiently proceeds in various solvents, at very low temperatures, in the absence of visible light, and in the presence of radical scavengers, quantum chemical analyses revealed that the ditelluride-ditelluride exchange likely proceeds via a non-synchronous concerted mechanism and not via an alternative radical mechanism. Integrated NMR spectroscopic and quantum chemical analyses also showed that the efficiency of the chalcogen-chalcogen exchange involving ditellurides follows the trend Te-Te > Te-Se > Te-S. Overall, the results indicate that ditellurides are excellent dynamic systems that may find applications across molecular sciences.

Radical Addition Reactions: Hierarchical Ab Initio Benchmark and DFT Performance Study

We performed a hierarchical ab initio benchmark study of the gas-phase radical addition reactions of X⋅+C2H2 and X⋅+C2H4 (X⋅ = CH3⋅, NH2⋅, OH⋅, SH⋅). The hierarchical series of ab initio methods (HF, MP2, CCSD, CCSD(T)) were paired with a hierarchal series of Dunning basis sets with and without diffuse functions ((aug)-cc-pVDZ, (aug)-cc-pVTZ, (aug)-cc-pVQZ). The HF ground-state wavefunctions were transformed into quasi-restricted orbital (QRO) reference wavefunctions to address spin contamination. Following extrapolation to the CBS limit, the energies from our highest- QRO-CCSD(T)/CBS+ level converged within 0.0-3.4 kcal mol-1 and 0.0-1.0 kcal mol-1 concerning the ab initio method and basis set, respectively. Our QRO-CCSD(T)/CBS+ reference data was used to evaluate the performance of 98 density functional theory (DFT) approximations. The MAE of the best functionals for reaction barriers and energies were: OLYP (1.9 kcal mol-1), BMK (1.0 kcal mol-1), M06-2X (0.9 kcal mol-1), MN12-SX (0.8 kcal mol-1) and CAM-B3LYP (0.8 kcal mol-1). These functionals also accurately reproduce key geometrical parameters of the stationary points within an average 2 % deviation from the reference QRO-CCSD(T)/cc-pVTZ level.

A Benchmark Study of DFT-Computed p-Block Element Lewis Pair Formation Enthalpies Against Experimental Calorimetric Data

The quantification of Lewis acidity is of fundamental and applied importance in chemistry. While the computed fluoride ion affinity (FIA) is the most widely accepted thermodynamic metric, only sparse experimental values exist. Accordingly, a benchmark of methods for computing Lewis pair formation enthalpies, also with a broader set of Lewis bases against experimental data, is missing. Herein, we evaluate different density functionals against a set of 112 experimentally determined Lewis acid/base binding enthalpies and gauge influences such as solvation correction in structure optimization. From that, we can recommend r2SCAN-3c for robust quantification of this omnipresent interaction.

Hypervalent Chalcogen Bonds Catalysis on the Intramolecular Aza-Michael Reaction of Aminochalcone: Catalytic Performance and Chalcogen Bond Properties

Chalcogen bond (ChB) catalysis, as a new type in the field of non-covalent bond catalysis, has become a hot research topic in the field of organocatalysis in recent years. In the present work, we investigated the catalytic performance of a series of hypervalent ChB catalysis based on the intramolecular Aza-Michael reaction of aminochalcone. The reaction includes the carbon-nitrogen bond coupling step (key step) and the proton transfer step. The catalytic performance of mono-dentate pentafluorophenyl chalcogen bond donor ChB1 was comparable to that of bis-dentate chalcogen bond donor ChB4, and stronger than that of mono-dentate chalcogen bond donors ChB2 and ChB3. The formation of the chalcogen bond between the catalyst and the carbonyl oxygen atom of the reactant, causing the charge rearrangement of the reactant and C(1) charge of the -C-Ph group to become more positive, thereby the ChB catalysis promoted the nucleophile reaction. The electron density of the chalcogen bond of the pre-complex, the most positive electrostatic potentials of the catalyst, and the NPA charge of the key atom are proportional to the Gibbs energy barrier of the C-N bond coupling process, which provides an idea to predict the catalytic activity of the ChB catalysis.

Chalcogen-Bond-Assisted Formation of the N→C Dative Bonds in the Complexes between Chalcogenadiazoles/Chalcogenatriazoles and Fullerene C60

The existence of the N→C dative bonds in the complexes between N-containing molecules and fullerenes have been verified both theoretically and experimentally. However, finding stable N→C dative bonds is still a highly challenging task. In this work, we investigated computationally the N→C dative bonds in the complexes formed by fullerene C60 with 1,2,5-chalcogenadiazoles, 2,1,3-benzochalcogenadiazoles, and 1,2,4,5-chalcogenatriazoles, respectively. It was found that the N→C dative bonds are formed along with the formation of the N-Ch···C (Ch = S, Se, Te) chalcogen bonds. In the gas phase, from S-containing complexes through Se-containing complexes to Te-containing complexes, the intrinsic interaction energies become more and more negative, which indicates that the N-Ch···C chalcogen bonds can facilitate the formation of the N→C dative bonds. The intrinsic interaction energies are compensated by the large deformation energy of fullerene C60. The total interaction energies of Te-containing complexes are negative, while both total interaction energies of the S-containing complexes and Se-containing complexes are positive. This means that the N→C dative bonds in the Te-containing complexes are more easily observed in experiments in comparison with those in the S-containing complexes and Se-containing complexes. This study provides a new theoretical perspective on the experimental observation of the N→C dative bonds in complexes involving fullerenes. Further, the formation of stable N→C dative bonds in the complexes involving fullerenes can significantly change the properties of fullerenes, which will greatly simulate and expand the application range of fullerenes.

Exploiting π and Chalcogen Interactions for the β-Selective Glycosylation of Indoles through Glycal Conformational Distortion

Harnessing unconventional noncovalent interactions (NCIs) is emerging as a formidable synthetic approach in difficult-to-access glycosidic chemical space. C-Glycosylation, in particular, has gained a flurry of recent attention. However, most reported methods are restricted to the relatively facile access to α-C-glycosides. Herein, we disclose a β-stereoselective glycosylation of indoles by employing a phosphonoselenide catalyst. The robustness of this protocol is exemplified by its amenability for reaction at both the indolyl C- and N- reactivity sites. In contrast to previous reports, in which the chalcogens were solely involved in Lewis acidic activation, our mechanistic investigation unraveled that the often neglected flanking aromatic substituents of phosphonoselenides can substantially contribute to catalysis by engaging in π-interactions. Computations and NMR spectroscopy indicated that the chalcogenic and aromatic components of the catalyst can be collectively exploited to foster conformational distortion of the glycal away from the usual half-chair to the boat conformation, which liberates the convex β-face for nucleophilic attack.

Computational Insight into the Nature and Strength of the π-Hole Type Chalcogen∙∙∙Chalcogen Interactions in the XO2∙∙∙CH3YCH3 Complexes (X = S, Se, Te; Y = O, S, Se, Te)

In recent years, the non-covalent interactions between chalcogen centers have aroused substantial research interest because of their potential applications in organocatalysis, materials science, drug design, biological systems, crystal engineering, and molecular recognition. However, studies on π-hole-type chalcogen∙∙∙chalcogen interactions are scarcely reported in the literature. Herein, the π-hole-type intermolecular chalcogen∙∙∙chalcogen interactions in the model complexes formed between XO2 (X = S, Se, Te) and CH3YCH3 (Y = O, S, Se, Te) were systematically studied by using quantum chemical computations. The model complexes are stabilized via one primary X∙∙∙Y chalcogen bond (ChB) and the secondary C-H∙∙∙O hydrogen bonds. The binding energies of the studied complexes are in the range of -21.6~-60.4 kJ/mol. The X∙∙∙Y distances are significantly smaller than the sum of the van der Waals radii of the corresponding two atoms. The X∙∙∙Y ChBs in all the studied complexes except for the SO2∙∙∙CH3OCH3 complex are strong in strength and display a partial covalent character revealed by conducting the quantum theory of atoms in molecules (QTAIM), a non-covalent interaction plot (NCIplot), and natural bond orbital (NBO) analyses. The symmetry-adapted perturbation theory (SAPT) analysis discloses that the X∙∙∙Y ChBs are primarily dominated by the electrostatic component.

Regularized and Opposite Spin-Scaled Functionals from Møller-Plesset Adiabatic Connection─Higher Accuracy at Lower Cost

Noncovalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller-Plesset adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation energies. This combination yields cosκos-SPL2, which exhibits superior accuracy for NCIs compared to any of the individual strategies. With the N4 formal scaling, cosκos-SPL2 is competitive or often outperforms more expensive dispersion-corrected double hybrids for NCIs. The accuracy of cosκos-SPL2 particularly shines for anionic halogen bonded complexes, where it surpasses standard dispersion-corrected DFT by a factor of 3 to 5.

The Key Role of Chalcogenurane Intermediates in the Reduction Mechanism of Sulfoxides and Selenoxides by Thiols Explored In Silico

Sulfoxides and selenoxides oxidize thiols to disulfides while being reduced back to sulfides and selenides. While the reduction mechanism of sulfoxides to sulfides has been thoroughly explored experimentally as well as computationally, less attention has been devoted to the heavier selenoxides. In this work, we explore the reductive mechanism of dimethyl selenoxide, as an archetypal selenoxide and, for the sake of comparison, the reductive mechanism of dimethyl sulfoxide to gain insight into the role of the chalcogen on the reaction substrate. Particular attention is devoted to the key role of sulfurane and selenurane intermediates. Moreover, the capacity of these system to oxidize selenols rather than thiols, leading to the formation of selenyl sulfide bridges, is explored in silico. Notably, this analysis provides molecular insight into the role of selenocysteine in methionine sulfoxide reductase selenoenzyme. The activation strain model of chemical reactivity is employed in the studied reactions as an intuitive tool to bridge the computationally predicted effect of the chalcogen on the chalcogenoxide as well as on the chalcogenol.

Hydrogen Bonding versus Halogen Bonding: Spectroscopic Investigation of Gas-Phase Complexes Involving Bromide and Chloromethanes

Hydrogen bonding and halogen bonding are important non-covalent interactions that are known to occur in large molecular systems, such as in proteins and crystal structures. Although these interactions are important on a large scale, studying hydrogen and halogen bonding in small, gas-phase chemical species allows for the binding strengths to be determined and compared at a fundamental level. In this study, anion photoelectron spectra are presented for the gas-phase complexes involving bromide and the four chloromethanes, CH₃ Cl, CH₂ Cl₂ , CHCl₃ , and CCl₄ . The stabilisation energy and electron binding energy associated with each complex are determined experimentally, and the spectra are rationalised by high-level CCSD(T) calculations to determine the non-covalent interactions binding the complexes. These calculations involve nucleophilic bromide and electrophilic bromine interactions with chloromethanes, where the binding motifs, dissociation energies and vertical detachment energies are compared in terms of hydrogen bonding and halogen bonding.

Involvement of Arsenic Atom of AsF3 in Five Pnicogen Bonds: Differences between X-ray Structure and Theoretical Models

Bonding within the AsF3 crystal is analyzed via quantum chemical methods so as to identify and quantify the pnicogen bonds that are present. The structure of a finite crystal segment containing nine molecules is compared with that of a fully optimized cluster of the same size. The geometries are qualitatively different, with a much larger binding energy within the optimized nonamer. Although the total interaction energy of a central unit with the remaining peripheral molecules is comparable for the two structures, the binding of the peripherals with one another is far larger in the optimized cluster. This distinction of much stronger total binding within the optimized cluster is not limited to the nonamer but repeats itself for smaller aggregates as well. The average binding energy of the cluster rises quickly with size, asymptotically approaching a value nearly triple that of the dimer.

Simulating chalcogen bonding using molecular mechanics: a pseudoatom approach to model ebselen

The organoselenium compound ebselen has recently been investigated as a treatment for COVID-19; however, efforts to model ebselen in silico have been hampered by the lack of an efficient and accurate method to assess its binding to biological macromolecules. We present here a Generalized Amber Force Field modification which incorporates classical parameters for the selenium atom in ebselen, as well as a positively charged pseudoatom to simulate the σ-hole, a quantum mechanical phenomenon that dominates the chemistry of ebselen. Our approach is justified using an energy decomposition analysis of a number of density functional theory–optimized structures, which shows that the σ-hole interaction is primarily electrostatic in origin. Finally, our model is verified by conducting molecular dynamics simulations on a number of simple complexes, as well as the clinically relevant enzyme SOD1 (superoxide dismutase), which is known to bind to ebselen.

Graphical Abstract

A Highly Efficient Neutral Anion Receptor in Polar Environments by Synergy of Anion-π Interactions and Hydrogen Bonding

Herein, it is shown how anion recognition in highly polar solvents by neutral metal-free receptors is feasible when multiple hydrogen bonding and anion-π interactions are suitably combined. A neutral aromatic molecular tweezer functionalized with azo groups is shown to merge these two kinds of interactions in a unique system and its efficiency as an anion catcher in water is evaluated using first-principles quantum methods. Theoretical calculations unequivocally prove the high thermodynamic stability in water of a model anion, bromide, captured within the tweezer's cavity. Thus, static calculations indicate anion-tweezer interaction energies within the range of covalent or ionic bonds and stability constants in water of more than 10 orders of magnitude. First-principles molecular dynamics calculations also corroborate the stability through the time of the anion-tweezer complex in water. It shows that the anion is always found within the tweezer's cavity due to the combination of the tweezer-anion interactions plus a hydrogen bond between the anion and a water molecule that is inside the tweezer's cavity.

Chalcogen S∙∙∙S Bonding in Supramolecular Assemblies of Cadmium(II) Coordination Polymers with Pyridine-Based Ligands

Two cadmium(II) coordination polymers, with thiocyanate and pyridine-based ligands e.g., 3-acetamidopyridine (3-Acpy) and niazid (nicotinic acid hydrazide, nia), namely one-dimensional {[Cd(SCN)2(3-Acpy)]}n (1) and two-dimensional {[Cd(SCN)2(nia)]}n (2), are prepared in the mixture of water and ethanol. The adjacent cadmium(II) ions in 1 are bridged by two N,S-thiocyanate ions and an N,O-bridging 3-Acpy molecule, forming infinite one-dimensional polymeric chains, which are assembled by the intermolecular N–H∙∙∙S hydrogen bonds in one direction and by the intermolecular S∙∙∙S chalcogen bonds in another direction. Within the coordination network of 2, the adjacent cadmium(II) ions are bridged by N,S-thiocyanate ions in one direction and by N,O,N’-chelating and bridging nia molecules in another direction. The coordination networks of 2 are assembled by the intermolecular N–H∙∙∙S and N–H∙∙∙N hydrogen bonds and S∙∙∙S chalcogen bonds. Being the only supramolecular interactions responsible for assembling the polymer chains of 1 in the particular direction, the chalcogen S∙∙∙S bonds are more significant in the structure of 1, whilst the chalcogen S∙∙∙S bonds which act in cooperation with the N–H∙∙∙S and N–H∙∙∙N hydrogen bonds are of less significance in the structure of 2.

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.

Chalcogen bonds: Hierarchical ab initio benchmark and density functional theory performance study

We have performed a hierarchical ab initio benchmark and DFT performance study of D2 Ch•••A- chalcogen bonds (Ch = S, Se; D, A = F, Cl). The ab initio benchmark study is based on a series of ZORA-relativistic quantum chemical methods [HF, MP2, CCSD, CCSD(T)], and all-electron relativistically contracted variants of Karlsruhe basis sets (ZORA-def2-SVP, ZORA-def2-TZVPP, ZORA-def2-QZVPP) with and without diffuse functions. The highest-level ZORA-CCSD(T)/ma-ZORA-def2-QZVPP counterpoise-corrected complexation energies (ΔECPC ) are converged within 1.1-3.4 kcal mol-1 and 1.5-3.1 kcal mol-1 with respect to the method and basis set, respectively. Next, we used the ZORA-CCSD(T)/ma-ZORA-def2-QZVPP (ΔECPC ) as reference data for analyzing the performance of 13 different ZORA-relativistic DFT approaches in combination with the Slater-type QZ4P basis set. We find that the three-best performing functionals are M06-2X, B3LYP, and M06, with mean absolute errors (MAE) of 4.1, 4.2, and 4.3 kcal mol-1 , respectively. The MAE for BLYP-D3(BJ) and PBE amount to 8.5 and 9.3 kcal mol-1 , respectively.

A Quantitative Molecular Orbital Perspective of the Chalcogen Bond

We have quantum chemically analyzed the structure and stability of archetypal chalcogen-bonded model complexes D2 Ch⋅⋅⋅A- (Ch = O, S, Se, Te; D, A = F, Cl, Br) using relativistic density functional theory at ZORA-M06/QZ4P. Our purpose is twofold: (i) to compute accurate trends in chalcogen-bond strength based on a set of consistent data; and (ii) to rationalize these trends in terms of detailed analyses of the bonding mechanism based on quantitative Kohn-Sham molecular orbital (KS-MO) theory in combination with a canonical energy decomposition analysis (EDA). At odds with the commonly accepted view of chalcogen bonding as a predominantly electrostatic phenomenon, we find that chalcogen bonds, just as hydrogen and halogen bonds, have a significant covalent character stemming from strong HOMO-LUMO interactions. Besides providing significantly to the bond strength, these orbital interactions are also manifested by the structural distortions they induce as well as the associated charge transfer from A- to D2 Ch.

Noncovalent Bonds through Sigma and Pi-Hole Located on the Same Molecule. Guiding Principles and Comparisons

Over the last years, scientific interest in noncovalent interactions based on the presence of electron-depleted regions called σ-holes or π-holes has markedly accelerated. Their high directionality and strength, comparable to hydrogen bonds, has been documented in many fields of modern chemistry. The current review gathers and digests recent results concerning these bonds, with a focus on those systems where both σ and π-holes are present on the same molecule. The underlying principles guiding the bonding in both sorts of interactions are discussed, and the trends that emerge from recent work offer a guide as to how one might design systems that allow multiple noncovalent bonds to occur simultaneously, or that prefer one bond type over another.