Na(n)(5<=n<=10)小团簇相变

采用紧密结合的分子动力学模型,对Na(n)(5<=n<=10)小团簇的键长涨落、势能、热容量等熔化性质在50K到1500K温区进行了模拟研究,结果表明：它们发生两次相变,一种在230K到300K的温度范围内,依次有块体玻璃态转变;一种在550K到870K温度段,依次经历了熔化相变.同时也得到随着团簇体系的减小,势能由下向上排列的曲线,即体系的势能由低变高.     

钠原子团簇的相变

采用紧密结合的分子动力学模型，对Ｎａ４、Ｎａ８、Ｎａ２０的结构特性进行了分析。通过对这些钠原子团簇的键长涨落和径向密度分布的分析，发现在２００Ｋ到２５０Ｋ的温度范围内有从类固相到类液相的改变     

钠原子团簇的碰撞动力学

采用紧密结合的分子动力学模型，对第一个满壳附近（ｎ＝８）钠原子团簇碰撞Ｎａｎ＋Ｎａｎ的动力学性质进行了系统性研究。对有限温度下的多次模拟事件，在各种碰撞参数与不同轰击能量时，双团簇（Ｎａｎ）２的稳定性进行了研究。并发现，当质心系单原子能量为０．０２５ｅＶ时，对中心碰撞，（Ｎａ８）２能够存在３０００ｆｓ，而对周边碰撞，（Ｎａ８）２能够存在１００００ｆｓ；不管对中心碰撞还是对周边碰撞，在较低的轰击能量下，这种暂态双团簇结构具有较长的寿命；在相同的单原子轰击能量和相同的碰撞参数，（Ｎａ８）２和（Ｎａ９）２的动力学寿命没有差别，只是（Ｎａ８）２比（Ｎａ９）２的温度更低。     

钠原子团簇物理学

发展了一个与距离相关的紧密结合分子动力学模型来系统研究钠原子团簇Ｎａｎ的结构、热力学和碰撞动力学性质．计算结果能够很好地符合从头计算法（ａｂｉｎｉｔｉｏｎ）的计算结果和实验结果．同时也探讨了这种有限数目团簇随温度升高可能发生的从类固相到类液相的相变和在碰撞过程中形成双团簇（Ｎａｎ）２的动力学稳定性． Based on the tight binding molecular dynamics, the structural, thermodynamical and collision properties are systematically studied for Nan. Present claculations are successful in reproducing some of the experimental data and in agreement very well with some other ab initio calculations. We also investigate the transition from solid like phase to liquid like phase with increasing temperature for this finite system and the stability of dimer (Nan)2 formed in cluster collisions.     

Nan(5≤n≤10)小团簇相变

采用紧束缚的分子动力学模型,对Nan(5≤n≤10)小团簇的键长涨落、势能、热容量等熔化性质在50 K～1500 K温区进行了模拟研究,结果表明:它们发生两次相变,一种在230 K～300 K的温度范围内,依次有块体玻璃态转变;一种在550 K～870 K温度段,依次经历了熔化相变.同时也得到随着团簇体系的减小,势能由下向上排列的曲线,即体系的势能由低变高.     

钠原子团簇相变的多极形变特征

运用距离相关紧密结合的分子动力学模型,通过无量纲的形变参数S1对简单金属原子团簇Na4、Na8、Na14和Na20相变时的多极形变特征进行了分析,发现S1可以用作表征团簇形变的灵敏探针.     

团簇碰音中双团簇分子（Na9）2的形成研究

采用紧密结合分子动力学模型,对Ｎａ９＋Ｎａ９碰撞系统在质心系中单原子的轰击能量为Ｅｃｍ／ｎ＝０．００６３ｅＶ的对心碰撞进行了研究。从动力学的观战对双团簇分子（Ｎａ９）２的形成给预了证明,结合所谓“「突然冷却”技术,鉴别出所形成的（Ｎａ９）２,并给出了其结构性质。     

理论研究C8H80-（H2O）n（n=1～5）团簇

应用密度泛函B3LYP／6—31＋G（d，p）方法对C8H80-（H2O）n（n=1～5）团簇这种弱相互作用体系进行垒自由度能量梯度优化，得到该系列团簇的稳定蛄构．计算结果表明。在该系列二元团簇中，一方面水分子数目的多少对苯基丙酮分子的结构影响很小，另一方面由于苯基丙酮分子的存在，破坏了团簇中水分子的对称性结构，在团簇内部极力形成O—H—O这样弯曲的有方向性的氢键．对苯基丙酮-水这样结构复杂的团簇，指认光谱的难度非常大，本文只讨论了与C=O有关的振动峰和水分子的对称伸缩振动的最强峰．     

团簇碰撞中双团簇分子(Na_9)_2的形成研究

采用紧密结合分子动力学模型 ,对Na9+Na9碰撞系统在质心系中单原子的轰击能量为Ecm/n =0 .0 0 63eV的对心碰撞进行了研究。从动力学的观点对双团簇分子 (Na9) 2 的形成给予了证明 ,结合所谓“突然冷却”技术 ,鉴别出了所形成的 (Na9) 2 ,并给出了其结构性质。     

硅团簇的结构及生长模式——紧束缚分子动力学：Si11—Si32

采用紧束缚分子动力学模拟硅团簇的结构，通过比较它们的结合能来确定基态结构，最后描绘出不同尺寸所对应的径向分布函数、角分布函数.模拟表明硅团簇在n=27处发生结构转变，从结构图上看，是由扁长结构向近球形结构转变.从径向分布函数图像、键角分布函数图像上也可以得到团簇结构在n=27处发生了变化，结构变得越来越紧密. 关键词： 硅团簇 紧束缚分子动力学 模拟退火 基态结构     

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

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

在小硅化物SiXm（X—H，F，Cl；m＝1─4）中平衡几何和振动频率的从头算研究

本文使用分子轨道从头算方法，在ＨＦ／６—３１Ｇ和ＭＰ２／６—３１Ｇ“水平下，计算了小硅化物ＳｉＸｍ的平衡几何和谐振子振动频率。理论值和实验值的比较表明，对于ＳｉＨｍ分子，实验的键长位于ＨＦ键长和ＭＰ２键长之间，而ＭＰ２振动频率更接近于实验的振动；率。对于ＳｉＣｌｍ分子，无论是键长还是振动频率ＨＦ值和ＭＰ２值都比较接近。对于ＳｉＦｍ分子，处于这两类的中间情况。关联作用随着Ｘ原子的原子序数的增加而减弱。     

Study of (Al2O3)n(Ox) clusters with n = 16 and x = 0, 1, 2 from first principles calculations

The ionic and electronic structure of (Al₂O₃)_n(O_x) clusters with n 16 and x = 0, 1, 2 is studied by means of first principles density functional calculations, norm-conserving pseudopotentials and a numerical atomic basis set. The equilibrium geometries have been determined by total energy minimization, starting with several initial geometries for each cluster size. The trends obtained for the atomic arrangements (structural isomers, coordination numbers, disordered versus ordered structures, etc.) and the electronic properties (binding energies, Homo-Lumo gap and dipole moments) are discussed. For most of the oxidized clusters studied here we find that the Homo-Lumo gap and the magnitude of dipole moment of isomeric species can vary drastically.     

(n,n)-(2n, 0)碳纳米管异质结的扭转力学特性

相近直径的锯齿型和扶手椅型碳纳米管可以共轴组合形成5-7碳环交替出现的柱形对称异质结. 本文利用分子动力学方法研究了直径相近且等长锯齿型和扶手椅型碳纳米管形成的(n, n)-(2n, 0)结在扭转过程中的扭矩和轴向应力随扭转角度的变化规律以及应力传递过程. 研究发现, (n, n)-(2n, 0)结扭转应变在达弹性限度内不会产生轴向应力, 该效应对基于碳纳米管扭转特性的纳米振荡器件的设计具有重要意义.     

Ab initio molecular dynamics study on Ag (n = 4, 5, 6)

Ab initio Molecular Dynamics (MD) method, based on density functional theory (DFT) with planewaves and pseudopotentials, was used to study the stability and internal motion in silver cluster Ag_n, with n =4-6. Calculations on the neutral, cationic and anionic silver dimer Ag₂ show that the bond distance and vibrational frequency calculated by DFT are of good quality. Simulations of Ag₄, Ag₅, and Ag₆ in canonical ensemble reveal distinct characteristics and isomerization paths for each cluster. At a temperature of 800 K, an Ag₄ has no definite structure due to internal motion, while for Ag₅ and Ag₆the clusters maintain the planar structure, with atomic rearrangement observed for Ag₅ but not for Ag₆. At a temperature of 200 K, Ag₄ can exist in two planar structures whilst Ag₅ is found to be stable only in the planar form. In contrast Ag₆ is stable in both planar trigonal and 3D pentagonal structures. Micro-canonical MD simulation was performed for all three clusters to obtain the vibrational density of states (DOS). Received 5 May 1999 and Received in final form 20 August 1999     

DFT study of structural, electronic and vibrational properties of pure (Al2O3)n (n = 9, 10, 12, 15) and Ni-doped (Al2O3)n (n = 9, 10) clusters

The geometrical, electronic and vibrational properties of pure (Al₂O₃)_n (n = 9, 10, 12, 15) clusters and Ni-doped (Al₂O₃)_9-10 clusters are investigated by density functional theory. There are four different Ni-doped (Al₂O₃)₉ clusters and one Ni-doped (Al₂O₃)₁₀ cluster taken into account. Compared with the pure clusters, the Ni-doped (Al₂O₃)_9-10 clusters have narrower HOMO-LUMO energy gaps. The results indicate that the impurity of Ni atom is mainly responsible for the reduction of the HOMO-LUMO energy gap. One characteristic vibration band at about 1030 cm⁻¹ is found in the vibrational frequencies of the Ni-doped (Al₂O₃)_9-10 clusters, which is caused by the asymmetric Al-O-Al stretching vibration. Another band at around 826 cm⁻¹ involving the characteristic vibration of Ni-O bond is in good agreement with experimental results.