The bifurcated σ-hole···σ-hole stacking interactions between organosulfur molecules, which are key components of organic optical and electronic materials, were investigated by using a combined method of the Cambridge Structural Database search and quantum chemical calculation. Due to the geometric constraints, the binding energy of one bifurcated σ-hole···σ-hole stacking interaction is in general smaller than the sum of the binding energies of two free monofurcated σ-hole···σ-hole stacking interactions. The bifurcated σ-hole···σ-hole stacking interactions are still of the dispersion-dominated noncovalent interactions. However, in contrast to the linear monofurcated σ-hole···σ-hole stacking interaction, the contribution of the electrostatic energy to the total attractive interaction energy increases significantly and the dispersion component of the total attractive interaction energy decreases significantly for the bifurcated σ-hole···σ-hole stacking interaction. Another important finding of this study is that the low-cost spin-component scaled zeroth-order symmetry-adapted perturbation theory performs perfectly in the study of the bifurcated σ-hole···σ-hole stacking interactions. This work will provide valuable information for the design and synthesis of novel organic optical and electronic materials.
Keywords: Cambridge structural database; bifurcated σ-hole···σ-hole stacking interaction; organosulfur molecules; quantum chemical calculation.