Symmetry breaking in the planar configurations of disilicon tetrahalides: Pseudo-Jahn–Teller effect parameters, hardness and electronegativity†
Abstract
CCSD(T), MP2, LC-BLYP, LC-ωPBE and B3LYP methods with the Def2-TZVPP basis set and natural bond orbital (NBO) interpretations were performed to investigate the correlations between the Pseudo-Jahn–Teller Effect (PJTE) parameters [i.e. vibronic coupling constant values (F), energy gaps between reference states (Δ) and the primary force constant (K0)], structural and configurational properties, global hardness, global electronegativity, natural bond orders, stabilization energies associated with electron delocalizations and natural atomic charges of disilicon tetrafluoride (1), disilicon tetrachloride (2), disilicon tetrabromide (3) and disilicon tetraiodide (4). All levels of theory showed the trans-bent (C2h) configurations as the energy minimum structures of compounds 1–4, and the flap angles between the X2Si planes and the SiSi bonds in the distorted (C2h) configurations decrease from compound 1 to compound 4. The negative curvatures of the ground state electronic configurations and the positive curvatures of the excited states of the adiabatic potential energy surfaces (APESs) which resulted from the mixing of the ground Ag and excited B2g states are due to the PJTE (i.e. PJT(Ag + B2g) ⊗ b2g problem). Contrary to the usual expectation, with the decrease of the energy gaps between reference states (Δ), the PJTE stabilization energy, EPJT, decreases from compound 1 to compound 4. The canonical molecular orbital (CMO) analysis revealed that the contributions of the ψHOMO(b3u) and ψLUMO(b1u) molecular orbitals in the vibronic coupling constant (F) decrease from compound 1 to compound 4. This fact clearly justifies the decrease of the vibronic coupling constant (F) and the primary force constant (the force constant without the PJTE) values on going from compound 1 to compound 4, leading to the decrease of the negative curvatures of the ground state electronic configuration curves of their corresponding APESs. The results obtained showed that the stabilization energies associated with the mixing of the distorted donor πSi–Si(bu) bonding and acceptor σSi–Si*(bu) antibonding orbitals along the b2g bending distortions decrease from compound 1 to compound 4. This fact reasonably explains the increase of the Si–Si natural bond orders (nbo) on going from compound 1 to compound 4. With the increase of the Si–Si natural bond orders, the corresponding EPJT decreases from compound 1 to compound 4. Importantly, the variations of the global hardness (η) differences (Δ[η(C2h) − η(D2h)]) do not correlate with the trend observed for their corresponding total energy differences, justifying that the configurational properties of compounds 1–4 do not obey the maximum hardness principle. Interestingly, the trans-bent (C2h) configurations of compounds 1–4 are more electronegative than their corresponding planar (D2h) forms and the variations of their global electronegativity (χ) differences (Δ[χ(C2h) − χ(D2h)]) succeed in accounting for the decrease of the EPJT stabilization energies for the D2h → C2h conversion processes on going from compound 1 to compound 4.