VB2n-(n=8~12)团簇几何结构、电子及热力学特性的理论研究

引用本文: 李成刚, 崔颍琦, 田浩, 张洁, 邵琴琴, 申梓刚, 任保增. VB_2n^-(n=8~12)团簇几何结构、电子及热力学特性的理论研究[J]. 无机化学学报, 2022, 38(8): 1523-1532. doi: 10.11862/CJIC.2022.150 shu

Citation: Cheng-Gang LI, Ying-Qi CUI, Hao TIAN, Jie ZHANG, Qin-Qin SHAO, Zi-Gang SHEN, Bao-Zeng REN. Theoretical Study on Geometric Structures, Electronic, and Thermodynamics Properties of VB_2n^- (n=8-12) Clusters[J]. Chinese Journal of Inorganic Chemistry, 2022, 38(8): 1523-1532. doi: 10.11862/CJIC.2022.150 shu

VB_2n^-(n=8~12)团簇几何结构、电子及热力学特性的理论研究

李成刚 ^{1, 2,} ,
崔颍琦 ^{1,
,} ,
田浩 ^1, ,
张洁 ^1, ,
邵琴琴 ^1, ,
申梓刚 ^1, ,
任保增 ^2,

1.
郑州师范学院物理与电子工程学院，郑州 450044
2.
郑州大学化工与能源学院，郑州 450001

通讯作者: 崔颍琦, E-mail：yqcui2007@163.com

基金项目:
国家自然科学基金 11904328

国家自然科学基金 12104416

河南省高等学校青年骨干教师培养计划 2018GGJS163

郑州师范学院青年骨干教师培养计划 QNGG-211361

郑州师范学院本科教学改革研究项目 JXGG-20773

郑州师范学院优秀基层教学组织建设项目(物理与电子工程学院物理学教研室)

摘要: 基于卡里普索结构预测程序和密度泛函理论的第一性原理计算，搜索确定了VB_2n^-(n=8~12)团簇的基态和亚稳态结构。结果发现，V原子的掺杂完全改变了原硼团簇的结构并提高了原体系的稳定性。掺杂体系基态结构分别呈现高对称性的鼓状(VB₁₆^-：C_2v)、管状(VB₁₈^-：C_2v和VB₂₀^-：C_s)及笼状(VB₂₂^-：C₂和VB₂₄^-：D_3h)结构。基于基态结构，研究了体系的电荷转移和极化率，拟合出了光电子能谱、红外和拉曼谱图，分析了流变键和芳香特性。最后，研究了体系的热力学特性，讨论了温度对热力学参数的影响。

关键词:

English

Theoretical Study on Geometric Structures, Electronic, and Thermodynamics Properties of VB_2n^- (n=8-12) Clusters

1.
College of Physics and Electronic Engineering, Zhengzhou Normal University, Zhengzhou 450044, China
2.
School of Chemical Engineering and Energy, Zhengzhou University, Zhengzhou 450001, China

Corresponding author: Ying-Qi CUI, yqcui2007@163.com

Received Date: 08 February 2022
Revised Date: 09 May 2022
Available Online: 10 August 2022

Abstract: The geometric structures, electronic, and thermodynamic properties of vanadium doped boron clusters, VB_2n^- (n=8 -12), were investigated systemically by using CALYPSO searching method and density functional theory. It is found that vanadium atom doping significantly modifies the structures of the boron clusters and strong the chem- ical activity of systems. A drum-shaped structure is the global minimum for VB₁₆^- cluster with C_2v point symmetry. Tubular-shaped VB₁₈^- and VB₂₀^- with C_2v and C_s symmetry exhibit a metal-centred tubular with a B ₂ unit over the B₁₆ and B ₁₈ drum, respectively. For VB₂₂^- and VB₂₄^- clusters, vanadium atom tends to encapsulate into boron cages. Based on the lowest energy structures, the charge transfer and polarizability were explored, the photoelectron spec- tra, infrared spectra, and Raman spectra were simulated, the fluxional bonds and aromatic properties were analyzed. At last, the thermodynamic properties were investigated, the thermodynamics parameters were discussed for the lowest energy structures of VB_2n^- (n=8-12) clusters.

Key words:

0. 引言

随着实验精度的不断提升和计算机技术的迅速发展，具有新颖几何结构、独特物理和化学特性的团簇引起了研究者广泛关注^[1-3]。B原子的特殊电子结构，使其易与其它原子成键，在构建多功能纳米材料及调控其性能时表现出独特的优势，成为研究关注的重点之一。近年来，基于实验光电子能谱(photoelectron spectroscopy，PES)和理论第一性原理计算，大量平面或准平面、管状硼环及双链硼纳米带结构的硼团簇被陆续发现^[4-11]。随后，利用金属掺杂纯硼团簇产生了一系列新颖结构和特殊性质的B掺杂体系，并以此结构为基元，构建了众多功能和纳米新型材料^[12-21]。例如：基于实验的PES和密度泛函理论的第一性原理计算，美国布朗大学的王来生教授课题组发现了一系列特殊结构的B掺杂体系，例如：RhB₉^-、IrB₉^-、RhB₉、IrB₉、CoB₈^-、RuB₉^-、TaB₁₀^-和NbB₁₀^-等团簇呈现掺杂原子位居中心的轮状平面结构^[12-14]，CoB₁₂^-和RhB₁₂^-的基态呈现半三明治夹心结构^[15]，AlB₇^- (C_7v)和AlB₈^- (C_8v)团簇的基态具有高对称性的伞状结构(Al原子位居伞的顶端)^[16]。最近，又发现了高配位数CoB₁₆^-具有高对称性(能量几乎简并D_8d和C_4v)鼓状结构^[17]。此外，其他课题组陆续发现了类平面(RhB₁₈^-)、双环状(LiB₂₀、NaB₂₀、KB₂₀、CaB₂₀、MgB₂₀)、笼状(TiB₂₂和ScB₂₂)及三环管状(Li₂B₂₄)等结构^[18-21]。

具有空d电子轨道的V元素在地壳中具有很高的丰度，易于同其它材料结合形成新颖结构材料。近年来，众多科研工作者对过渡金属V原子掺杂硼团簇进行了深入研究^[22-24]。Li课题组研究了VB₁₀^-团簇的结构及芳香性，结果发现VB₁₀^-基态具有船型结构^[22]。Chen等研究发现VB_n⁺ (n=3~6)同甲烷反应的过程中，很容易使甲烷脱氢，该研究为甲烷的催化和转换指明了方向^[23]。2019年，Pham课题组基于第一性原理计算，研究得知V₂B₇⁰/V₂B₇⁺/V₂B₇^-具有陀螺型结构^[24]。基于此结构，对体系的磁性、芳香性及与二氧化碳分子的作用进行了深入探究。此外，基于量子化学计算，Tran等研究V₂B₇⁰/V₂B₇⁺/V₂B₇^-团簇的几何和电子结构^[25]，优化结构发现，中性和阳性体系的基态具有C_s点群对称的平面结构。此外，分析了V₂B₇⁺团簇同甲烷反应的机理。最近，我们课题组对过渡金属V掺杂硼团簇的结构和性质进行了一系列研究^[26-28]。例如，结合密度泛函理论和卡里普索(crystal structure analysis by particle swarm optimization，CALYPSO)结构预测程序，优化得到了拥有高对称性(C_2v点群对称)管状结构的VB₁₆^-^[26]；紧接着，在VB₁₆^-高对称结构基础上，通过堆积的方式获得了高对称V₂B₂₄^-和V₃B₃₂^- (C_4h和D_8d点群对称)团簇^[27]。最近，基于相同的方法，系统搜索了VB₂₀^-团簇的几何结构，优化得到了高对称性管状分子马达^[28]。在此结构基础上，研究了电子、光谱特性，揭示了其结构流变的本质。目前为止，对不同尺寸下V掺杂硼团簇几何结构、电子和热力学特性的研究未见报道。本文中，基于前期工作基础，利用卡里普索结构预测程序，并结合密度泛函理论的第一性原理计算，首先搜索了不同尺寸下VB_2n^-(n=8、9、10、11、12)团簇的基态和亚稳态结构。基于基态结构，分析了不同尺寸下的团簇的相对稳定性，研究了掺杂体系的电荷转移和极化率，拟合出了PES、红外和拉曼谱图，分析了流变键和芳香特性。最后讨论了体系的热力学特性和温度对热力学参数的影响。

1. 计算方法

研究可知，合理基态结构的确定是研究其物理和化学特性的前提和基础。采用卡里普索结构预测程序对VB_2n^- (n=8~12)团簇的结构进行了系统搜索^[29-31]。具体搜索过程如下：首先，基于不同自旋多重态(2S+1=1、3、5、7)，在PBE/3-21G水平下对初始结构进行搜索，根据设置的种群数(20个)和演化的代数(50代)，可以得到相应数量的初始结构(总结构数=种群数×代数)，从上述初始结构中筛选出能量差超过0.3 eV的50个结构。随后，在高精度的PBE0/6- 311+G(d) 水平下对筛选的50个结构进行重新优化^[32-34]，根据优化后能量大小关系，即可得到研究对象的基态和亚稳态结构。对于此预测程序的具体描述可以参考相关已发表的文献^{[26-28, 35-47]}。此外，为了考虑不同泛函对优化结果的影响，在相同的基组下，利用TPSSH泛函对50个能量差超过0.3 eV的结构进行了优化^[48]。众多研究已证明TPSSH泛函可以应用到过渡金属原子掺杂硼团簇的研究^{[24, 49-51]}。所有优化过程中不设置对称性限制并考虑了自旋多重度的影响，上述运算均在Gaussian 09程序下运行^[52]。

2. 结果与讨论

2.1 VB_2n^-(n=8~12)团簇结构和稳定性

通过卡里普索结构预测程序，在PBE0/6-311+G(d)和TPSSH/6-311+G(d)水平下对纯硼B_2n^-和VB_2n^-(n=8~12)团簇的结构进行了搜索优化，优化后基态及亚稳态结构见图 1。首先，优化得到的纯硼团簇的基态和美国布朗大学王来生教授课题组通过PES实验与理论研究所得的结果完全吻合^[52]，这充分证明所采用泛函基组的合理性。对于掺杂体系，本课题组已对VB₁₆^-和VB₂₀^-的基态结构做过详细报道^{[26, 28]}。其次，对于掺杂体系，2种泛函下的基态和亚稳态完全相同，说明优化的基态结构是合理的。最后，对比掺杂前后结构发现，V掺杂后完全改变了原团簇的结构，形成了具有高对称性的鼓状(VB₁₆^-)、管状(VB₁₈^-和VB₂₀^-)和笼式(VB₂₂^-和VB₂₄^-)结构，而且，掺杂的V原子都位居内部，每个团簇基态对应的点群对称性分别为C_2v、C_2v、C_s、C₂和D_3h。此外，图 2给出了VB_2n^- (n=8~12)团簇的分子轨道能级，对应的HOMO-LUMO能隙分别为2.076、3.039、2.987、3.283和4.059 eV；而对应未掺杂体系B_2n^-(n=8~12)团簇的HOMO-LUMO能隙分别为1.671、1.493、1.578、1.542和1.544 eV，数据对比发现，掺杂后体系的HOMO- LUMO能隙明显增大，说明V掺杂硼团簇VB_2n^-(n=8~12)的稳定性提高了。

图 1

图 1. VB_2n^- (n=8~12)团簇的基态及亚稳态结构、点群对称、电子组态及PBE0/6-311+G(d)和TPSSH/6-311+G(d)水平下的相对能量

Figure 1. Lowest energy and low-lying structures VB_2n^- (n=8-12) clusters along with the point group symmetry, electronic state, and relative energies at PBE0/6-311+G(d) and TPSSH/6-311+G(d) levels

V and B atoms are marked in blue and green, respectively; All units in eV

下载: 全尺寸图片幻灯片

图 2

图 2. VB_2n^- (n=8~12)团簇基态结构的分子轨道能级

Figure 2. Molecular orbital energy levels of the lowest energy structures of VB_2n^- (n=8-12) clusters

下载: 全尺寸图片幻灯片

为了研究体系的稳定性，在PBE0和TPSSH泛函下计算了体系的平均束缚能(E_b)：

$ {{E}_{\text{b}, n}}=\frac{{{E}_{\text{V}}}+n{{E}_{{{\text{B}}^{-}}}}-{{E}_{\text{V}{{\text{B}}_{n}}}}^{-}}{n+1} $

$ {{E}_{\text{b}, n+1}}=\frac{n{{E}_{{{\text{B}}^{-}}}}-{{E}_{{{\text{B}}_{n}}-}}}{n+1} $

其中，E为对应原子和团簇的能量。图 3给出了平均束缚能随尺寸变化的关系，从中可以看出，掺杂体系的平均束缚能大于纯硼团簇的平均束缚能，说明掺杂体系的稳定性强于对应的纯硼团簇的稳定性，该结论同HOMO-LUMO能隙分析结论完全相同。同时，表 1给出了不同尺寸团簇中V—B和B— B键的Wiberg键级的大小，从中可以看出，每个掺杂体系中V—B键级大于B—B键级，说明V—B键比B—B键更稳定。

图 3

图 3. PBE0和TPSSH泛函下VB_2n^- (n=8~12)团簇基态结构平均束缚能随尺寸变化

Figure 3. Clusters size dependence on the binding energy of the lowest energy structures of VB_2n^- (n=8-12) clusters

下载: 全尺寸图片幻灯片

表 1

表 1 VB_2n^- (n=8~12)团簇基态结构中V原子的净电荷(Q)、NEC、键长和Wiberg能级

Table 1. Atomic net charges (Q) of V atoms, NEC, bond lengths, and Wiberg bond order of the lowest energy structures of VB_2n^- (n=8-12) clusters

下载: 导出CSV

VB_2n^-	Q / V		NEC(B): 2s22p1	NEC(V): 3d34s2	Average Wiberg bond order
VB_2n^-	NPA	Hirshfeld	NEC(B): 2s22p1	NEC(V): 3d34s2	V—B	B—B
VB₁₆^-	2.759	0.253	2_s0.67-0.702_p2.17-2.37	3_d3.944_s0.024_p0.01	0.521	0.245
VB₁₈^-	4.293	0.242	2_s0.61-0.702_p2.17-2.37	3_d4.044_s0.014_p0.01	0.488	0.273
VB₂₀^-	2.972	0.257	2_s0.67-0.732_p2.17-2.33	3_d3.944_s0.024_p0.01	0.424	0.199
VB₂₂^-	5.454	0.240	2_s0.53-0.792_p1.38-2.48	3_d4.264_s0.014_p0.01	0.391	0.183
VB₂₄^-	4.925	0.279	2_s0.61-0.632_p2.24-2.54	3_d4.064_s0.224_p0.01	0.345	0.167

2.2 电荷转移

为了了解团簇内部电荷转移信息，利用不同的电荷分析方法对电子在原子间及内层轨道上转移特征进行了分析，具体结果见表 1。对比发现，基于自然布居分布(natural population analysis，NPA) 和Hirshfeld分析方法得到的V原子的净电荷均为正值，说明在VB_2n^-(n=8~12)体系中，电子的转移方向为V→B，这与V的电负性(1.63 eV)小于B原子的电负性(2.04 eV)相吻合^[53]。众所周知，自由B原子的电子排布为1s²2s²2p¹，自由V原子的电子排布为[Ar]3d³4s²。基于自然电子组态(natural electron configu-ration，NEC)方法计算发现，B原子的B2s轨道失去电子，B2p轨道得到电子。然而，V3d轨道得到电子，V4s轨道失去电子。数据说明，电子从B2s和V4s轨道转移到B2p和V3d轨道，B原子和掺杂的V原子之间存在强烈的spd轨道杂化。此外，自由V原子的价电子排布为3d³4s²，当V原子掺杂硼团簇后，3d轨道得到大约1个电子，而4s轨道上的2个价电子几乎全部转移到其它轨道，即3d轨道得到一个电子，4s轨道失去2个电子，所以，VB_2n^-(n=8~12)体系中V均约为+1价。

2.3 极化率

为了考察VB_2n^-(n=8~12)对外场的响应，计算了各向同性平均极化率α、团簇平均极化率α、各向异性极化率Δα，其中α和Δα定义如下：

$ \begin{array}{c} & \bar{\alpha }=\frac{1}{3N}\left( {{\alpha }_{xx}}+{{\alpha }_{yy}}+{{\alpha }_{zz}} \right) \\ & \Delta \alpha =\sqrt{\frac{1}{2}\left[ {{\left( {{\alpha }_{xx}}-{{\alpha }_{yy}} \right)}^{2}}+{{\left( {{\alpha }_{yy}}-{{\alpha }_{zz}} \right)}^{2}}+{{\left( {{\alpha }_{z}}-{{\alpha }_{xx}} \right)}^{2}}+6\left( \alpha _{xy}^{2}+\alpha _{xz}^{2}+\alpha _{yz}^{2} \right) \right]} \\ \end{array} $

式中，N为团簇的总原子数，α_xx、α_xy、α_xz、α_yy、α_yz和α_zz为张量对角元，具体计算结果见表 2。对比数据发现，VB₁₆^-团簇的平均极化率最大，说明VB₁₆^-团簇的电子结构相对稳定性差，容易被外加场极化，非线性光学效应强。VB₂₄^-团簇平均极化率最小，其相对稳定性强，最不容易被外加场极化，非线性光学效应弱。其次，VB₁₆^-团簇各向异性极化率最大，说明对外场的各向异性响应较强；VB₂₄^-团簇各向异性极化率最小，说明对外场的各向异性响应最弱。最后，对比HOMO-LUMO能隙、平均极化率和各向异性极化率发现，三者之间存在反比关系，较大的HOMO-LUMO能隙对应较小的平均极化率和各向异性极化率，原因可能是在微弱外电场下，HOMO- LUMO能隙越大，电子跃迁时的阻力越大，因此，能转移的电子数越少，形成的平均极化率和各向异性极化率越小。

表 2

表 2 VB_2n^- (n=8~12)团簇基态对应的HOMO-LUMO能隙和极化率

Table 2. HOMO-LUMO energy gaps and polarizabilities of lowest energy structures of VB_2n^- (n=8-12) clusters

下载: 导出CSV

Parameter	VB₁₆^-	VB₁₈^-	VB₂₀^-	VB₂₂^-	VB₂₄^-
ΔE^* / eV	2.08	3.04	2.99	3.28	4.06
α_xx / nm³	2.75×10^-1	2.77×10^-1	2.29×10^-1	2.99×10^-1	3.08×10^-1
α_xy / nm³	0	0	0	0	0
α_xz / nm³	0	0	0	0	0
α_yy / nm³	2.75×10^-1	2.77×10^-1	3.19×10^-1	3.04×10^-1	3.08×10^-1
α_yz / nm³	0	0	0	0	0
α_zz / nm³	1.89×10^-1	2.11×10^-1	3.15×10^-1	2.83×10^-1	3.17×10^-1
α / nm³	2.46×10^-1	2.55×10^-1	2.88×10^-1	2.95×10^-1	3.11×10^-1
α/ nm³	4.34×10^-2	4.03×10^-2	4.11×10^-2	3.85×10^-2	3.73×10^-2
Δα / nm³	8.64×10^-2	6.64×10^-2	8.79×10^-2	1.93×10^-2	9.11×10^-2
^*Energy gap of HOMO-LUMO.

2.4 PES、红外、拉曼谱图

光电子技术是研究物质电子结构的重要方法，是鉴别其结构的一个重要手段。我们利用Multiwfn软件拟合出VB_2n^- (n=8~12)团簇的PES(图 4)^[54]，希望理论拟合的PES谱图为实验上鉴别VB_2n^-(n=8~12)团簇的基态结构提供有力的支持。

图 4

图 4. (a) VB₁₆^-、(b) VB₁₈^-、(c) VB₂₀^-、(d) VB₂₂^-及(e) VB₂₄^-团簇基态结构的PES谱图

Figure 4. PES spectra of the lowest energy structures for (a) VB₁₆^-, (b) VB₁₈^-, (c) VB₂₀^-, (d) VB₂₂^-, and (e) VB₂₄^- clusters

下载: 全尺寸图片幻灯片

根据已知的基态结构，计算并拟合出VB_2n^- (n= 8~12)团簇的红外和拉曼谱图(图 5和6)。首先，红外和拉曼谱图的主要特征峰集中在200~1 400 cm^-1范围内。对于鼓状结构VB₁₆^-，其红外谱图的主要特征峰位于461 cm^-1处，表现为所有B5、B7、B11和B17的左右摇摆振动，而剩余的B原子和掺杂V原子受迫做水平振动。拉曼谱图的主要特征峰在715 cm^-1处，表现为所有B原子围绕管中心V原子径向向里的呼吸振动模式，而中心的V原子几乎不动。VB₁₈^-的红外和拉曼谱图的主要特征峰在495和718 cm^-1处，分别为B₂单元中的B16和B17的左右摇摆振动和所有B原子径向向外的呼吸振动模式。对于VB₂₀^-，红外谱图的主要特征峰在386 cm^-1处，表现为中心V原子沿Z轴平面方向的左右水平振动，而剩余的B原子做受迫振动。拉曼谱图的主要特征峰在665 cm^-1处，所有B原子做径向向外的呼吸振动，其中，B11和B18沿着V—B键方向，上下环中的B原子沿着上下环的B—B键方向。VB₂₂^-的红外和拉曼谱图的主要特征峰分别在368和661 cm^-1处，分别对应所有原子的左右摇摆振动(掺杂V原子沿Z轴方向运动)和所有B原子沿不同径向方向的呼吸振动模式。对于VB₂₄^-，红外谱图的主要特征峰在332 cm^-1处，表现为中心V原子沿Z轴方向的左右水平振动模式，拉曼谱图的特征峰(624 cm^-1)表现为所有B原子沿B—V径向方向的呼吸振动模式。

图 5

图 5. (a) VB₁₆^-、(b) V B₁₈^-、(c) VB₂₀^-、(d) VB₂₂^-及(e) VB₂₄^-团簇基态结构的红外光谱和最强峰的振动模式图(插图)

Figure 5. IR spectra and vibrational mode of the strongest peak (Inset) for (a) VB₁₆^-, (b) VB₁₈^-, (c) VB₂₀^-, (d) VB₂₂^-, and (e) VB₂₄^- clusters

下载: 全尺寸图片幻灯片

图 6

图 6. (a) VB₁₆^-、(b) VB₁₈^-、(c) VB₂₀^-、(d) VB₂₂^-及(e) VB₂₄^-团簇基态结构的拉曼光谱和最强峰的振动模式图(插图)

Figure 6. Raman spectra and vibrational mode of the strongest peak (Inset) for (a) VB₁₆^-, (b) VB₁₈^-, (c) VB₂₀^-, (d) VB₂₂^-, and (e) VB₂₄^- clusters

下载: 全尺寸图片幻灯片

2.5 流变键

流变键(fluxional bonds)是由山西大学李思殿教授课题组首次提出^[55-57]，用以描述一定条件下通过协同机制不断断裂和生成的化学键。基于第一性原理计算和适应性自然密度划分(adaptive natural density partitioning，AdNDP)法，该课题组对准平面结构B₁₉^-、管状TaB₂₀^-、硼球烯B₃₉^-^[55]、半三明治结构MB₁₈^- (M=K、Rb、Cs)^[56]及NiB₁₁^-^[57]进行了详细的流变键分析。其中，研究管状TaB₂₀^-发现，上方的B₂单元与管底座[TaB₁₈]间形成的2个σ键在一定温度下均不断地断裂和生成，致使该团簇在基态和过渡态之间发生流变。鉴于VB₂₀^-和TaB₂₀^-团簇具有相同的价电子数(33对)，以VB₂₀^-团簇为例，通过对基态和过渡态进行结构的流变和AdNDP键合分析，研究了其流变的机理。对于基态结构，管上方B₂单元的B—B键长为0.159 7 nm，比与管连接的B—B键长(0.171 5和0.187 3 nm)短，说明B₂单元中B—B键能大于B₂单元和管连接的B—B键能，当发生流变时，B₂单元和管连接的B—B易断裂。首先，通过优化的基态结构确定了其过渡态结构，图 7为VB₂₀^-团簇在势能面上的结构流变图。从中可以看出，对于基态，管上方的B₂单元与管顶部B₉环之间有2个4c-2e σ键。当上方的B₂单元克服8.46 kJ·mol^-1的流变能垒后，可以绕底座[VB₁₈]开始顺时针旋转，进而流变到过渡态，此时2个4c-2e σ键就变成了2个3c-2e σ键。因此，VB₂₀^-团簇的基态与过渡态之间发生流变。

图 7

图 7. VB₂₀^-团簇基态的结构流变键图

Figure 7. Fluxional bonds of the lowest energy structure of VB₂₀^- cluster

3c-2e σ bonds and 4c-2e σ bonds formed atop the B₁₈ double-ring tube are highlighted in orange

下载: 全尺寸图片幻灯片

为了深入研究其流变机理，对VB₂₀^-的基态和过渡态进行了AdNDP分析(图 8)。其中，nc-2e代表n中心双电子键，红色和蓝色代表不同的相。为了方便对比VB₂₀^-的基态和过渡态中流变键的变化，图 8仅提供了前22对价电子数的AdNDP。对于基态和过渡态，首先，位于[VB₁₈]管上端的B—B形成了1个2c - 2e的双中心双电子σ键，占据数为1.82e和1.81e。其次，管上方B₂单元和中心V原子之间形成一个3c-2e σ键，对应的占据数为1.97e和1.96e。对于B₁₈的管状而言，每相邻的3个B原子之间构成了18个3c - 2e σ键，其占据数分别为1.91e~1.93e和1.90e~1.96e。对于过渡态，管顶部B₉中的B原子和B₂单元B原子之间构成了2个3c-2e σ键；对于基态，继续搜索可以得到管顶部B₉中的B原子和B₂单元B原子之间构成了2个4c-2e σ键，基态和过渡态中的占据数分别为1.79e和1.88e。总之，二者的不同之处在于，在基态中是2个4c-2e σ键，而在过渡态中是2个3c-2e σ键，该分析同上文中对流变键的分析相符。

图 8

图 8. VB₂₀^-团簇AdNDP成键图

Figure 8. AdNDP bonding patterns of VB₂₀^- cluster

ON is indicated occupation number and the green ball represents B atoms

下载: 全尺寸图片幻灯片

2.6 芳香性

核独立化学位移(nucleus independent chemical shifts，NICS)方法是一种分子芳香性的判断依据，广泛应用于有机无机化合物及团簇芳香性的研究中，NICS为负值表示芳香性，NICS为正值表示反芳香性，当NICS值接近零时为非芳香性。本工作中，计算NICS时，管状和笼状结构的参考点分别为管状和笼状结构的中心，试探原子Bq分别在其中心的0.05、0.10、0.15和0.20 nm处，具体计算结果见表 3。结果发现，VB₂₀^-和VB₂₄^-具有芳香性。

表 3

表 3 VB_2n^- (n=8~12)团簇基态结构的NICS值

Table 3. NICS values of the lowest energy structures of VB_2n^- (n=8-12) clusters

下载: 导出CSV

Cluster	NICS(0.5)_zz /ppm	NICS(1.0)_zz /ppm	NICS(1.5)_zz /ppm	NICS(2.0)_zz /ppm
VB₁₆^-	94.69	-127.38	-68.03	-31.80
VB₁₈^-	75.23	-6.68	-22.05	-9.65
VB₂₀^-	-157.12	-123.02	-86.06	-66.68
VB₂₂^-	25.18	-50.21	-27.63	-7.72
VB₂₄^-	-83.47	-66.61	-46.45	-15.92

2.7 热力学性质

在298.15 K、1.013×10⁵ Pa下，计算了VB_2n^-(n=8~ 12)团簇的热力学参数：定容热容C_v、标准熵S和标准生成焓ΔH，其中ΔH可表示如下：

$ \Delta H={{E}_{\text{M}{{\text{B}}_{2n}}}}-{{E}_{M}}-2n{{E}_{{{\text{B}}^{-}}}} $

计算结果见表 4，从中可以看出，所有体系的生成焓均为负值，说明生成团簇的反应都是放热反应，热力学上都是稳定的，而且，随着团簇尺寸的增大，稳定性逐渐提高。为了探究温度对C_v和S的影响，计算了不同温度下的C_v和S值，具体结果见表 5和图 9。结果发现，热力学参数C_v和温度之间存在二次函数关系，S和温度之间存在近似线性关系。

表 4

表 4 计算得到的VB_2n^- (n=8~12)团簇基态结构热力学参数

Table 4. Calculated thermodynamics parameters of the lowest energy structure of VB_2n^-(n=8-12) clusters

下载: 导出CSV

Parameter	ΔH /eV	C_v /(J·mol^-1·K^-1)	S /(J·mol^-1·K^-1)
VB₁₆^-	-84.14	212.82	420.86
VB₁₈^-	-94.98	224.85	406.81
VB₂₀^-	-104.98	270.09	465.57
VB₂₂^-	-115.47	280.06	454.79
VB₂₄^-	-125.48	299.41	460.61

表 5

表 5 不同温度下的定容热容和标准熵

Table 5. Constant heat capacities and standard entropies at different temperatures J·mol^-1·K^-1

下载: 导出CSV

T / ℃	VB₁₆^-		VB₁₈^-		VB₂₀^-		VB₂₂^-		VB₂₄^-
T / ℃	C_v	S	C_v	S	C_v	S	C_v	S	C_v	S
100	54.89	278.94	44.16	269.55	64.65	287.04	57.71	280.8	61.99	274.93
200	139.96	347.61	136.66	331.99	180.62	372.70	176.03	361.17	188.09	360.87
300	212.82	420.86	224.85	406.81	270.09	465.57	280.06	454.79	299.41	460.61
400	265.57	493.42	289.74	484.70	333.17	556.40	356.39	550.53	383.00	563.05
500	300.64	558.29	332.47	555.82	375.37	637.09	407.10	637.36	439.34	656.43
600	324.66	616.64	361.39	620.41	404.67	709.49	441.75	716.03	478.10	741.34
700	341.45	669.09	381.42	678.76	425.42	774.53	465.92	787.04	505.23	818.16
800	353.48	716.42	395.68	731.56	440.43	833.24	483.21	851.29	524.68	887.78

图 9

图 9. VB_2n^- (n=8~12)团簇的(a) C_v和(b) S随温度变化的曲线

Figure 9. Temperature dependence curves of (a) C_v and (b) S for VB_2n^- (n=8-12) clusters

下载: 全尺寸图片幻灯片

3. 结论

结合密度泛函理论和卡里普索结构预测程序，系统研究了过渡金属V原子掺杂硼的VB_2n^-(n=8~12)团簇的几何结构、稳定性、电子和热力学特性。结构优化发现，V的掺杂完全改变了纯硼团簇的结构，掺杂体系的基态结构分别呈现鼓状(VB₁₆^-：C_2v)、管状(VB₁₈^-：C_2v和VB₂₀^-：C_s)及笼状(VB₂₂^-：C₂和VB₂₄^-：D_3h)构型，掺杂提高了硼团簇的稳定性。电荷转移分析发现，在原子之间，电荷从V原子转移到硼原子；在原子内部轨道上，电子从B2s和V4s轨道向B2p和V3d轨道转移，B和V原子之间存在强烈的spd轨道杂化。极化率分析发现，易被外加场极化的VB₁₆^-拥有最大的平均极化率和各向异性极化率，不容易被外加场极化的VB₂₄^-拥有最小的平均极化率和各向异性极化率。而且，HOMO-LUMO能隙、平均极化率和各向异性极化率之间存在反比关系。拟合得到的光电子能谱、红外和拉曼谱图有望为以后的实验研究提供帮助。最后，对管状VB₂₀^-的成键分析发现，基态结构上方B₂单元分别与管顶部的3个B原子之间形成2个4c-2e σ键；过渡态结构上方B₂单元分别与管顶部的2个B原子之间形成2个3c-2e σ键，基态与过渡态之间发生流变。此外，芳香性研究发现，VB₂₀^-和VB₂₄^-具有芳香性。热力学性质分析表明，所有研究体系在热力学上都是稳定的，而且定容热容和温度之间存在二次函数关系，标准熵和温度之间存在近似线性关系。

1. [1]
  李林, 王佳伟, 霍俞锦, 孙彩凤, 张瀚月, 张彩凤. 聚乙烯吡咯烷酮保护荧光铜纳米团簇制备及其对乙醇的快速准确检测[J]. 无机化学学报, 2021,37,(12): 2113-2124. doi: 10.11862/CJIC.2021.258LI L, WANG J W, HUO Y J, SUN C F, ZHANG H Y, ZAHGN C F. Preparation of Polyvinylpyrrolidone-Protected Fluorescent Copper Nanoclusters for Rapid and Accurate Detection of Ethanol[J]. Chinese J. Inorg. Chem., 2021, 37(12): 2113-2124. doi: 10.11862/CJIC.2021.258
2. [2]
  郑明, 张红星, 袁福龙, 潘清江. 金红石TiO₂纳米团簇与铀酰相互作用的相对论密度泛函理论计算[J]. 无机化学学报, 2018,34,(5): 874-882. ZHNEG M, ZHANG H X, YUAN F L, PAN Q J. Relativistic DFT Calculations of Interaction between Rutile TiO₂ Nanoparticle Clusters and Uranyl Species[J]. Chinese J. Inorg. Chem., 2018, 34(5): 874-882.
3. [3]
  常贯儒, 鲁信勇, 裴春, 陈龙, 李召. 介孔钛纳米复合团簇应用于光热-光动力疗法[J]. 无机化学学报, 2016,32,(7): 1141-1148. CHANG G R, LU X Y, PEI C, CHEN L, LI Z. Application of Mesoporous Titanium Composite Nanoclusters in Photothermal/Photodynamics Therapy[J]. Chinese J. Inorg. Chem., 2016, 32(7): 1141-1148.
4. [4]
  Alexandorova A N, Zhai H J, Wang L S, Boldyrev A I. Molecular Wheel B₈^2- as a New Inorganic Ligand[J]. Photoelectron Spectroscopy and Ab Initio Characterization of LiB₈^-. Inorg. Chem., 2004, 43: 3552-3554.
5. [5]
  Piazza Z A, Li W L, Romanescu C, Sergeeva A P, Wang L S, Boldyrev A I. A Photoelectron Spectroscopy and Ab Initio Study of B₂₁^-: Negatively Charged Boron Clusters Continue to be Planar at 21[J]. J. Chem. Phys., 2012, 136: 104310. doi: 10.1063/1.3692967
6. [6]
  Piazza Z A, Hu H S, Li W L, Zhao Y F, Li J, Wang L S. Planar Hexagonal B 36 as a Potential Basis for Extended Single-Atom Layer Boron Sheets[J]. Nat. Commun., 2014, 5: 3113. doi: 10.1038/ncomms4113
7. [7]
  Zhai H J, Wang L S, Alexandrova A N, Boldyrev A I, Zakrzewski V G. Photoelectron Spectroscopy and Ab Initio Study of B₃^- and B₄^- Anions and Their Neutrals[J]. J. Phys. Chem. A, 2003, 107: 9319-9328.
8. [8]
  Piazza Z A, Popov I A, Li W L, Pal R, Zeng X C, Boldyrev A I, Wang L S. A Photoelectron Spectroscopy and Ab Initio Study of the Structures and Chemical Bonding of the B₂₅ Cluster[J]. J. Chem. Phys., 2014, 141: 034303. doi: 10.1063/1.4879551
9. [9]
  Chen Q, Zhang S Y, Bai H, Tian W J, Gao T, Li H R, Miao C Q, Mu Y W, Lu H G, Zhai H J, Li S D. Cage-like B₄₁⁺ and B₄₂²⁺: New Chiral Members of the Borospherene Family[J]. Angew. Chem. Int. Ed., 2015, 54: 8160-8164. doi: 10.1002/anie.201501588
10. [10]
  Alexandorova A N, Boldyrev A I, Zhai H J, Wang L S. Photoelectron Spectroscopy and Ab Initio Study of the Doubly Antiaromatic B₆^2- Dianion in the LiB₆^- Cluster[J]. J. Chem. Phys., 2005, 122: 054313. doi: 10.1063/1.1839575
11. [11]
  Alexandorova A N, Boldyrev A I, Zhai H J, Wang L S, Steiner E, Fowler P W. Structures and Bonding in B₆^- and B₆: Planarity and Antiaromaticity[J]. J. Phys. Chem. A, 2003, 107: 1359-1369. doi: 10.1021/jp0268866
12. [12]
  Li W L, Romanescu C, Galeev T R, Piazza Z A, Boldyrev A I, Wang L S. Transition-Metal-Centered Nine-Membered Boron Rings: MB₉ and MB₉^-(M=Rh, Ir)[J]. J. Am. Chem. Soc., 2012, 134: 165-168. doi: 10.1021/ja209808k
13. [13]
  Romanescu C, Galeev T R, Li W L, Boldyrev A I, Wang L S. Aromatic Metal-Centered Monocyclic Boron Rings: Co©B₈^- and Ru©B₉^-[J]. Angew. Chem. Int. Ed., 2011, 123: 9506-9509. doi: 10.1002/ange.201104166
14. [14]
  Galeev T R, Romanescu C, Li W L, Wang L S, Boldyrev A I. Observation of the Highest Coordination Number in Planar Species: Decacoordinated TaB₁₀^- and NbB₁₀^- Anions[J]. Angew. Chem. Int. Ed., 2012, 51: 2143-2147.
15. [15]
  Popov I A, Li W L, Piazza Z A, Boldyrev A I, Wang L S. Complexes between Planar Boron Clusters and Transition Metals: A Photoelectron Spectroscopy and Ab Initio Study of CoB₁₂^- and RhB₁₂^-[J]. J. Phys. Chem. A, 2014, 118: 8098-8105. doi: 10.1021/jp411867q
16. [16]
  Galeev T R, Romanescu C, Li W L, Wang L S, Boldyrev A I. Valence Isoelectronic Substitution in the B₈^- and B₉^- Molecular Wheels by an Al Dopant Atom: Umbrella-like Structures of AlB₇^- and AlB₈^-[J]. J. Chem. Phys., 2011, 135: 104301. doi: 10.1063/1.3625959
17. [17]
  Popov I A, Jian T, Lopez G V, Boldyrev A I, Wang L S. Cobalt-Centred Boron Molecular Drums with the Highest Coordination Number in the CoB₁₆^- Cluster[J]. Nat. Commun., 2015, 6: 8654. doi: 10.1038/ncomms9654
18. [18]
  Pham H T, Nguyen M T. Formation of the MB₁₈ Teetotum Boron Clusters with 4d and 5d Transition Metals, M=Rh, Pd, Ir, Pt[J]. J. Phys. Chem. A, 2019, 123(38): 8170-8178. doi: 10.1021/acs.jpca.9b04078
19. [19]
  Ariyarathna I R. Ground and Electronically Excited States of Main-Group-Metal-Doped B₂₀ Double Ring[J]. J. Phys. Chem., 2022, 126(4): 506-512. doi: 10.1021/acs.jpca.1c08631
20. [20]
  Celaya C A, Buendia F, Miralrio A, Borhon L O P, Beltran M, Nguyen M T, Sansores L E. Structures, Stabilities and Aromatic Properties of Endohedrally Transition Metal Doped Boron Clusters MB₂₂, M=Sc and Ti: A Theoretical Study[J]. Phys. Chem. Phys. Chem., 2020, 22: 8077-8087. doi: 10.1039/D0CP00307G
21. [21]
  Dong X, Jalife S, Espinal A V, Barroso J, Orozco M, Ravell E, Cabellos J L, Liang W Y, Cui Z H, Merino G. Li₂B₂₄: The Simplest Combination for a Three-Ring Boron Tube[J]. Nanoscale, 2019, 11: 2143-2147. doi: 10.1039/C8NR09173K
22. [22]
  Li W L. Romanescu C, Piazza Z A, Wang L S[J]. Geometrical Requirements for Transition-Metal-Centered Aromatic Boron Wheels: The Case of VB₁₀^-. Phys. Chem. Phys. Chem., 2012, 14: 13663-13669.
23. [23]
  Chen Q, Zhao Y X, Jiang L X, Li H F, Chen J J, Zhang T, Liu Q Y, He S G. Thermal Activation of Methane by Vanadium Boride Cluster Cations VB_n⁺ (n=3-6)[J]. Phys. Chem. Phys. Chem., 2018, 20: 4641-4645. doi: 10.1039/C8CP00071A
24. [24]
  Pham H T, Phan P M H, Nguyen M T. Impressive Capacity of the B₇^- and V₂B₇ Clusters for CO₂ Capture[J]. Chem. Phys. Lett., 2019, 728: 186-194. doi: 10.1016/j.cplett.2019.04.087
25. [25]
  Tran V T, Tran Q T. Geometric and Electronic Structures of VB₄^0/+/- Clusters and Reactivity of the Cationic Cluster with Methane from Quantum Chemical Calculations[J]. J. Phys. Chem. A, 2019, 123(42): 9223-9233. doi: 10.1021/acs.jpca.9b08536
26. [26]
  Li C G, Shen Z G, Zhang J, Cui Y Q, Li J J, Xue H Y, Li H F, Ren B Z, Hu Y F. Analysis of the Structures, Stabilities and Electronic Properties of MB₁₆^-(M=V, Cr, Mn, Fe, Co, Ni) Clusters and Assemblies[J]. New. J. Chem., 2020, 44: 5109-5119. doi: 10.1039/C9NJ06335H
27. [27]
  李文静, 李果磊, 马宇, 周丽鸽, 王愿香, 刘珊珊, 李成刚, 崔颖琦, 任保增, 胡燕飞. V_nB_8(n+1)^-1 (n=1-3)团簇的几何结构、电子和光谱特性[J]. 原子与分子物理学报, 2021,38,(4): 042006. LI W J, LI G L, MA Y, ZHOU L G, WANG Y X, LIU S S, LI C G, CUI Y Q, REN B Z, HU Y F. Geometric Structures, Electronic and Spectral Properties of V_nB_8(n+1)^-1 (n=1-3) Clusters[J]. Journal of Atomic and Molecular Physics, 2021, 38(4): 042006.
28. [28]
  Li C G, Li H J, Cui Y Q, Tian H, Shao Q Q, Zhang J, Zhao G, Ren B Z, Hu Y F. A Density Functional Investigation on the Structures, Electronic, Spectral and Fluxional Properties of VB₂₀^- Cluster[J]. J. Mol. Liq., 2021, 339: 116764. doi: 10.1016/j.molliq.2021.116764
29. [29]
  Wang Y C, Lv J, Zhu L, Ma Y M. Crystal Structure Prediction via Particle Swarm Optimization[J]. Phys. Rev. B, 2010, 82: 094116. doi: 10.1103/PhysRevB.82.094116
30. [30]
  Wang Y C, Lv J, Zhu L, Ma Y M. CALYPSO: A Method for Crystal Structure Prediction[J]. Comput. Phys. Commun., 2012, 183: 2063-2070. doi: 10.1016/j.cpc.2012.05.008
31. [31]
  Lv J, Wang Y C, Zhu L, Ma Y M. Particle Swarm Structure Prediction on Clusters[J]. J. Chem. Phys., 2012, 137: 84104. doi: 10.1063/1.4746757
32. [32]
  Perdew J P, Burke K, Ernzerhof M. Generalize Gradient Approximation Made Simple[J]. Phys. Rev. Lett., 1996, 77: 3865-3868. doi: 10.1103/PhysRevLett.77.3865
33. [33]
  Adamo C, Barone V. Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model[J]. J. Chem. Phys., 1999, 110: 6158-6170. doi: 10.1063/1.478522
34. [34]
  Tao J, Perdew J P, Staroverov V N, Scuseria G E. Climbing the Density Functional Ladder: Nonempirical Meta-Generalized Gradient Approximation Designed for Molecules and Solids[J]. Phys. Rev. Lett., 2003, 91: 146401-146405. doi: 10.1103/PhysRevLett.91.146401
35. [35]
  李成刚, 张洁, 申梓刚, 崔颍琦, 任保增, 袁玉全, 胡燕飞. NiB₂₀^-团簇的结构、电子和光谱特性研究[J]. 光学学报, 2020,40,(20): 144-149. LI C G, ZHANG J, SHEN Z G, CUI Y Q, REN B Z, YUAN Y Q, HU Y F. Investigation of Structure, Electronic and Spectral Properties of NiB₂₀^- Cluster[J]. Acta Optica Sinica, 2020, 40(20): 144-149.
36. [36]
  Li C G, Cui Y Q, Tian H, Shao Q Q, Zhang J, Ren B Z, Yuan Y Q. Systematic Investigation of Geometric Structures and Electronic Properties of Lithium Doped Magnesium Clusters[J]. Comput. Mater. Sci., 2021, 200: 110800. doi: 10.1016/j.commatsci.2021.110800
37. [37]
  Li C G, Cui Y Q, Li J X, Guo J S, Cheng L, Ren B Z, Yuan Y Q. Probing the Structural, Electronic and Spectral Properties of a NbB₂₀^- Cluster[J]. Mol. Phys., 2021, 119(10): 1910744. doi: 10.1080/00268976.2021.1910744
38. [38]
  Lu C, Miao M S, Ma Y M. Structural Evolution of Carbon Dioxide Under High Pressure[J]. J. Am. Chem. Soc., 2013, 135: 14167-14171. doi: 10.1021/ja404854x
39. [39]
  Zhao Y R, Xu Y Q, Chen P, Yuan Y Q, Qian Y, Li Q. Structural and Electronic Properties of Medium-Sized Beryllium Doped Magnesium BeMg_n Clusters and Their Anions[J]. Results. Phys., 2021, 26: 104341. doi: 10.1016/j.rinp.2021.104341
40. [40]
  Zhao Y. R, Bai T T, Jia L N, Xin W, Hu Y H, Zheng X S[J]. Probing the Structural and Electronic Properties of Neutral and Anionic Lanthanum-Doped Silicon Clusters. J. Phys. Chem. C, 2019, 123: 28561-28568.
41. [41]
  Jin S Y, Chen B L, Kuang X Y, Lu C, Sun W G, Xia X X, Gutsev G L. Structural and Electronic Properties of Medium-Sized Aluminum-Doped Boron Clusters AlBn and Their Anions[J]. J. Phys. Chem. C, 2019, 123: 6276-6283. doi: 10.1021/acs.jpcc.9b00291
42. [42]
  Peng J F, Jin W Y, Kuang X Y, Lu C. Two-Dimensional Fe₈N Nanosheets: Ferromagnets and Nitrogen Diffusion[J]. J. Phys. Chem. Lett., 2021, 12(35): 8453-8459. doi: 10.1021/acs.jpclett.1c02242
43. [43]
  Chen B L, Sun W G, Kuang X Y, Lu C, Xia X X, Shi H X, Maroulis G. Structural Stability and Evolution of Medium-Sized Tantalum -Doped Boron Clusters: A Half-Sandwich-Structured TaB₁₂^- Cluster[J]. Inorg. Chem., 2018, 57: 343-350. doi: 10.1021/acs.inorgchem.7b02585
44. [44]
  Li Q Y, Xi S G, Hu Y F, Yuan Y Q, Zhao Y R, Li M C, Yuan J J, Yang Y J. Probing the Structural and Electronic Properties of Neutral and Anionic Strontium-Doped Magnesium Clusters[J]. Comput. Mater. Sci., 2021, 197: 110605. doi: 10.1016/j.commatsci.2021.110605
45. [45]
  Lu C, Gong W G, Li Q, Chen C F. Elucidating Stress-Strain of ZrB₁₂ from First-Principles Studies[J]. J. Phys. Chem. Lett., 2020, 11: 9165-9170. doi: 10.1021/acs.jpclett.0c02656
46. [46]
  Shen X Y, Die D, Yang J P, Guo J J. The Ground State Structures and Spectra of Ag₂₀ Clusters and the Adsorption to Carbon Monoxide[J]. Mater. Chem. Phys., 2021, 273: 125134. doi: 10.1016/j.matchemphys.2021.125134
47. [47]
  Tan L P, Die D, Zheng B X. Growth Mechanism, Electronic Properties and Spectra of Aluminum Clusters[J]. Spectrochim. Acta Part A, 2022, 267: 120545. doi: 10.1016/j.saa.2021.120545
48. [48]
  Staroverov V N, Scuseria G E, Tao J, Perdew J P. Comparative Assessment of a New Nonempirical Density Functional: Molecules and Hydrogen-Bonded Complexes[J]. J. Chem. Phys., 2003, 119: 12129. doi: 10.1063/1.1626543
49. [49]
  Yuan Y, Cheng L J. Ferrocene Analogues of Sandwich B₁₂CrB₁₂: A Theory Study[J]. J. Chem. Phys., 2013, 138: 024301. doi: 10.1063/1.4773281
50. [50]
  Nguyen M T, Hung T P, Duong L V, Pham-Ho M P, Nguyen M T. Fullerene-like Boron Clusters Stability by an Endohedrally Doped Iron Atom: B_nFe with n=14, 16, 18 and 20[J]. Phys. Chem. Chem. Phys., 2015, 17: 3000-3003. doi: 10.1039/C4CP04279D
51. [51]
  Pham H T, Nguyen M T. Effects of Bimetallic Doping on Small Cyclic and Tubular Boron Clusters: B₇M₂ and B₁₄M₂ Structures with M=Fe, Co[J]. Phys. Chem. Chem. Phys., 2015, 17: 17335-17345. doi: 10.1039/C5CP01650A
52. [52]
  Wang L S. Photoelectron Spectroscopy of Size-Selected Boron Clusters: From Planar Structures to Borophenes and Borospherenes[J]. Int. Rev. Phys. Chem., 2016, 35: 69-142. doi: 10.1080/0144235X.2016.1147816
53. [53]
  Zhou G D, Duan L Y. Structural Chemistry Basis, Beijing: Peking University Press, 2002: 1-10
54. [54]
  Lu T, Chen F W. Multiwfn: A Multifunctional Wavefunction Analyzer[J]. J. Comput. Chem., 2012, 33: 580-592. doi: 10.1002/jcc.22885
55. [55]
  Yan M, Li H R, Zhao X Y, Lu X Q, Mu Y W, Lu H G, Li S D. Fluxional Bonds in Planar B₁₉^-, Tubular TaB₂₀^-, and Cage-like B₃₉^-[J]. J. Comput. Chem., 2019, 40: 966-970. doi: 10.1002/jcc.25728
56. [56]
  Yan M, Li H R, Tian X X, Mu Y W, Lu H G, Li S D. Fluxional Bonds in Quasi-Planar B₁₈^2- and Half-Sandwich MB₁₈^- (M=K, Rb, Cs)[J]. J. Comput. Chem., 2019, 40: 1227-1232. doi: 10.1002/jcc.25782
57. [57]
  Zhao X Y, Luo X M, Tian X X, Lu H G, Li S D. NiB₁₀, NiB₁₁^-, NiB₁₂ and NiB₁₃⁺ : Half-Sandwich Complexes with the Universal Coordination Bonding Pattern of σ Plus π Double Delocalization[J]. J. Cluster. Sci., 2019, 30: 115-121. doi: 10.1007/s10876-018-1457-4

图 1 VB_2n^- (n=8~12)团簇的基态及亚稳态结构、点群对称、电子组态及PBE0/6-311+G(d)和TPSSH/6-311+G(d)水平下的相对能量

Figure 1 Lowest energy and low-lying structures VB_2n^- (n=8-12) clusters along with the point group symmetry, electronic state, and relative energies at PBE0/6-311+G(d) and TPSSH/6-311+G(d) levels

V and B atoms are marked in blue and green, respectively; All units in eV

下载: 全尺寸图片幻灯片

图 2 VB_2n^- (n=8~12)团簇基态结构的分子轨道能级

Figure 2 Molecular orbital energy levels of the lowest energy structures of VB_2n^- (n=8-12) clusters

下载: 全尺寸图片幻灯片

图 3 PBE0和TPSSH泛函下VB_2n^- (n=8~12)团簇基态结构平均束缚能随尺寸变化

Figure 3 Clusters size dependence on the binding energy of the lowest energy structures of VB_2n^- (n=8-12) clusters

下载: 全尺寸图片幻灯片

图 4 (a) VB₁₆^-、(b) VB₁₈^-、(c) VB₂₀^-、(d) VB₂₂^-及(e) VB₂₄^-团簇基态结构的PES谱图

Figure 4 PES spectra of the lowest energy structures for (a) VB₁₆^-, (b) VB₁₈^-, (c) VB₂₀^-, (d) VB₂₂^-, and (e) VB₂₄^- clusters

下载: 全尺寸图片幻灯片

图 5 (a) VB₁₆^-、(b) V B₁₈^-、(c) VB₂₀^-、(d) VB₂₂^-及(e) VB₂₄^-团簇基态结构的红外光谱和最强峰的振动模式图(插图)

Figure 5 IR spectra and vibrational mode of the strongest peak (Inset) for (a) VB₁₆^-, (b) VB₁₈^-, (c) VB₂₀^-, (d) VB₂₂^-, and (e) VB₂₄^- clusters

下载: 全尺寸图片幻灯片

图 6 (a) VB₁₆^-、(b) VB₁₈^-、(c) VB₂₀^-、(d) VB₂₂^-及(e) VB₂₄^-团簇基态结构的拉曼光谱和最强峰的振动模式图(插图)

Figure 6 Raman spectra and vibrational mode of the strongest peak (Inset) for (a) VB₁₆^-, (b) VB₁₈^-, (c) VB₂₀^-, (d) VB₂₂^-, and (e) VB₂₄^- clusters

下载: 全尺寸图片幻灯片

图 7 VB₂₀^-团簇基态的结构流变键图

Figure 7 Fluxional bonds of the lowest energy structure of VB₂₀^- cluster

3c-2e σ bonds and 4c-2e σ bonds formed atop the B₁₈ double-ring tube are highlighted in orange

下载: 全尺寸图片幻灯片

图 8 VB₂₀^-团簇AdNDP成键图

Figure 8 AdNDP bonding patterns of VB₂₀^- cluster

ON is indicated occupation number and the green ball represents B atoms

下载: 全尺寸图片幻灯片

图 9 VB_2n^- (n=8~12)团簇的(a) C_v和(b) S随温度变化的曲线

Figure 9 Temperature dependence curves of (a) C_v and (b) S for VB_2n^- (n=8-12) clusters

下载: 全尺寸图片幻灯片

表 1 VB_2n^- (n=8~12)团簇基态结构中V原子的净电荷(Q)、NEC、键长和Wiberg能级

Table 1. Atomic net charges (Q) of V atoms, NEC, bond lengths, and Wiberg bond order of the lowest energy structures of VB_2n^- (n=8-12) clusters

VB_2n^-	Q / V		NEC(B): 2s22p1	NEC(V): 3d34s2	Average Wiberg bond order
VB_2n^-	NPA	Hirshfeld	NEC(B): 2s22p1	NEC(V): 3d34s2	V—B	B—B
VB₁₆^-	2.759	0.253	2_s0.67-0.702_p2.17-2.37	3_d3.944_s0.024_p0.01	0.521	0.245
VB₁₈^-	4.293	0.242	2_s0.61-0.702_p2.17-2.37	3_d4.044_s0.014_p0.01	0.488	0.273
VB₂₀^-	2.972	0.257	2_s0.67-0.732_p2.17-2.33	3_d3.944_s0.024_p0.01	0.424	0.199
VB₂₂^-	5.454	0.240	2_s0.53-0.792_p1.38-2.48	3_d4.264_s0.014_p0.01	0.391	0.183
VB₂₄^-	4.925	0.279	2_s0.61-0.632_p2.24-2.54	3_d4.064_s0.224_p0.01	0.345	0.167

下载: 导出CSV

表 2 VB_2n^- (n=8~12)团簇基态对应的HOMO-LUMO能隙和极化率

Table 2. HOMO-LUMO energy gaps and polarizabilities of lowest energy structures of VB_2n^- (n=8-12) clusters

Parameter	VB₁₆^-	VB₁₈^-	VB₂₀^-	VB₂₂^-	VB₂₄^-
ΔE^* / eV	2.08	3.04	2.99	3.28	4.06
α_xx / nm³	2.75×10^-1	2.77×10^-1	2.29×10^-1	2.99×10^-1	3.08×10^-1
α_xy / nm³	0	0	0	0	0
α_xz / nm³	0	0	0	0	0
α_yy / nm³	2.75×10^-1	2.77×10^-1	3.19×10^-1	3.04×10^-1	3.08×10^-1
α_yz / nm³	0	0	0	0	0
α_zz / nm³	1.89×10^-1	2.11×10^-1	3.15×10^-1	2.83×10^-1	3.17×10^-1
α / nm³	2.46×10^-1	2.55×10^-1	2.88×10^-1	2.95×10^-1	3.11×10^-1
α/ nm³	4.34×10^-2	4.03×10^-2	4.11×10^-2	3.85×10^-2	3.73×10^-2
Δα / nm³	8.64×10^-2	6.64×10^-2	8.79×10^-2	1.93×10^-2	9.11×10^-2
^*Energy gap of HOMO-LUMO.

下载: 导出CSV

表 3 VB_2n^- (n=8~12)团簇基态结构的NICS值

Table 3. NICS values of the lowest energy structures of VB_2n^- (n=8-12) clusters

Cluster	NICS(0.5)_zz /ppm	NICS(1.0)_zz /ppm	NICS(1.5)_zz /ppm	NICS(2.0)_zz /ppm
VB₁₆^-	94.69	-127.38	-68.03	-31.80
VB₁₈^-	75.23	-6.68	-22.05	-9.65
VB₂₀^-	-157.12	-123.02	-86.06	-66.68
VB₂₂^-	25.18	-50.21	-27.63	-7.72
VB₂₄^-	-83.47	-66.61	-46.45	-15.92

下载: 导出CSV

表 4 计算得到的VB_2n^- (n=8~12)团簇基态结构热力学参数

Table 4. Calculated thermodynamics parameters of the lowest energy structure of VB_2n^-(n=8-12) clusters

Parameter	ΔH /eV	C_v /(J·mol^-1·K^-1)	S /(J·mol^-1·K^-1)
VB₁₆^-	-84.14	212.82	420.86
VB₁₈^-	-94.98	224.85	406.81
VB₂₀^-	-104.98	270.09	465.57
VB₂₂^-	-115.47	280.06	454.79
VB₂₄^-	-125.48	299.41	460.61

下载: 导出CSV

表 5 不同温度下的定容热容和标准熵

Table 5. Constant heat capacities and standard entropies at different temperatures J·mol^-1·K^-1

T / ℃	VB₁₆^-		VB₁₈^-		VB₂₀^-		VB₂₂^-		VB₂₄^-
T / ℃	C_v	S	C_v	S	C_v	S	C_v	S	C_v	S
100	54.89	278.94	44.16	269.55	64.65	287.04	57.71	280.8	61.99	274.93
200	139.96	347.61	136.66	331.99	180.62	372.70	176.03	361.17	188.09	360.87
300	212.82	420.86	224.85	406.81	270.09	465.57	280.06	454.79	299.41	460.61
400	265.57	493.42	289.74	484.70	333.17	556.40	356.39	550.53	383.00	563.05
500	300.64	558.29	332.47	555.82	375.37	637.09	407.10	637.36	439.34	656.43
600	324.66	616.64	361.39	620.41	404.67	709.49	441.75	716.03	478.10	741.34
700	341.45	669.09	381.42	678.76	425.42	774.53	465.92	787.04	505.23	818.16
800	353.48	716.42	395.68	731.56	440.43	833.24	483.21	851.29	524.68	887.78

下载: 导出CSV

扫一扫看文章

计量

PDF下载量: 2
文章访问数: 1523
HTML全文浏览量: 147

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142