Cui Li-Juan, Liu Yu-Qian, Pan Sudip, Cui Zhong-Hua, Frenking Gernot
Institute of Atomic and Molecular Physics, Jilin University Changchun 130023 China
Key Laboratory of Physics and Technology for Advanced Batteries (Ministry of Education), Jilin University Changchun 130023 China.
Chem Sci. 2024 Jul 31;15(36):14705-20. doi: 10.1039/d4sc01979b.
Quantum chemical calculations are reported for the complexes of alkaline earth metals AeOLi (Ae = Be-Ba) at the BP86-D3(BJ)/def2-QZVPP and CCSD(T)/def2-QZVPPQZVPP levels. The nature of the Ae-OLi bond has been analyzed with a variety of methods. The AeOLi molecules exhibit an unprecedented σ donor bond Ae→OLi where the ()s lone-pair electrons of the Ae atom are donated to vacant O-Li antibonding orbitals having the largest coefficient at lithium. This is a covalent bond where the accumulation of the associated electronic charge is located at two positions above and below the Ae-OLi axis. The bifurcated component of orbital interactions is structurally related to the recently proposed collective bonding model, but exhibits a completely different type of bonding. The most stable isomer of AeOLi has a geometry and a singlet (A) electronic ground state. The bond dissociation energy (BDE) of the Ae-OLi bonds exhibits a zig-zag trend from BeOLi to BaOLi, with BeOLi having the largest BDE ( = 73.0 kcal mol) and MgOLi possessing the lowest BDE ( = 42.3 kcal mol) at the CCSD(T) level. The calculation of the atomic partial charges by the Hirshfeld and Voronoi methods suggests that Be and Mg carry small negative charges in the lighter molecules whereas the heavier atoms Ca-Ba have small positive charges. In contrast, the NBO and QTAIM methods give positive charges for all Ae atoms that are larger for Ca-Ba than that calculated by the Hirshfeld and Voronoi approaches. The molecules AeOLi have large dipole moments where the negative end is at the Ae atom with the polarity Ae→OLi. The largest dipole moments are predicted for the lighter species BeOLi and MgOLi and the smallest value is calculated for BaOLi. The calculation of the vibrational spectra shows a significant red-shift toward lower wave numbers for the Ae-OLi stretching mode in comparison to diatomic AeO. Besides the Ae→OLi σ-donor bonds there are also three dative bonds due to Ae←OLi backdonation which consist of one σ bond and two π bonds. The appearance of strong Ae→OLi σ donation leads to quadruple bonds AeOLi in all systems AeOLi, even for the lightest species with Ae = Be, Mg. The valence orbitals of Ca, Sr, and Ba, which are involved in the dative interactions, are the ()s and (-1)d AOs whereas Be and Mg use their ()s and ()p AOs. The EDA-NOCV results are supported by the AdNDP calculations which give four 2c-2e bonding orbitals. Three bonding orbitals have occupation numbers ∼2. One σ orbital has smaller occupation numbers between 1.32 and 1.73 due to the delocalization to the lithium atoms. The analysis of the electronic structure with the ELF method suggests multicenter bonds with mainly trisynaptic and tetrasynaptic basins, which also support the results of the EDA-NOCV calculations.
报道了在BP86-D3(BJ)/def2-QZVPP和CCSD(T)/def2-QZVPPQZVPP水平下碱土金属AeOLi(Ae = Be - Ba)配合物的量子化学计算。已用多种方法分析了Ae - OLi键的性质。AeOLi分子呈现出前所未有的σ给体键Ae→OLi,其中Ae原子的()s孤对电子被捐赠给在锂原子处具有最大系数的空O - Li反键轨道。这是一种共价键,相关电子电荷的积累位于Ae - OLi轴上方和下方的两个位置。轨道相互作用的分叉成分在结构上与最近提出的集体键合模型相关,但呈现出完全不同类型的键合。AeOLi最稳定的异构体具有 几何结构和单重态(A)电子基态。在CCSD(T)水平下,Ae - OLi键的键解离能(BDE)从BeOLi到BaOLi呈现出锯齿形趋势,BeOLi具有最大的BDE( = 73.0 kcal mol),而MgOLi具有最低的BDE( = 42.3 kcal mol)。通过Hirshfeld和Voronoi方法计算的原子部分电荷表明,在较轻的分子中Be和Mg带有小的负电荷,而较重的原子Ca - Ba带有小的正电荷。相比之下,NBO和QTAIM方法给出所有Ae原子的正电荷,对于Ca - Ba,其正电荷比通过Hirshfeld和Voronoi方法计算的更大。分子AeOLi具有大的偶极矩,负极在Ae原子处,极性为Ae→OLi。预测较轻的物种BeOLi和MgOLi具有最大的偶极矩,而计算得出BaOLi的偶极矩最小。振动光谱的计算表明,与双原子AeO相比,Ae - OLi伸缩模式向更低波数有明显的红移。除了Ae→OLi σ给体键外,由于Ae←OLi反馈捐赠还存在三个配位键,其中包括一个σ键和两个π键。强Ae→OLi σ捐赠的出现导致在所有AeOLi体系中形成四重键AeOLi,即使对于最轻的物种Ae = Be、Mg也是如此。参与配位相互作用的Ca、Sr和Ba的价轨道是()s和(-1)d原子轨道,而Be和Mg使用它们的()s和()p原子轨道。EDA - NOCV结果得到AdNDP计算的支持,AdNDP计算给出四个2c - 2e键合轨道。三个键合轨道的占据数约为2。由于向锂原子的离域,一个σ轨道的占据数较小,在1.32和1.73之间。用ELF方法对电子结构的分析表明存在主要具有三中心和四中心盆地的多中心键,这也支持了EDA - NOCV计算的结果。