Correlations between hardness, electronegativity, anomeric effect associated with electron delocalizations and electrostatic interactions in 1,4,5,8-tetraoxadecalin and its analogs containing S and Se atoms

Hybrid density functional theory (B3LYP/6-311++G∗∗) based method and natural bond orbital (NBO) interpretation were used to investigate the correlations between the global hardness (η), global electronegativity (χ), anomeric effect associated with electron delocalization and electrostatic model associated with dipole–dipole interactions (four important conceptual descriptor in chemistry and physics) in 1,4,5,8-tetraoxadecalin (1), 1,4,5,8-tetrathiadecalin (2) and 1,4,5,8-tetraselenadecalin (3). B3LYP/6-311++G∗∗ results showed that the calculated Gibbs free energy and corrected electronic energy differences between the cis- and trans-stereoisomers [i.e. ΔG = Gtrans − Gcis, ΔEo = Eo(trans) − Eo(cis)] decrease from compound 1 to compound 3. On the other hand, the global hardness (η) differences between the cis- and trans-stereoisomers (i.e. Δ[η(cis) − η(trans)]) decrease from compound 1 to compound 2 but increase from compound 2 to compound 3. Accordingly, the variation of ΔG = Gtrans − Gcis, ΔEo = Eo(trans) − Eo(cis) parameters from compound 1 to compound 3 cannot be justified by the variation of their corresponding global hardness. Therefore, compounds 2 and 3 do not obey the maximum hardness principle. Interestingly, the calculated global electronegativity (χ) differences between the cis- and trans-stereoisomers (i.e. Δ[χ(cis) − χ(trans)]) decrease from compound 1 to compound 3. NBO results showed that the anomeric effect associated with the electron delocalization is in favor of the cis stereoisomers. Effectively, with the decrease of the anomeric effect, ΔG = Gtrans − Gcis, ΔEo = Eo(trans) − Eo(cis) parameters decrease from compound 1 to compound 3. Therefore, the variations of the anomeric effect and global electronegativity explain the variations of ΔG = Gtrans − Gcis, ΔEo = Eo(trans) − Eo(cis) parameters from compound 1 to compound 3. It should be noted that the variation of the calculated dipole moment differences between the cis and trans stereoisomers [i.e. Δ(μ(cis) − μ(trans)] does not justify the variation of ΔG = Gtrans − Gcis, ΔEo = Eo(trans) − Eo(cis) parameters from compound 1 to compound 3. The total steric exchange energy differences between the cis- and trans-stereoisomers (Δ[TSEEcis − TSEEtrans]) increase from compound 1 to compound 3 which in the trend observed for the anomeric effect and global electronegativity. The correlations between the global hardness, global electronegativity, anomeric effect, zero point energies, electrostatic model, structural parameters, corrected electronic energies (ΔEo), and thermodynamic parameters [ΔH, ΔG and ΔS] of compounds 1–3 have been investigated.