预测Au_（13-75）团簇基态结构的启发式算法

针对具有NP难度的团簇结构预测问题,提出启发式求解算法——TP-ISDO作算法.该算法包括两阶段局部搜索、内部操作、表面操作和扰动操作.利用TP-ISDO算法预测了Aun（13≤n≤75）团簇的基态结构,其中Au团簇采用Sutton-Chen势能函数模型描述.实验结果表明,该算法能快速地得到Aun（13≤n≤75）团簇的当前已知最低能量结构.特别是对于Au58团簇,得到了两种新构型,这两种构型都是10面体结构,它们的势能值分别为-15648.5689和-15648.8754能量单位,小于当前已知的最低势能值.     

小尺寸Au-Ag合金团簇稳定结构的模拟研究

采用微正则分子动力学方法模拟了原子数N=80的7种Au-Ag合金团簇（按原子数比例不同），得到团簇势能随温度变化关系．结果表明，降温速度对团簇结构的势能影响很大，冷却速度较慢的情况下，形成的团簇势能较低，结构相对稳定．采用相同的降温速度模拟了N=144的7种Au-Ag合金团簇的这种稳定结构，得出各种原子数比例不同的合金团簇的稳定结构具有一定的相似性．     

Gupta势结合多种算法求解Con(n=2～60)团簇基态的效率

采用Gupta势结合完全随机数方法、无规网格点随机数方法、分子动力学模拟淬火法、分子动力学模拟退火法和遗传算法分别求解了Con(n=2～60)团簇的基态几何结构与能量,通过与标准数据库结果对比分析,结果表明：在所讨论的尺寸范围内,遗传算法的求解效率最好,能找到全部的基态结构,分子动力学模拟退火法效率次之,只在个别团簇尺寸时没找到基态结构; 而对于含原子数在35以下的中小团簇,随机数方法、无规网格点随机数法也都有很好的效率;同时,文中也讨论了运用不同随机数、不同网格点及不同淬火温度对结果的影响.     

Gupta势结合多种算法求解Con（n=2-6O）团簇基态的效率

采用Gupta势结合完全随机数方法、无规网格点随机数方法、分子动力学模拟淬火法、分子动力学模拟退火法和遗传算法分别求解了Con（n=2-60）团簇的基态几何结构与能量，通过与标准数据库结果对比分析，结果表明：在所讨论的尺寸范围内，遗传算法的求解效率最好，能找到全部的基态结构，分子动力学模拟退火法效率次之，只在个别团簇尺寸时没找到基态结构；而对于含原子数在35以下的中小团簇，随机数方法、无规网格点随机数法也都有很好的效率；同时，文中也讨论了运用不同随机数、不同网格点及不同淬火温度对结果的影响．     

小金团簇的基态结构和能量

应用遗传算法和Gupta势对小金团簇的结构进行了优化计算，分析讨论了基态结构和能量随团簇尺寸的变化规律，发展Au13为二十面体以及Ann在n≥16时呈现非晶无序结构的特征。     

14粒子Gupta团簇的3次融化现象

利用徽正则经典分子动力学方法,对具有Gupta势的14粒子团簇的融化过程进行了研究.发现在Gupta势的势参量p和q取值合适的情况下,14粒子团簇的融化表现出了一种新的特性,即团簇的3次融化过程.利用原子等价系数σi(t)对这种过程中的团簇结构变化特征进行了研究.     

利用密度泛函理论计算Li2BeN(N=1～12)团簇的最低能量结构及其电子性质

从第一性原理出发,利用密度泛函理论(DFT)计算了Li2BeN(N=1~12)团簇的最低能量结构及其电子性质.计算结果表明,Li2BeN(N=1~12)团簇最低能量结构的对称性与单一成分的铍团簇相比有所下降,能量二阶有限差分随团簇尺寸的增加出现了奇偶振荡现象.从结构稳定性上来看,N=9是Li2BeN(N=1~12)团簇的一个幻数.当N≥10时Li2BeN团簇已经开始具有部分固相的特点.同时发现了团簇中从Be到Li的电荷转移现象.     

第一性原理研究Au_n(n=2-11)团簇的结构和稳定性

从第一性原理出发,利用密度泛函理论中广义梯度近似对Aun(n=2-11)团簇进行了结构优化、能量和振动频率分析,得到了金团簇最低能量结构及亚稳态结构.计算结果表明,Aun(n=2-11)团簇的基态结构均为平面结构,Aun可通过Aun-1边戴帽一个原子生长而得,归咎于相对论效应引起的sd杂化.综合Aun团簇基态的平均结合能、离解能及二阶能量差分可知,n为偶数的Aun团簇的稳定性相对较高,并且Au6团簇的稳定性在所有团簇中是最突出的.     

AlnN+m(n+m=2～5)团簇的结构与稳定性

用自洽场理论 (HF)在 STO- 3G水平上对 Aln N+m(n+m=2～ 5)团簇的几何构型进行优化和频率计算 .通过对得到的各团簇基态结构的键型、键长以及键角等几何参数进行的分析和讨论表明 Al- N键为主要键型 .对其系统的能量和原子化能的比较和研究发现 ,Al N+m (m>1 )基态结构均为多重度较高的电子态对应的结构 ;在相同原子数时 Aln N+结构最稳定 .     

MgPn(n=1～10)团簇的最低能量结构与稳定性的密度泛函研究

利用杂化密度泛函B3LYP方法在6-311G基组水平上对MgPn(n=1～10)团簇进行了几何构型优化和频率分析,得到了各个尺寸团簇的最低能量结构,并对体系的最邻近键长、平均结合能、裂化能、二阶能量差分进行了计算和分析. 结果表明: MgPn(n=1～10)团簇中Mg原子没有嵌入P团簇中,在n≥4时,MgPn(n=1～10)团簇最低能量结构的几何构型可以看作一个Mg原子直接替代纯磷团簇的P原子获得. 在优化得到的最低能量结构的基础上,发现可以用平均结合能和Mg-P间最邻近键长来共同描述MgPn(n=1～10)团簇稳定性的变化,MgPn(n=1～10)团簇在n=2,4具有较高的稳定性. 与纯磷团簇相比,MgPn(n=1～10)团簇的对称性与稳定性都有所降低.     

平面B_nCP(n=1-6)团簇的结构和稳定性的理论研究

利用密度泛函理论(DFT)对平面B_nCP(n=1-6)的几何结构、相对能量、稳定性以及势能面进行了计算和分析.得到了B_nCP(n=1-6)团簇的最稳定结构.除了BCP团簇是直线型结构,其余的B_nCP团簇为相应的n元环状结构.通过分析平均原子成键能(BE)、能量二次差分(△~2E)和成键能增量(IBE)与B原子数的关系后发现,BCP、B_3CP和B_5CP团簇具有较高的稳定性.此外,B_nCP(n=1-6)团簇的稳定性与离域的π-分子轨道、σ-正切分子轨道和σ-径向分子轨道的相互作用有关.价轨道、ADNDP和ELF分析揭示了B_5CP具有σ-芳香性和π-反芳香性.     

Gupta团簇的2次融化过程及其与结构变化的关系

利用微正则经典分子动力学方法,对具有Gupta势的13粒子团簇的融化过程进行了研究.发现在Gupta势的势参量p和q取值合适的情况下,13粒子团簇的融化过程具有2步融化现象.进一步利用原子等价系数σi(t)的研究还发现,团簇的第1步融化仅涉及表面粒子的异构化过程,而第2步融化中的异构化过程则涉及团簇中所有粒子.     

Mg掺杂硅团簇MgSin(n=1~10)的结构和电子性质研究

运用杂化密度泛函B3LYP理论方法,在6-311+G(d,p)基组水平上对MgSin(n=1~10)多种低能异构体进行优化,获得了各个尺寸下团簇的最低能量结构,并研究MgSin(n=1~10)的结构特点和电子性质.结果表明单个镁原子掺杂硅可得到稳定的二元掺杂团簇.     

中等尺寸ZnSe团簇结构稳定性和电子性质研究

文章研究了中等尺寸的(ZnSe)n(n=24,28,36,48)团簇的结构稳定性和电子性质.通过手工搭建和从体材料中切割得到团簇构型,并参考了其它Ⅱ-Ⅳ族团簇的一些低能结构,使用DMol软件包进行结构优化和能量计算.研究结果表明,中等尺寸的(ZnSe)n团簇,当n=24,28和36时,笼状和管状结构的能量最低;当n=48时,纤锌矿(WZ)体材料结构的能量最低;笼状结构、管状结构具有相对较大的HOMO-LUMO能隙,而洋葱结构和WZ体材料结构的能隙相对较小.     

密度泛函理论对NiBen(n=2-12)团簇结构和性质的研究

采用密度泛函理论中的广义梯度近似对NiBen(n=2-12)团簇的几何构型进行优化,并对能量、频率和磁性进行了计算,同时考虑了电子的自旋多重度。结果表明:团簇的几何结构与自旋多重度有关,最低能量结构的自旋多重度是3和1;n≥4时,较强的p,d杂化导致团簇的磁矩淬灭;与Ben 1(n=2-12)团簇相比,NiBen(n=2-12)团簇的稳定性增强;能隙明显减小,化学活性增强,得出了掺杂Ni原子可以增强团簇稳定性和化学活性的结论;二阶能量差分表明,4,10是团簇的幻数。     

平面的BC_nN(n=1～5)团簇的结构和稳定性

在CCSD(T)/6-311+G(d)∥B3LYP/6-311+G(d)理论水平上系统研究了平面BCnN(n=1～5)团簇的几何结构和稳定性。BCnN(n=1～5)团簇的最低能量结构是线形结构。在BCnN(n=1～5)团簇的较低能量结构中,线形结构比环状结构以及含有支链的环状结构在能量上更有利。BCnN(n=1～5)团簇的最低能量结构的平均成键能、成键能增量、垂直电离能和垂直亲和能的研究结果显示,异构体BCN、BC3B和BC5N比它们相邻的团簇更稳定。     

平面B_nCPt(n=1-5)团簇的结构、稳定性和成键性质的密度泛函理论研究

采用密度泛函理论方法研究了团簇B_nCPt(n=1-5)的结构、相对稳定性及势能面.B_nCPt(n=1-5)的最低能量结构都是环状结构.平均成键能(BE)、成键能增量(IBE)、能量二次差分(Δ~2E)和前线轨道能隙(E_(HOMO)-E_(LUMO))的结果表明B_2CPt团簇相对更稳定.最高占据轨道与最低未占据轨道能级差(E_(HOMO)-E_(LUMO))呈现奇偶振荡现象.基于Mayer键级(MBO)、适应性自然密度分割(AdNDP)和核独立化学位移(NICS(1)_zz)分析揭示了B_4CPt团簇的成键特征和芳香性.另外,通过势能面分析考察了B_4CPt异构体的动力学稳定性.     

(HgTe)_n(n=2～8)团簇的基态结构和电子性质

采用密度泛函理论计算,研究了(HgTe)n(n=2～8)团簇基态结构、结合能、能隙、垂直亲和势和二阶能量差分等性质。在广义梯度近似下,对(HgTe)n各种可能的几何结构进行了优化,预测了各团簇的最稳定结构,并对其电子性质进行了计算。研究结果表明:当n=2、n=3时,团簇的最低能量构型是平面结构;当n≥4时,团簇的最低能量结构是三维结构,并且可以看成是由(HgTe)2和(HgTe)3单元组成。(HgTe)n团簇的能隙、垂直亲和势和二阶能量差分表明:(HgTe)3和(HgTe)6是团簇的稳定结构,即n=3和n=6是团簇的幻数。     

Al_n(n=2～10)团簇结构和性质的密度泛函理论研究

在密度泛函理论框架下,采用广义梯度近似泛函BPW91研究了Al_n(n=2～10)团簇的基态结构,系统计算了平均结合能、能量二阶差分、HOMO-LUMO能隙及Al_n(n=2～10)团簇的磁性,并与已有的相关数据进行了比较.结果表明,Al_n团簇的最稳定结构在原子数为n=6时从平面结构转变为立体结构,Al7团簇较相邻尺寸团簇的稳定性要高,磁性分析得到Al_n团簇的平均每原子磁矩在n=2,4,6时为2μB,当n=1,3,5,7,9时为1μB,与实验结果一致.     

全金属阳离子团簇X+3（X=Au,Ag,Cu）芳香性的理论研究

运用密度泛函理论(B3LYP)、从头算方法(HF,MP2),对全金属阳离子团簇X3+(X=Au,Ag,Cu)的稳定结构、振动频率、原子化能(atomization energy,AE)、最高占据态和最低空轨道能隙(homo-lumo gap)与电子总能量(考虑了零点能ZPE)作了理论计算.全金属阳离子团簇X3+(X=Au,Ag,Cu)的稳定结构是D3h.在此基础上,研究了全金属阳离子团簇X3+(X=Au,Ag,Cu)环状平面结构D3h的两种磁性质:各向异性磁化率(χanis)和核独立位移(NICS).计算结果表明:环状平面结构的全金属阳离子团簇X3+(X=Au,Ag,Cu)具有很强的芳香性特征.对它们的分子轨道分析表明,环状平面全金属阳离子团簇X3+(X=Au,Ag,Cu)具有多重芳香性,两个σ分子轨道和一个π分子轨道对全金属阳离子团簇X3+(X=Au,Ag,Cu)的芳香性起着重要的作用.