Chalcogen Bonds in Small-Organic Molecule Compounds Derived from the Cambridge Structural Database (CSD) ()
1. Introduction
There are numerous computational studies reporting on interactions of divalent chalcogen atoms (acting as σ-hole bond donors), mostly in small-inorganic molecule compounds (Pal, D. & Chakrabarti, P., 2001 [1]; Politzer, P. et al., 2007 [2]; Bleiholder, C. et al., 2006 [3]; Clark, T. et al., 2007 [4]; Ramos, L. A. et al., 2010 [5]; Politzer, P. et al., 2010 [6]; Politzer, P. et al., 2013 [7]; Politzer, P. & Murray, J. S., 2013 [8]; Politzer, P. et al., 2014 [9]; Pandiyan, B. V. et al., 2016 [10]; Zabardasti, A. et al., 2017 [11]; Bauza, A. & Frontera A., 2018 [12]; 2019 [13]; 2020 [14] ). Computational studies also showed that positive electrostatic potentials as well as σ-hole interactions are strongly affected by the environment of the covalently bonded chalcogen atom (Murray et al., 2007 [15]; Murray & Politzer, 2011 [16]; Politzer & Murray, 2017 [17] ).
To understand the intermolecular interaction between chalcogen atoms (O, S, Se, Te) quantum chemistry calculations using some model systems have been performed. Previous investigation on the intermolecular interactions between two molecules containing elements of group VI showed that the chalcogen-chalcogen interaction dominates and hydrogen bonding plays only a minor role in aggregates with heavier atoms (Bleiholder, C. et al., 2006 [3] ). The strength of the chalcogen interaction increased steadily from oxygen via sulphur to selenium and peaked around 6 kcal/mol for tellurium. The bond strength also increased when adding an electron withdrawing substituent such as the cyanide group CN.
However, experimental data supporting chalcogen-bonded complexes are rare and scarce (Goettel, J. T. & Gerken, M., 2016 [18]; Karjalainen, M. M. et al., 2016 [19] ). Analysis of consistent experimental data dealing with chalcogen interactions in organic molecules may be crucial due to their nature, importance and implications in the chemical, biological and medicinal areas. Theoretical study in addition to the data analysis is useful to provide a better understanding of the interaction and to guide its application in a wide variety of fields, like crystal engineering and drug design. The present study aims at using crystallographic data derived from the Cambridge Structural Database (CSD) to investigate on σ-hole interactions in small-organic molecule compounds containing divalent chalcogen atoms in conjunction with ab initio calculations in order to provide consistent insights into the nature of chalcogen bonding at divalent S, Se and Te atoms acting as σ-hole bond donors. This is of some interest as previous studies of inorganic chalcogen-bonded complexes reported that the chalcogen bond was a relatively strong and nearly linear interaction. They also showed the dependence of the interaction strength on the environment of chalcogen atoms [15] [16] [17].
2. Materials and Methods
2.1. CSD Analysis
Work CSD version 5.41 plus November 2019 update was used to perform the analysis of sulfur, selenium and tellurium compounds, this version contains a total of 1,034,174 structural entries. The CSD-System program ConQuest (Bruno et al., 2002 [20] ) was used to question the database substructure searches to match selenium containing systems, as well as relevant non-bonded interactions. The following secondary search criteria were employed: 1) atomic coordinates present in the entry; 2) error-free after CSD evaluation; 3) no disorder in the crystal structure; 4) no polymeric (catena) bonding; 5) no powder studies; 6) crystallographic R-factor < 0.10; 7) only organic structures (according to standard CSD definitions). Any specific variations of these criteria are noted in the text. In the results and discussion section, atoms and groups will be referred to by their coordination number as a superscript, with N3, for example, denoting a three-coordinate nitrogen atom. The CSD-System program Mercury (Macraeet al., 2006 [21], 2008 [22] ) was used for structure visualization, while numerical analysis of geometrical data retrieved by ConQuest was carried out using the Data Analysis Module implemented within Mercury (Sykes et al., 2011 [23] ).
The geometrical parameters δ (intermolecular distance) and θ (angle) of the chalcogen interaction are defined according to Scheme 1. The differences between sums of van der Waals radii of interacting atoms (rCh + rQA) and distances δ represent the “normalized” distances Δ (Wood et al., 2008 [24] ). The values of van der Waals radii reported by Bondi et al. (1964) [25] were used for calculations.
![]()
Scheme 1. Ch = chalcogen atom S, Se or Te; R1, R2 = substituents as illustrated in Table 1; QA = nucleophile; θ = interaction angle; δ = interaction distance.
2.2. Computational Study
The second-order Moller-Plesset perturbation theory MP2 and the density functional theories B3LYP and B3PW91 at basis set 6-311G (3df,2p) were used for calculations, when investigating on S- and Se-containing chalcogen bond donors. Only HF/ chemical model was used for Te chalcogen donors containing Te atoms.
We use Equations (1) and (2) (Murray et al., 2011 [26]; Lundemba S. et al., 2020 [27] ) to perform electrostatic potential (V) in the space around a chalcogen-bond donor, written in atomic units (au), and interaction energy (Eint) between a chalcogen donor and an acceptor, respectively. Energies Eint are energy minima at 0 K. Large basis sets are used to minimize the problem of correction for the basis set superposition error (Grimme S., 2006 [28]; Politzer P. et al., 2013 [7] ).
(1)
(2)
Figure 1 reports examples of molecular models of chalcogen bond used for ab initio calculations.
![]()
Figure 1. Molecular models of chalcogen-bond donors derived from CSD used for ab initio calculation.
3. Results and Discussion
3.1. CSD Surveys
3.1.1. Frequency of Occurrence (FoO)
CSD analysis reveals that many organic compounds containing divalent chalcogen atoms R1-X-R2 (X = S, Se or Te) interact with various nucleophile atoms (QA) like N, O, F, P, S, Cl, As, Se, Br, Sb, Te or I. Occurrence frequencies of molecular fragments containing divalent chalcogens R1-X-R2 (X = S, Se or Te) involved in interactions with σ-hole acceptors are listed in Table 1(a)-(c).
(a) R1-S-R2
(b) R1-Se-R2
(c) R1-Te-R2![]()
Table 1. Number (Nf) and occurrence frequency (FoO in %) of fragments forming intermolecular chalcogen-bonds.
According to those results shows, C, C-X compounds represent the most important group forming intermolecular interactions. It is noteworthy that the hybridization of carbon atoms affects the frequency of occurrence. Analysis of the carbon atom hybridization reveals that the occurrence frequency increases in the sequence X-Csp < X-Csp3 < X-Csp2.
3.1.2. Interaction Geometry
Figure 2 shows the propancy of the nucleophile bonding at chalcogen atoms, mainly along the extensions of the covalent bonds X-R.
![]()
Figure 2. Intermolecular interactions of the C, C-S chalcogen bond donor with nucleophiles.
It can be seen that nucleophiles mostly interact with chalcogen atoms of the σ-hole donors R1-X-R2 along the extension of the X-R covalent bond.
The geometrical parameters θ and δ as described in scheme 1 as well as the length of the chalcogen interaction Δ are summarized in Table 2.
![]()
Table 2. Values of the geometrical parameters θ, δ and Δ of the chalcogen interaction.
The results in Table 2 report mean bond angles θ > 160˚, which support nearly linear interactions (Politzer et al. 2013 [7] ) as illustrated by the pattern of O bonding at S atom in Figure 3.
![]()
Figure 3. Bond pattern of S-bonded complex.
Nucleophile atom bonding at chalcogen atom with such θ values is one of main σ-hole bond characteristics (Politzer et al. 2009 [29], 2013 [7] ).
Geometric data also show that mean interaction distances δ and lengths Δ have values less than the sums of van der Waals radii of the interacting atoms, demonstrating that the chalcogen interactions are relatively strong. Comparison of bond Δ values reveals that the ability of chalcogen atoms to interact with nucleophiles increases in the sequence S < Se < Te for a similar environment.
3.2. Ab Initio Calculations
3.2.1. Molecular Electrostatic Potential (MEP) and Lowest Molecular Orbital (LUMO) of Chalcogen Bond Donor
The lowest molecular orbital (LUMO) of chalcogen bond donors shows a typical pattern depicted in Figure 5 for NC-X-CN (X = S, Se and Te), while Figure 4 illustrates calculated molecular electrostatic potential surfaces.
As it can be seen, the patterns of MEP and LUMO predict a directional bonding of nucleophiles at divalent chalcogen atoms. Indeed, MEP in Figure 4 features two positive regions, which are located on the chalcogen atoms, and consequently, allow electrostatic interactions with nucleophiles.
![]()
![]()
![]()
Figure 4. Molecular electrostatic potential surfaces calculated for typical chalcogen S, Se and Te bond donor.
It is noteworthy that positive region are more or less broad depending on the environment of the chalcogen atom. Values of the most positive electrostatic potential (VS,max) calculated for some examples of chalcogen bond donors are listed in Table 3.
VS,max values show that the magnitude of the most positive electrostatic potential depends on chalcogen bond donor and substituents R1 and R2, enhanced by electron withdrawing substituents as previously reported by Murray et al. (2009 [30], 2014 [31] ) and Politzer et al. (2007 [2], 2013 [7], 2017 [17] ). It increases in the sequence S < Se < Te for similar chalcogen environments. One may then expect an electrostatic interaction between chalcogen-bond donors and nucleophiles (Politzer et al. 2017 [17] ).
The lowest molecular orbital (LUMO) of chalcogen bond donors shows a typical pattern depicted in Figure 5 for NC-X-CN (X = S, Se and Te).
![]()
Figure 5. LUMO’s examples for chalcogen-bond donors: CN-S-CN, CN-Se-CN and CN-Te-CN.
Interestingly, the lowest unoccupied molecular orbital also lies along extension of the covalent bond X-R like MEP does. Thus, both MEP and LUMO of chalcogen-bond donor predict directional bonding of nucleophiles at divalent chalcogen atoms as previously reported (Pal et al. 2001 [1]; Lundemba et al. 2020 [27] ).
3.2.2. Interaction Energetics
Values of interaction energy (Eint) calculated using different theories are listed in Table 4.
Energy values are negative demonstrating attractive interactions between chalcogen donors and acceptors. Energies calculated with MP2 have values close to those calculated using B3LYP with dispersion correction. but markedly different to those without dispersion correction. That finding reveals that dispersion forces play an important role in stabilizing chalcogen-bonded complexes. Eint depends on chalcogen atoms and their environments similarly to bond length (Δ) and Vs.max. It increases according to Vs.max in the sequence S < Se for interactions involving σ–hole donors containing S and Se divalent atoms.
![]()
Table 4. Values of Eint (kcal·mol−1), θ (˚), δ (Å) and number of imaginary frequencies (Nf).
Table 5 reports results of NBO analysis for some chalcogen-bonded complexes.
![]()
Table 5. Charge transfer and stabilizing energy E(2) of chalcogen interactions.
According to the results reported in Table 4. The chalcogen-bond consists of charge transfer (CT) from a nucleophile lone pair (LP) acting as σ-hole bond acceptor to an antibonding orbital (BD*) X-R of the chalcogen bond donor (Guo et al. 2015 [32] ). Furthermore, orbital contributions E(2) associated with electron delocalisation between depend on chalcogen atoms increasing as expected in the sequence S < Se.
4. Conclusions
We investigated on the chalcogen bonding at divalent S, Se and Te atoms in organic molecules, based on CSD analysis in conjunction with ab initio calculations using some chemical models. The work aimed to enable a better understanding of the interaction nature, geometry and strength. Several conclusions can be drawn from the results obtained:
1) CSD surveys showed that nucleophile atoms N. O. F. P. Cl. As. Se. Te and I interacted with the divalent chalcogen atoms. A significant frequency of occurrence (FoO) was found for 7 subgroups with (C; C), (C. N) and (N. N) as the most common σ-hole bond donors and the FoO decreasing in the sequence (C; C) > (N. N) > (C. N).
2) FoO and geometry were dependent on the R1 and R2 substituents. Values of the Δ interaction distance ranged from −0.052 Å to −0.512 Å for sulphur, from −0.086 Å to −1.017 Å for selenium and from −0.160 Å to −0.951 Å for tellurium, while those of the angle θ were close to 170˚ suggesting an attractive and nearly linear chalcogen interaction.
3) Ab initio calculations showed that the most positive electrostatic potential surfaces (VS.max) were located on chalcogens (S. Se; Te) in the extensions X-R1 and X-R2 bonds predicting directional and electrostatic interactions between the chalcogen atoms X and lone pairs of nucleophile atoms. LUMO orbitals of donors, with energy values decreasing in the sequence Te > Se > S, were involved in interactions nucleophiles.
4) All values of the interaction energy Eint were negative confirming the attractive nature of interaction. Comparison of Eint values calculated with the different quantum mechanical methods suggested that dispersion interactions played an important role in stabilizing σ-hole-bonded complexes. NBO analysis revealed that the chalcogen interactions consisted of transfer of charge from a nucleophile lone pair to an X-R1 and X-R2 antibonding orbital.
The present study provides a reliable contribution to the structural chemistry of chalcogen-bonded complexes. The relatively marked strength of the chalcogen bond and its characteristic directionality should guide the use of divalent S. Se and Te atoms in varied application areas like crystal engineering, drug design and supramolecular chemistry.
Acknowledgements
Albert L. Lundemba gratefully acknowledges support from the Cambridge Crystallographic Data Centre.