模拟进化法优化Si和Be的团簇结构

模拟生物进化的机理，建立了优化团簇结构的理论程序，并对Ｓｉ和Ｂｅ团簇结构进行了优化计算。并成功地得到了用其它方法已发现的团簇基态结构，还发现了用其它方法尚未得到的团簇基态结构；所得到的团簇基态能量等于或小于已发表的用其它方法所得到的相应团簇的基态能量，说明此方法对于优化团簇基态结构是可行和有效的。     

模拟进化法Si和Be的团簇结构

模拟生物进化的机理，建立了优化团簇结构的理论程序，并对Ｓｉ和Ｂｅ团簇结构进行了优化计算。并成功地得到了用其它方法已发现的团簇基态结构，还发现了用其它方法尚未得到的团簇基态结构；所得到的团簇基态能量等于或小于已发表的用其它方法所得到的相团簇的基有量，说明此方法对于优化团簇基态结构是可行和有效的。     

铍团簇Be_n(n=2-55)的基态能量、结构及其拓扑性质

采用一种第一性原理多体势来描述团簇Ben中原子之间的相互作用,用模拟自然界生物进化的方法系统地优化计算了Ben团簇(n=2-55)的基态能量和几何结构,分析了其拓扑性质,总结了能量和结构随原子数目的变化规律,并将得到的Ben的结果和长程Morse势所描述的金属团簇Mn,以及LennardJones势所描述的惰性气体团簇LJn进行了比较     

ConCu(38-n)(n=0～38)团簇基态结构和熔化行为

基于半经验的Gupta多体势采用遗传算法和分子动力学方法并结合模淬火技术,系统研究了ConCu(38-n)(n=0～38)团簇的基态结构与熔化行为．结果表明：除Co7Cu31团簇基态为类Ih结构外,其它ConCu(38-n)混合团簇的基态结构均是在单质Cu38（Co38）的Oh基态结构基础上的畸变;对于混合团簇的基态结构,随Co原子增加均表现出Co原子先占中心后占表面的特征;通过分析基态团簇二阶差分能和混合能表明Co7Cu31和Co14Cu24为幻数结构团簇;ConCu(38-n)混合团簇在熔化过程中均表现出热容曲线无明显宽峰的异常熔化行为;通过对团簇原子等价指数和淬火结构势能分布图对团簇异常熔化行为进行了分析表征,指出在团簇熔化过程中两种动力学稳定结构（类Oh结构与类Ih结构）之间的相互竞争对团簇熔化行为的重要影响．     

ConCu（38-n）（n=0-38）团簇基态结构和熔化行为北大核心CSCD

基于半经验的Gupta多体势采用遗传算法和分子动力学方法并结合模拟淬火技术,系统研究了ConCu（38-n）（n=0-38）团簇的基态结构与熔化行为．结果表明：除Co7Cu31团簇基态为类IA结构外,其它ConCu（38-n）混合团簇的基态结构均是在单质Cuss（Co38）的Oh基态结构基础上的畸变;对于混合团簇的基态结构,随Co原子增加均表现出Co原子先占中心后占表面的特征;通过分析基态团簇二阶差分能和混合能表明Co7Cu31和Co14Cu24为幻数结构团簇;ConCu（38-n）混合团簇在熔化过程中均表现出热容曲线无明显宽峰的异常熔化行为;通过对团簇原子等价指数和淬火结构势能分布图对团簇异常熔化行为进行了分析表征,指出在团簇熔化过程中两种动力学稳定结构（类Oh结构与类IA结构）之间的相互竞争对团簇熔化行为的重要影响．     

MgnLa(n=2-6)掺杂团簇的密度泛函研究

采用密度泛函(DFT)中的B3LYP方法,选择CEP-31g基组,优化并得到了MgnLa(n=2- 6)小团簇的基态平衡结构,计算出了掺杂团簇的基态结构的平均键长、对称性、费米能级、能级间隙、束缚能、总能的二阶差分．结果表明,随着团簇原子数的增加,La原子更容易趋于团簇表面位置. 与未掺杂Mg团簇相比,La原子多采用取代其稳定基态结构的一个Mg原子形成掺杂团簇的稳定结构. 从有限的结果看,其中 Mg3La和 Mg6La的结构更为稳定．     

铍团簇Bem（n—2—55）的基态能量,结构及其拓扑性质

采用一种第一性原理多体势来描述团簇Ｂｅｎ中原子之间的相互作用,用模拟自然界生物进化方法系统地优化计算了Ｂｅｎ团簇（ｎ＝２－５５）的基态能量和几何结构,分析了其拓扑性质,总结了能量和结构随原子数目的变化规律,并将得到Ｂｅｎ的结果和长程Ｍｏｒｓｅ势所获这的金属团簇Ｍｎ,以及Ｌｅｎｎａｒｄ－Ｊｏｎｅｓ势所描述的惰性气体团簇ＬＪｎ进行了比较     

WnSi02,±(n=1～5)团簇结构与光谱的理论研究

利用密度泛函理论(DFT)在B3LYP/Lanl2dz水平上对WnSi0,2±(n=1～5)团簇的各种可能构型进行了几何结构优化,预测了各团簇的基态结构.并对基态构型的平均结合能、能隙、红外光谱及拉曼光谱进行了系统的理论分析.结果表明:WnSi0,2±(n=2～5)团簇的基态结构及亚稳态结构中的Si-Si原子不成键,W5Si0,2±团簇的基态构型与W7团簇的基态构型相似,W3Si2+团簇的基态构型与W0,5±团簇的基态构型相似;W5Si2-团簇是比较理想的复合材料;WnSi0,2±(n=1～5)团簇的所有IR吸收峰都属于红外光谱中的指纹区的特征吸收区域;除W4Si2-团簇外,其它团簇都具有非刚性特性.     

模拟淬火算法结合Gupta势研究Rhn(n=2~100)团簇的基态特性

采用Gupta多体势结合分子动力学模拟淬火方法及遗传算法、分别求解了Rhn（n=2～100）团簇的最低能几何结构及能量,结果表明：与模拟淬火方法所得最低能结构（视为基态结构）相比,当铑团簇尺寸为60以下时,遗传算法基本可以找到全部基态(除Rh50以外）;但随团簇尺寸增大遗传算法寻找基态结构效率明显下降．通过系统分析淬火结构势能分布图得出模拟淬火算法寻找团簇基态能量（结构）的有效温度区间,同时也进一步说明了该方法在寻找团簇基态结构时与团簇尺寸的非简单依赖关系．     

WnSi20,±（n=1-5）团簇结构与光谱的理论研究

利用密度泛函理论（DFT）在B3LYP/Lanl2dz水平上对WnSi20,±（n=1-5）团簇的各种可能构型进行了几何结构优化,预测了各团簇的基态结构。并对基态构型的平均结合能、能隙、红外光谱及拉曼光谱进行了系统的理论分析。结果表明：WnSi20,±（n=2-5）团簇的基态结构及亚稳态结构中的Si-Si原子不成键,W5Si20,±团簇的基态构型与W7团簇的基态构型相似,W3Si2＋团簇的基态构型与W50,±团簇的基态构型相似;W5Si2－团簇是比较理想的复合材料;WnSi20,±（n=1-5）团簇的所有IR吸收峰都属于红外光谱中的指纹区的特征吸收区域;除W4Si2－团簇外,其它团簇都具有非刚性特性。     

WnSi20,±（n=1-5）团簇结构与光谱的理论研究

利用密度泛函理论（DFT）在B3LYP/Lanl2dz水平上对WnSi20,±（n=1-5）团簇的各种可能构型进行了几何结构优化,预测了各团簇的基态结构。并对基态构型的平均结合能、能隙、红外光谱及拉曼光谱进行了系统的理论分析。结果表明：WnSi20,±（n=2-5）团簇的基态结构及亚稳态结构中的Si-Si原子不成键,W5Si20,±团簇的基态构型与W7团簇的基态构型相似,W3Si2＋团簇的基态构型与W50,±团簇的基态构型相似;W5Si2－团簇是比较理想的复合材料;WnSi20,±（n=1-5）团簇的所有IR吸收峰都属于红外光谱中的指纹区的特征吸收区域;除W4Si2－团簇外,其它团簇都具有非刚性特性。     

Cu13团簇几何和电子结构的密度泛函研究

利用遗传算法结合Gupta原子间相互作用势,采用密度泛函理论,系统研究了带电 团簇的基态与低激发态的几何结构和电子结构并与中性Cu13团簇进行了对比．计算结果表明：对 Cu13（中性及带电）团簇,高对称性几何构型在众多异构中无能量竞争性优势,团簇基态结构皆为非紧致低对称性结构,对 负电Cu13找到一种新的低对称性最低能结构;带电明显影响团簇结构稳定性,带电Cu13团簇与中性Cu13团簇的结构稳定性序列显著不同;基态 Cu13（中性及带电）团簇的具有磁矩最小化效应,而其高对称性结构则有较大磁矩;计算所得Cu13团簇电离能及电子亲和势与实验结果相符．     

第一性原理研究SinB（n=1～12）团簇的稳定性

基于密度泛函理论（DFT），我们研究了SinB（n=1～12）团簇的稳定性．结果表明：SinB的基态构型是在Sin-1B的基态或亚稳态构型上带帽一个Si原子而得到；随着团簇尺寸的增大，B原子逐渐从吸附在Sin团簇的表面位置移动到Sin团簇笼内；掺杂B原子提高了纯硅团簇的稳定性；电子总是从Si向B转移，B原子所带的电荷数不仅与B原子的配住数有关，还与SinB团簇的基态结构密切相关.     

铍小团簇的基态能量、结构及其拓扑性质

采用一种1988年得到的第一性原理多体势来描述团簇Be_n中原子之间的相互作用,用模拟自然生物进化的方法系统地优化计算了Be_n团簇(n=2—25)的基态能量和几何结构,分析了其拓扑性质,总结了能量和结构随原子数目的变化规律,并将得到的Be_n的结构和长程Morse势所描述的金属团簇M_n,以及Lennard-Jones势所描述的惰性气体团簇LJ_n的结构进行了比较分析.结果表明：在n=2—25的尺寸范围内 关键词：     

FeBN (N≤3)团簇的结构与磁性

基于第一性原理,利用密度泛函理论中的广义梯度近似(GGA)对FeBN(N≤3)团簇进行了结构优化、能量和频率的计算,得到了FeBN(N≤3)团簇在不同自旋多重度下的平衡结构并确定了团簇的基态构型.结果表明,FeBN(N≤3)基态团簇的自旋多重度分别为4、3、2.我们对FeBN基态团簇的磁性做了系统的研究,得出了Fe原子磁矩和团簇总磁矩随团簇尺寸增大而减小的结论.     

混合团簇Cu11Fe2的的密度泛函研究

以Cu13非紧致低对称性基态结构及四种13原子高对称性（Ih、Oh、D5h、D3h）密堆积结构为初始构型，通过不等价位原子替换，利用密度泛函理论系统研究了Cu11Fe2混合团簇的稳定性及磁性．结果表明：对四种高对称性结构相应中心位替代构型稳定性均高于双表面位替代构型；Cu11Fe2团簇的基态结构为含一个中心位替代和一个表面位替代的Ih替换结构，这与Cu13的基态结构完全不同；Cu11Fe2无明显近基态同分异构现象；相较Cu13团簇、Cu11Fe2团簇的稳定性及磁性均得到增强．     

四碘甲状腺素团簇结构与光谱性质的密度泛函理论研究

在B3LYP/Lanl2mb基组水平上，利用密度泛函理论(DFT)优化了四碘甲状腺素团簇的几何结构.基于该团簇的几何结构下，其吸收和发射光谱的研究使用相同的基组水平并采用极化连续介质模型(PCM)下用含时密度泛函理论(TDDFT).研究结果表明，优化所得甲状腺素团簇的几何结构对称性为C_1;在基态稳定结构基础上，得出其输运性质，即甲状腺素团簇为p型输运材料;通过含时密度泛函理论，在优化好的基态结构基础上，又计算了它的溶剂效应，进一步得出该分子在水溶剂中的吸收光谱和发射光谱特性.     

(CoMn)n(n=1~5)合金团簇的结构和磁性

采用密度泛函理论中的局域自旋密度近似和广义梯度近似对(CoMn)n(n=1～5)团簇的几何构型进行优化、能量、频率和磁性计算,确定了团簇的基态,对其基态的磁性和电子结构进行了系统研究,并与相对应的一元团簇进行了结构和磁性比较.研究表明,两种方法确定的基态构型基本一致,当n=1～4时,等比CoMn二元合金团簇的几何形状仍与一元团簇相同;(CoMn)3和(CoMn)4团簇出现了磁性双稳态,显示铁磁性和反铁磁性耦合;二元CoMn合金团簇中Co、Mn原子磁性仍能保持一元Co、Mn团簇基态的磁性.     

Morse团簇Mn（n=2-100）基态结构特性的模拟淬火算法研究

本文基于Morse两体势函数采用微正则系综分子动力学模拟淬火方法求解了Morse团簇Mn（n=2~100）的基态几何结构和能量．计算结果准确重复了剑桥数据库中已给出的Mn(n=5~80)团簇的最低能量值,体现出基于微正则系综的分子动力学模拟淬火方法在寻找团簇基态结构上的有效性．通过分析Morse团簇的平均束缚能、二阶差分能和一阶差分能、平均最近邻原子间距及平均配位数,得到n=13,19,23,26,39,46,55,71为团簇幻数数列,发现相较平均配位数的影响、团簇平均最近邻原子间距对Morse团簇基态结构的稳定性基本不产生作用．     

密度泛函理论研究Bn Ni(n=6-12)团簇的结构和磁性

基于第一性原理,用密度泛函理论中的广义梯度近似(generalized gradient approximation,GGA)方法,在充分考虑自旋多重度的前提下,优化并得到了Bn(n=6一12)和BnNi(n=6-12)团簇的平衡构型,按照能量最低原理确定其基态结构.Bn团簇的计算结果与已有的理论结果相一致.当Ni原子掺杂在Bn团簇中,B12Ni团簇的基态结构为平面结构,其余均为三维结构.基态结构的自旋多重度除了n=8以外呈现2,l交替的规律.计算团簇基态结构的平均结合能(Eb)、团簇能量的二阶差分(△2 E)和能隙(HOMO.I.UMO,gap)均表明,l=8为B.Ni(,l=6一12)团簇的幻数,即B8Ni团簇较相邻团簇稳定.计算团簇的磁矩表明B8Ni团簇磁矩最大(2цB),团簇总磁矩和平均磁矩随团簇尺寸增大呈现奇偶振荡趋势且磁矩主要由Ni原子的3d轨道提供.