Exploring the influence of (n - 1)d subvalence correlation and of spin-orbit coupling on chalcogen bonding

This article presents a comprehensive computational investigation into chalcogen bonding interactions, focusing specifically on elucidating the role of subvalence (n - 1)d and (n - 1)sp correlation. The incorporation of inner-shell (n - 1)d correlation leads to a decrease in interaction energies for chalcogen-bonded systems (at least those studied herein), contradicting the observations regarding halogen bonding documented by Kesharwani et al. in J. Phys. Chem. A, 2018, 122 (8), 2184-2197. The significance of (n - 1)sp subvalence correlation appears to be lower by an order of magnitude. Notably, among the various components of interaction energies computed at the PNO-LCCSD(T) or DF-CCSD levels, we identify the PNO-LMP2 or DF-MP2 component of the (n - 1)d correlation as predominant. Furthermore, we delve into the impact of second-order spin-orbit coupling (SOC2) on these interactions. The SOC2 effects appear to be less significant than the (n - 1)d correlation; however, they remain non-trivial, particularly for Te complexes. For the Se complexes, SOC2 is much less important. Generally, SOC2 stabilizes monomers more than dimers, resulting in reduced binding of the latter. Notably, at equilibrium and stretched geometries, SOC2 and (n - 1)d destabilize the complex; however, at compressed geometries, they exhibit opposing effects, with (n - 1)d becoming stabilizing.

Exploring the Effects of Se Basicity on a Te···Se Interaction Supported by a Rigid Indazolium Backbone

With an interest in chalcogen bonding, we use a rigid indazolium backbone to install a formally zero-valent Se center next to a divalent Te center, allowing us to investigate the effects of oxidation of the Se center on the observed Te···Se interaction. Through spectroscopic and computational comparison of the Se(0) species with its Se(II) counterpart and their monochalcogen analogues, we experimentally and computationally investigate the effect of modulating Se basicity on the resulting Te···Se interaction. Comparison with well-studied naphthalene and acenaphthene variants indicates that the increased basicity of the Se(0) center allows for a comparably strong Te···Se interaction despite longer peri distances and a larger splay angle. Finally, our study illuminates the potential non-innocence of cationic organic substituents in chalcogen-bonding catalysis of the transfer hydrogenation of quinolines.

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.

Nature of the E⋯E′ interactions (E, E′ = O, S, Se, and Te) at naphthalene 1,8-positions with fine details of the structures: experimental and theoretical investigations

The nature of E⋯E′ in 1-RE–8-R′E′C₁₀H₆ (E/E′ = O, S, Se and Te) is clarified with the QTAIM approach and NBO analysis, after structural determinations.

The nature of G⋯E–Y σ(3c–4e) in o-MenGCH2C6H4EY (MenG = Me2N and MeE; E = O, S, Se and Te; Y = F, Cl, Br, EMe and Me) with contributions from CT and compliance constants in noncovalent G⋯E interactions

The intrinsic dynamic and static nature of G–*–E–*–Y σ(3c–4e) interactions was elucidated with the quantum theory of atoms in molecules dual functional analysis (QTAIM-DFA), employing o-Me_nGCH₂C₆H₄EY (Me_nG = Me₂N and MeE; E = O, S, Se and Te; Y = F, Cl, Br, I, EMe and Me). Asterisks (*) are employed to emphasize the existence of bond critical points (BCPs) on the bond paths (BPs), corresponding to the interactions in question. Data from the fully optimized structure correspond to the static nature of interactions. The dynamic nature is called the intrinsic dynamic nature if the perturbed structures are generated using the coordinates derived from the compliance constants. Basis sets of the Sapporo-TZP type with diffusion functions are employed for the heteroatoms at the MP2 level. The noncovalent G–*–E interactions in GEY σ(3c–4e) are predicted to demonstrate van der Waals bonding to CT-TBP (trigonal bipyramidal adduct formation through charge transfer) nature, while the E–*–Y bonds have the covalent nature. Some E–F bonds show strong ionic character when G–*–E is predicted to be stronger than E–*–Y. The contributions of the CT terms to the G–*–E interactions, evaluated with NBO, are discussed in relation to the predicted nature. The E(2) values based on NBO are strongly correlated to the compliance constants for the G–*–E interactions if suitably treated separately.

The nature of E⋯E′ in 1-RECH₂-2-R′E′C₆H₄ (E/E′ = O, S, Se and Te) is clarified with QTAIM approach and NBO analysis, after structural determinations.

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.

Intermolecular interactions of oligothienoacenes: Do S⋯S interactions positively contribute to crystal structures of sulfur-containing aromatic molecules?

Intermolecular interactions in the crystals of tetra- and penta-thienoacene were studied using ab initio molecular orbital calculations for evaluating the magnitude of characteristic S⋯S interactions with great attention paid to their origin. The interactions between the π-stacked neighboring molecules are significantly greater than those between the neighboring molecules exhibiting the S⋯S contact, although it has sometimes been claimed that the S⋯S interactions play important roles in adjusting the molecular arrangement of sulfur-containing polycyclic aromatic molecules in the crystals owing to short S⋯S contacts. The coupled cluster calculations with single and double substitutions with noniterative triple excitation interaction energies at the basis set limit estimated for the π-stacked and S⋯S contacted neighboring molecules in the tetrathienoacene crystal are -11.17 and -4.27 kcal/mol, respectively. Those for π-stacked molecules in the pentathienoacene crystal is -14.38 kcal/mol, while those for S⋯S contacted molecules are -7.02 and -6.74 kcal/mol. The dispersion interaction is the major source of the attraction between the π-stacked and S⋯S contacted molecules, while the orbital-orbital interactions are repulsive: The orbital-orbital interactions, which are significant for charge carrier transport properties, are not much more than the results of the short S⋯S contact caused by the strong dispersion interactions. Besides, the intermolecular interaction energy calculated for a trithienoacene dimer has strong orientation dependence.

Isolation and reversible dimerization of a selenium-selenium three-electron σ-bond

Three-electron σ-bonding that was proposed by Linus Pauling in 1931 has been recognized as important in intermediates encountered in many areas. A number of three-electron bonding systems have been spectroscopically investigated in the gas phase, solution and solid matrix. However, X-ray diffraction studies have only been possible on simple noble gas dimer Xe∴Xe and cyclic framework-constrained N∴N radical cations. Here, we show that a diselena species modified with a naphthalene scaffold can undergo one-electron oxidation using a large and weakly coordinating anion, to afford a room-temperature-stable radical cation containing a Se∴Se three-electron σ-bond. When a small anion is used, a reversible dimerization with phase and marked colour changes is observed: radical cation in solution (blue) but diamagnetic dimer in the solid state (brown). These findings suggest that more examples of three-electron σ-bonds may be stabilized and isolated by using naphthalene scaffolds together with large and weakly coordinating anions.

Probing interactions through space using spin–spin coupling

A series of eight 5-(TeAr)-6-(SePh)acenaphthenes (Ar = aryl) were prepared and structurally characterised by X-ray crystallography, solution and solid-state NMR spectroscopy and DFT/B3LYP calculations.

Electrochemically informed synthesis: oxidation versus coordination of 5,6-bis(phenylchalcogeno)acenaphthenes

Chalcogen dications: Facile synthesis of E--E bonded dications can be readily achieved. Radical cations are identified as the intermediates.

Investigating silver coordination to mixed chalcogen ligands

Six silver(I) coordination complexes have been prepared and structurally characterised. Mixed chalcogen-donor acenaphthene ligands L1-L3 [Acenap(EPh)(E'Ph)] (Acenap = acenaphthene-5,6-diyl; E/E' = S, Se, Te) were independently treated with silver(I) salts (AgBF₄/AgOTf). In order to keep the number of variables to a minimum, all reactions were carried out using a 1:1 ratio of Ag/L and run in dichloromethane. The nature of the donor atoms, the coordinating ability of the respective counter-anion and the type of solvent used in recrystallisation, all affect the structural architecture of the final silver(I) complex, generating monomeric, silver(I) complexes {[AgBF₄(L)₂] (1 L = L1; 2 L = L2; 3 L = L3), [AgOTf(L)₃] (4 L = L1; 5 L = L3), [AgBF₄(L)₃] (2a L = L1; 3a L = L3)} and a 1D polymeric chain {[AgOTf(L3)](n) 6}. The organic acenaphthene ligands L1-L3 adopt a number of ligation modes (bis-monodentate μ₂-η²-bridging, quasi-chelating combining monodentate and η⁶-E(phenyl)-Ag(I) and classical monodentate coordination) with the central silver atom at the centre of a tetrahedral or trigonal planar coordination geometry in each case. The importance of weak interactions in the formation of metal-organic structures is also highlighted by the number of short non-covalent contacts present within each complex.

Onset of three-centre, four-electron bonding in peri-substituted acenaphthenes: a structural and computational investigation

Two series of sterically crowded peri-substituted acenaphthenes have been prepared, containing mixed halogen-chalcogen functionalities at the 5,6-positions in A1-A6 (Acenap[X][EPh] (Acenap = acenaphthene-5,6-diyl; X = Br, I; E = S, Se, Te) and chalcogen-chalcogen moieties in A7-A12 (Acenap[EPh][E'Ph] (Acenap = acenaphthene-5,6-diyl; E/E' = S, Se, Te). The related dihalide compounds A13-A16 Acenap[XX'] (XX' = BrBr, II, IBr, ClCl) have also been prepared. Distortion of the acenaphthene framework away from the ideal was studied as a function of the steric bulk of the interacting halogen and chalcogen atoms occupying the peri-positions. The acenaphthene series experiences a general increase in peri-separation for molecules accommodating heavier congeners and maps the trends observed previously for the analogous naphthalene compounds N1-N12 (Nap[X][EPh], Nap[EPh][E'Ph] (X = Br, I; E/E' = S, Se, Te). The conformation of the aromatic ring systems and subsequent location of p-type lone-pairs dominates the geometry of the peri-region. The differences in peri-separations observed for compounds adopting differing conformations of the peri-substituted phenyl group can be correlated to the ability of the frontier orbitals of the halogen or chalcogen atoms to take part in attractive or repulsive interactions. Density-functional studies have confirmed these interactions and suggested the onset of formation of three-centre, four-electron bonding under appropriate geometric conditions.

A bidentate Lewis acid with a telluronium ion as an anion-binding site

The search for receptors that can selectively capture small and potentially toxic anions in protic media has sparked a renewed interest in the synthesis and anion-binding properties of polydentate Lewis acids. Seeking new paradigms to enhance the anion affinities of such systems, we synthesized a bidentate Lewis acid that contains a boryl and a telluronium moiety as Lewis acidic sites. Anion-complexation studies indicate that this telluronium borane displays a high affinity for fluoride in methanol. Structural and computational studies show that the unusual fluoride affinity of this bidentate telluronium borane can be correlated with the formation of a B-F → Te chelate motif supported by a strong lone-pair(F) → σ*(Te-C) donor-acceptor interaction. These results, which illustrate the viability of heavier chalcogenium centres as anion-binding sites, allow us to introduce a novel strategy for the design of polydentate Lewis acids with enhanced anion affinities.

Synthetic and structural studies of 1,8-chalcogen naphthalene derivatives

Four novel 1,8-disubstituted naphthalene derivatives 4-7 that contain chalcogen atoms occupying the peri positions have been prepared and fully characterised by using X-ray crystallography, multinuclear NMR spectroscopy, IR spectroscopy and MS. Molecular distortion due to noncovalent substituent interactions was studied as a function of the bulk of the interacting chalcogen atoms and the size and nature of the alkyl group attached to them. X-ray data for 4-7 was compared to the series of known 1,8-bis(phenylchalcogeno)naphthalenes 1-3, which were themselves prepared from novel synthetic routes. A general increase in the E...E' distance was observed for molecules containing bulkier atoms at the peri positions. The decreased S...S distance from phenyl-1 and ethyl-4 analogues is ascribed to a weaker chalcogen lone pair-lone pair repulsion acting in the ethyl analogue due to the presence of two equatorial S(naphthyl) ring conformations. Two novel peri-substituted naphthalene sulfoxides of 1, Nap(O=SPh)(SPh) 8 and Nap(O=SPh)(2) 9, which contain different valence states of sulfur, were prepared and fully characterised by using X-ray crystallography and multinuclear NMR spectroscopy, IR spectroscopy and MS. Molecular structures were analysed by using naphthalene ring torsions, peri-atom displacement, splay angle magnitude, S...S interactions, aromatic ring orientations and quasi-linear O=S...S arrangements. The axial S(naphthyl) rings in 8 and 9 are unfavourable for S...S contacts due to stronger chalcogen lone pair-lone pair repulsion. Although quasi-linear O=S...S alignments suggest attractive interaction is conceivable, analysis of the B3LYP wavefunctions affords no evidence for direct bonding interactions between the S atoms.

Evidence for effective p(Z)-pi(Ar) conjugations (Z = S, Se, and Te, as well as Z = O) in 9-(arylchalcogenyl)triptycenes: experimental and theoretical investigations

The p(Z)-pi(Ar) conjugations must operate fully in the ground states of 9-(arylchalcogenyl)triptycenes (p-YC(6)H(4)ZTpc:1 (Z = O), 2 (Z = S), 3 (Z = Se), and 4 (Z = Te)), where the p-YC(6)H(4) group is placed in the bisected area between two phenyl planes of the triptycyl group with the parallel orientation. The ground-state geometries, which we call (A: pl), are confirmed by X-ray analysis. However, the conjugations never operate in the transition states between (A: pl) and/or the topomeric structures (A': pl'), where the Z-C(Tpc) bond is perpendicular to the plane. The site-exchange processes correlate to the conjugations. Temperature-dependent (1)H NMR spectra are analyzed for 2 and 3 to demonstrate the effective p(Z)-pi(Ar) conjugations. The activation energies for the interconversion between (A: pl) and (A': pl') (GR: gear process) were obtained for 2 (DeltaG(GR)(2)) and 3 (DeltaG(GR)(3)). DeltaG(GR)(3) correlate well with DeltaG(GR)(2), and DeltaG(GR)(2) are well analyzed by the Hammett-type dual parameters. DeltaG(GR)(2) and DeltaG(GR)(3) are demonstrated to be controlled by the resonance interaction of the p(Z)-pi(C(6)H(4))-p(Y) conjugations. QC calculations are performed on the ground and exited states of 1-4, which clarify the effective p(Z)-pi(C(6)H(4))-p(Y) conjugations for Z of heavier atoms.