Chalcogen bonding of two ligands to hypervalent YF4 (Y = S, Se, Te, Po)

σ-Hole Site-Based Interactions within Hypervalent Pnicogen, Halogen, and Aerogen-Bearing Molecules with Lewis Bases: A Comparative Study

σ-Hole site-based interactions in the trigonal bipyramidal geometrical structure of hypervalent pnicogen, halogen, and aerogen-bearing molecules with pyridine and NCH Lewis bases (LBs) were comparatively examined. In this respect, the ZF5···, XF3O2···, and AeF2O3···LB complexes (where Z = As, Sb; X = Br, I; Ae = Kr, Xe; and LB = pyridine and NCH) were investigated. The electrostatic potential (EP) analysis affirmations outlined the occurrence of σ-holes on the systems under consideration with disparate magnitudes that increased according to the following order: AeF2O3 < XF3O2 < ZF5. In line with EP outcomes, the proficiency of σ-hole site-based interactions increased as the atomic size of the central atom increased with a higher favorability for the pyridine-based complexes over NCH-based ones. The interaction energy showed the most favorable negative values of -35.97, -44.53, and -56.06 kcal/mol for the XeF2O3···, IF3O2···, and SbF5···pyridine complexes, respectively. The preferentiality pattern of the studied interactions could be explained as a consequence of (i) the dramatic rearrangement of ZF5 molecules from the trigonal bipyramid geometry to the square pyramidal one, (ii) the significant and tiny deformation energy in the case of the interaction of XF3O2 molecules with pyridine and NCH, respectively, and (iii) the absence of geometrical deformation within the AeF2O3···pyridine and ···NCH complexes other than the XeF2O3···pyridine one. Quantum theory of atoms in molecules and noncovalent interaction index findings reveal the partially covalent nature of most of the investigated interactions. Symmetry-adapted perturbation theory affirmations declared that the electrostatic component was the driving force beyond the occurrence of the considered interactions. The obtained findings will help in improving our understanding of the effect of geometrical deformation on intermolecular interactions.

Hole interactions of aerogen oxides with Lewis bases: an insight into σ-hole and lone-pair-hole interactions

σ-Hole and lone-pair (lp)-hole interactions of aerogen oxides with Lewis bases (LB) were comparatively inspected in terms of quantum mechanics calculations. The ZOn ⋯ LB complexes (where Z = Kr and Xe, n = 1, 2, 3 and 4, and LB = NH3 and NCH) showed favourable negative interaction energies. The complexation features were explained in light of σ-hole and lp-hole interactions within optimum distances lower than the sum of the respective van der Waals radii. The emerging findings outlined that σ-hole interaction energies generally enhanced according to the following order: KrO4 ⋯ < KrO⋯ < KrO3⋯ < KrO2⋯LB and XeO4⋯ < XeO⋯ < XeO2⋯ < XeO3⋯LB complexes with values ranging from -2.23 to -12.84 kcal mol-1. Lp-hole interactions with values up to -5.91 kcal mol-1 were shown. Symmetry-adapted perturbation theory findings revealed the significant contributions of electrostatic forces accounting for 50-65% of the total attractive forces within most of the ZOn⋯LB complexes. The obtained observations would be useful for the understanding of hole interactions, particularly for the aerogen oxides, with application in supramolecular chemistry and crystal engineering.

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.

Excited states of polonium(IV): electron correlation and spin-orbit coupling in the Po4+ free ion and in the bare and solvated [PoCl5]- and [PoCl6]2- complexes

Polonium (Po, Z = 84) is a main-block element with poorly known physico-chemical properties. Not much information has been firmly acquired since its discovery by Marie and Pierre Curie in 1898, especially regarding its speciation in aqueous solution and spectroscopy. In this work, we revisit the absorption properties of two complexes, [PoCl5]- and [PoCl6]2-, using quantum mechanical calculations. These complexes have the potential to exhibit a maximum absorption at 418 nm in HCl medium (for concentrations of 0.5 mol L-1 and above). Initially, we examine the electronic spectra of the Po4+ free ion and of its isoelectronic analogue, Bi3+, in the spin-orbit configuration interaction (SOCI) framework. Our findings demonstrate that the SOCI matrix should be dressed with correlated electronic energies and that the quality of the spectra is largely improved by decontracting the reference states at the complete active space plus singles (CAS + S) level. Subsequently, we investigate the absorption properties of the [PoCl5]- and [PoCl6]2- complexes in two stages. Firstly, we perform methodological tests at the MP2/def2-TZVP gas phase geometries, indicating that the decontraction of the reference states can be skipped without compromising the accuracy significantly. Secondly, we study the solution absorption properties by means of single-point calculations performed at the solvated geometries, obtained by an implicit solvation treatment or a combination of implicit and explicit solvation. Our results highlight the importance of saturating the first coordination sphere of the PoIV ion to obtain a qualitatively correct picture. Finally, we conclude that the known-for-decades 418 nm peak could be attributed to a mixture of both the [PoCl5(H2O)]- and [PoCl6]2- complexes. This finding not only aligns with the behaviour of the analogous BiIII ion under similar conditions but also potentially provides an explanation for previous discrepancies in the literature.

Sigma-Hole and Lone-Pair-Hole Site-Based Interactions of Seesaw Tetravalent Chalcogen-Bearing Molecules with Lewis Bases

For the first time, sigma (σ)- and lone-pair (lp)-hole site-based interactions of SF4 and SeF4 molecules in seesaw geometry with NH3 and FH Lewis bases were herein comparatively investigated. The obtained findings from the electrostatic potential analysis outlined the emergence of sundry holes on the molecular entity of the SF4 and SeF4 molecules, dubbed the σ- and lp-holes. The energetic viewpoint announced splendid negative binding energy values for σ-hole site-based interactions succeeded by lp-hole analogues, which were found to be -9.21 and -0.50 kcal/mol, respectively, for SeF4···NH3 complex as a case study. Conspicuously, a proper concurrence between the strength of chalcogen σ-hole site-based interactions and the chalcogen's atomic size was obtained, whereas a reverse pattern was proclaimed for the lp-hole counterparts. Further, a higher preference for the YF4···NH3 complexes with elevated negative binding energy was promulgated over the YF4···FH ones, indicating the eminent role of Lewis basicity. The indications of the quantum theory of atoms in molecules generally asserted the closed-shell nature of all the considered interactions. The observation of symmetry-adapted perturbation theory revealed the substantial contributing role of the electrostatic forces beyond the occurrence of σ-hole site-based interactions. In comparison, the dispersion forces were specified to govern the lp-hole counterparts. Such emerging findings would be a gate for the fruitful forthcoming applications of chalcogen bonding interactions in crystal engineering and biological systems.

Structure and Bonding of Halonium Compounds

The geometrical parameters and the bonding in [D···X···D]⁺ halonium compounds, where D is a Lewis base with N as the donor atom and X is Cl, Br, or I, have been investigated through a combined structural and computational study. Cambridge Structural Database (CSD) searches have revealed linear and symmetrical [D···X···D]⁺ frameworks with neutral donors. By means of density functional theory (DFT), molecular electrostatic potential (MEP), and energy decomposition analyses (EDA) calculations, we have studied the effect of various halogen atoms (X) on the [D···X···D]⁺ framework, the effect of different nitrogen-donor groups (D) attached to an iodonium cation (X = I), and the influence of the electron density alteration on the [D···I···D]⁺ halonium bond by variation of the R substituents at the N-donor upon the symmetry, strength, and nature of the interaction. The physical origin of the interaction arises from a subtle interplay between electrostatic and orbital contributions (σ-hole bond). Interaction energies as high as 45 kcal/mol suggest that halonium bonds can be exploited for the development of novel halonium transfer agents, in asymmetric halofunctionalization or as building blocks in supramolecular chemistry.

Competition Between the Two σ-Holes in the Formation of a Chalcogen Bond

A chalcogen atom Y contains two separate σ-holes when in a R₁ YR₂ molecular bonding pattern. Quantum chemical calculations consider competition between these two σ-holes to engage in a chalcogen bond (ChB) with a NH₃ base. R groups considered include F, Br, I, and tert-butyl (tBu). Also examined is the situation where the Y lies within a chalcogenazole ring, where its neighbors are C and N. Both electron-withdrawing substituents R₁ and R₂ act cooperatively to deepen the two σ-holes, but the deeper of the two holes consistently lies opposite to the more electron-withdrawing group, and is also favored to form a stronger ChB. The formation of two simultaneous ChBs in a triad requires the Y atom to act as double electron acceptor, and so anti-cooperativity weakens each bond relative to the simple dyad. This effect is such that some of the shallower σ-holes are unable to form a ChB at all when a base occupies the other site.

Geometries, interaction energies and bonding in [Po(H2O)n]4+ and [PoCln]4-n complexes

Polonium (Z = 84) is one of the rarest elements on Earth. More than a century after its discovery, its chemistry remains poorly known and even basic questions have not yet been satisfactorily addressed. In this work, we perform a systematic study of the geometries, interactions energies and bonding in basic polonium(IV) species, namely the hydrated [Po(H₂O)_n]⁴⁺ and chlorinated [PoCl_n]^4-n complexes by means of gas-phase electronic structure calculations. We show that while up to nine water molecules can fit in the first coordination sphere of the polonium(IV) ion, its coordination sphere can already be filled with eight chloride ligands. Capitalising on previous theoretical studies, a focused methodological study based on interaction energies and bond distances allows us to validate the MP2/def2-TZVP level of theory for future ground-state studies. After discussing the similarities and differences between complexes with the same number of ligands, we perform topological analyses of the MP2 electron densities in the quantum theory of atoms in molecules (QTAIM) fashion. While the water complexes display typical signatures of closed-shell interactions, we reveal large Po-Cl delocalisation indices, especially in the hypothetical [PoCl]³⁺ complex. This "enhanced" covalency opens the way for a significant spin-orbit coupling (SOC) effect on the corresponding bond distance, which has been studied using two independent approaches (i.e. one a priori and one a posteriori). We finally conclude by stressing that while the SOC may not affect much the geometries of high-coordinated polonium(IV) complexes, it should definitely not be neglected in the case of low-coordinated ones.

Various Sorts of Chalcogen Bonds Formed by an Aromatic System

The chalcogen Y atom in the aromatic ring of thiophene and its derivatives YC₄H₄ (Y = S, Se, Te) can engage in a number of different interactions with another such unit within the homodimer. Quantum calculations show that the two rings can be oriented perpendicular to one another in a T-shaped dimer in which the Y atom accepts electron density from the π-system of the other unit in a Y···π chalcogen bond (ChB). This geometry best takes advantage of attractions between the electrostatic potentials surrounding the two monomers. There are two other geometries in which the two Y atoms engage in a ChB with one another. However, instead of a simple interaction between a σ-hole on one Y and the lone pair of its neighbor, the interaction is better described as a pair of symmetrically equivalent Y···Y interactions, in which charge is transferred in both directions simultaneously, thereby effectively doubling the strength of the bond. These geometries differ from what might be expected based simply on the juxtaposition of the electrostatic potentials of the two monomers.

Principles Guiding the Square Bonding Motif Containing a Pair of Chalcogen Bonds between Chalcogenadiazoles

The bonding motif adopted by a dimer of chalcogenadiazole molecules is characterized by a pair of equivalent Ch···N chalcogen bonds. Quantum calculations show that the interaction energy is substantial, varying between 4 kcal/mol for Ch = S and 17 kcal/mol for Te. The interaction is cooperative in that the total bond strength is greater than either chalcogen bond individually. Neither the addition of a phenyl ring nor the addition of a pair of cyano substituents to the diazole ring has much influence on this binding. Removal of one N from the diazole weakens the binding, and addition of two nitrogens has little effect. The largest perturbation arises with three N atoms in each ring, for which the binding energy increases by some 25%. The ring size plays a minor role in most cases, although a near doubling of bond strength occurs if there are two N atoms present on a four-membered ring.

Square tetravalent chalcogen bonds in dimeric aggregates: a joint crystallographic survey and theoretical study

Square tetravalent chalcogen bonds were systematically investigated through a combination of crystal structure analysis and DFT calculations.

Chalcogen Bond Involving Zinc(II)/Cadmium(II) Carbonate and Its Enhancement by Spodium Bond

Carbonate MCO3 (M = Zn, Cd) can act as both Lewis acid and base to engage in a spodium bond with nitrogen-containing bases (HCN, NHCH2, and NH3) and a chalcogen bond with SeHX (X = F, Cl, OH, OCH3, NH2, and NHCH3), respectively. There is also a weak hydrogen bond in the chalcogen-bonded dyads. Both chalcogen and hydrogen bonds become stronger in the order of F > Cl > OH > OCH3 > NH2 > NHCH3. The chalcogen-bonded dyads are stabilized by a combination of electrostatic and charge transfer interactions. The interaction energy of chalcogen-bonded dyad is less than -10 kcal/mol at most cases. Furthermore, the chalcogen bond can be strengthened through coexistence with a spodium bond in N-base-MCO3-SeHX. The enhancement of chalcogen bond is primarily attributed to the charge transfer interaction. Additionally, the spodium bond is also enhanced by the chalcogen bond although the corresponding enhancing effect is small.

The role of electric field, peripheral chains, and magnetic effects on significant 1H upfield shifts of the encapsulated molecules in chalcogen-bonded capsules

The chalcogen-bonded homo-cavitand and hetero-cavitand AY+AY' capsules (Y, Y' = Se, Te), as well as their encapsulated complexes with one or two guest molecules have been studied theoretically via density functional theory (DFT), while the ¹H NMR spectra of the homo-cavitand encapsulated complexes (in ASe+ASe) have been measured experimentally. There is excellent agreement between theoretical and experimental spectra. In all cases, we found significant ¹H upfield shifts which are more intense in the ASe+ASe cage compared to the ATe+ATe and ASe+ATe cages. The non-uniform electron distribution which gives rise to an inherent electric field and a non-zero electric dipole moment of the encapsulated complexes, the induced electric field effects, the magnetic anisotropy which is enhanced due to the polarizability of chalcogen atoms, and the peripheral chains, which are responsible for the solubility of the cages, increase the upfield shifts of ¹H of the encapsulated molecules; the peripheral chains lead to an increase of the upfield shifts by up to 1.8 ppm for H of the rim and up to 1.2 ppm for the terminal H in the interior of the cage. Hence, substantial ¹H upfield chemical shifts of the guests in these capsules are consequences of (i) the enhanced aromaticity of the walls of the capsules due to the polarizability of chalcogen atoms, (ii) the induced and inherent electric field effects, and (iii) the peripheral chains.

Competition between Inter and Intramolecular Tetrel Bonds: Theoretical Studies Complemented by CSD Survey

Crystal structures document the ability of a TF₃ group (T=Si, Ge, Sn, Pb) situated on a naphthalene system to engage in an intramolecular tetrel bond (TB) with an amino group on the adjoining ring. Ab initio calculations evaluate the strength of this bond and evaluate whether it can influence the ability of the T atom to engage in a second, intermolecular TB with another nucleophile. A very strong CN^- anionic base can approach the T either along the extension of a T-C or T-F bond and form a strong TB with an interaction energy approaching 100 kcal/mol, although this bond is weakened a bit by the presence of the internal T⋅⋅⋅N bond. The much less potent NCH base engages in a correspondingly longer and weaker TB, less than 10 kcal/mol. Such an intermolecular TB is weakened by the presence of the internal TB, to the point that it only occurs for the two heavier tetrel atoms Sn and Pb.

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.

Comparison of Bifurcated Halogen with Hydrogen Bonds

Bifurcated halogen bonds are constructed with FBr and FI as Lewis acids, paired with NH₃ and NCH bases. The first type considered places two bases together with a single acid, while the reverse case of two acids sharing a single base constitutes the second type. These bifurcated systems are compared with the analogous H-bonds wherein FH serves as the acid. In most cases, a bifurcated system is energetically inferior to a single linear bond. There is a larger energetic cost to forcing the single σ-hole of an acid to interact with a pair of bases, than the other way around where two acids engage with the lone pair of a single base. In comparison to FBr and FI, the H-bonding FH acid is better able to participate in a bifurcated sharing with two bases. This behavior is traced to the properties of the monomers, in particular the specific shape of the molecular electrostatic potential, the anisotropy of the orbitals of the acid and base that interact directly with one another, and the angular extent of the total electron density of the two molecules.

Experimental and Theoretical Studies of Dimers Stabilized by Two Chalcogen Bonds in the Presence of a N···N Pnicogen Bond

The structure of the 5,6-dichloro-2,1,3-benzoselenadiazole homodimer, obtained by adding the ligand, 4,5-dichloro-o-phenylenediamine, to the methanolic solution of SeCl₄, was determined by X-ray crystallography, augmented by Fourier transform infrared, Raman, and NMR spectroscopy. The binding motif involves a pair of Se···N chalcogen bonds, with a supplementary N···N pnicogen bond. Quantum calculations provide assessments of the strengths of the individual interactions as well as their contributing factors. All together, these three bonds compose a total interaction energy between 5.4 and 16.8 kcal/mol, with the larger chalcogen atom associated with the strongest interactions. Replacement of the Se atoms by S and Te analogues allows analysis of the dependence of these forces on the nature of the chalcogen atom. Calculations also measure the importance to the binding of the presence of a second N atom on each diazole unit as well as the substituted phenyl ring to which it is fused.

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.

Versatility of the Cyano Group in Intermolecular Interactions

Several cyano groups are added to an alkane, alkene, and alkyne group so as to construct a Lewis acid molecule with a positive region of electrostatic potential in the area adjoining these substituents. Although each individual cyano group produces only a weak π-hole, when two or more such groups are properly situated, they can pool their π-holes into one much more intense positive region that is located midway between them. A NH3 base is attracted to this site, where it forms a strong noncovalent bond to the Lewis acid, amounting to as much as 13.6 kcal/mol. The precise nature of the bonding varies a bit from one complex to the next but typically contains a tetrel bond to the C atoms of the cyano groups or the C atoms of the linkage connecting the C≡N substituents. The placement of the cyano groups on a cyclic system like cyclopropane or cyclobutane has a mild weakening effect upon the binding. Although F is comparable to C≡N in terms of electron-withdrawing power, the replacement of cyano by F substituents substantially weakens the binding with NH3.

Coordination of a Central Atom by Multiple Intramolecular Pnicogen Bonds

A central pnicogen Z atom (Z = Sb, As) is covalently attached to the O atom of three -O(CH₂)_nX chains where X represents either an aldehyde or amine group. The chain can fold around so that the basic X group can engage in a noncovalent pnicogen bond with the central Z. The formation of up to three pnicogen bonds is energetically favored. The amine appears to engage in stronger pnicogen bonds than does the aldehyde, and bonds to Sb are favored over As, but there is little dependence on the length of the chain. The formation of each successive pnicogen bond reduces the magnitude of the σ-holes surrounding the Z atom, which tends to weaken the attraction for the basic end of the chain.

How Many Pnicogen Bonds can be Formed to a Central Atom Simultaneously?

A central ZF₃ molecule (Z = P, As, Sb, Bi) is allowed to interact with a number of nucleophiles exemplified by NCH, NH₃, and NC^– anion. The Z···N pnicogen bond (ZB) to a single base grows stronger for heavier Z atom: P < A < Sb < Bi and follows the NCH < NH₃ < NC^– order for the three bases. The maximum number of ZBs depends on both the nature of the base and pnicogen atom. PF₃ and AsF₃ can pnicogen bond with only a single CN^–; SbF₃ and BiF₃ can interact with two anions but only weakly. The weak NCH nucleophile can engage in a maximum of two ZBs, while three ZBs occur for NH₃. The latter NH₃ maximum can be extended to four ZBs but only for BiF₃. The fourth ZB is somewhat longer and weaker than the others, and the entire (H₃N)₄···BiF₃ complex relies partially on secondary interactions for its stability.