（AgI）n团簇熔化行为的分子动力学模拟

用遗传算法结合经验势搜索了 (AgI)n(n=3-15) 团簇的可能稳定结构,并用微正则分子动力学方法研究了它们的熔化行为.(AgI)n 团簇的稳定结构主要以四叫元环和六元环相接的笼状结构为主.大多数(AgI)n会在一个较大的温度范围内随温度升高结构不断扭曲,原子间距涨落及动能涨落不断增大,直到在某个温度下熔化,结构变得完全无序.(AgI)6 的结构具有很高的对称性,熔化发生在一个较窄的温度范围.对于(AgI)5,会在熔化前较大的温度范围内发生最稳定结构与能量较高的环状异构体之间的转化,并可能出现负热容现象.     

中介尺度Au纳米团簇熔化的分子动力学模拟

采用分子动力学模拟技术，研究了原子个数为16～8628的 Au纳米团簇的熔化过程.采用 Johnson的EAM (embedded atom method) 模型,模拟结果表明，金属纳米团簇存在一中介尺度区域.对Au纳米团簇而言，当原子个数N >456时，团簇的热力学性质与团簇尺寸呈线性关系，熔化首先从表面开始，逐步向中心区域推进，且满足Tmb－Tmc(N)＝aN(－1/3)的关系.另外，计算了中介区域的团簇的尺寸、熔化温度、表面能、熵、焓等热力学量以及均方根位移（RMSD）等动力学量，为研究纳米团簇提供定量数据.     

小尺寸金属团簇熔化过程的分子动力学模拟

运用分子动力学方法模拟了小尺寸金属团簇的熔化过程, 原子之间的作用采用嵌入原子法(EAM)模型, 计算了均方根键长涨落δ随温度的变化, 以及升温过程中团簇热容的变化. 包含55、56个原子的面心立方(FCC)结构Au团簇的熔化过程是基本相同的. 而同样结构和数目Cu团簇的熔化过程却呈现出不同的趋势. Cu55、Cu56在模拟过程中都出现了FCC结构到二十面体结构的转变. 但由于表面多出了一个原子, Cu56的热容曲线比Cu55多了一个峰, 体系出现了预熔化现象. 这表明小尺寸团簇的固液转变的过程与团簇的原子类型、几何结构和原子数目密切相关.     

(BCO)+n(n=1-12)团簇的结构与稳定性

用密度泛函理论(DFT)的B3LYP方法, 在6-31G*水平上对(BCO)⁺_n(n=1-12)团簇的几何结构、电子结构、振动频率等性质进行了理论研究. 结果表明, (BCO)⁺_n团簇的基态结构均为羰基端配位(μ₁-CO)结构, 且含三元环和五元环数目越多或四元环和六元环的数目越少, 相应的结构越稳定. 能量分析得到, n 为奇数的(BCO)⁺_n团簇比n为偶数的稳定.     

Monte Carlo模拟负载型Au团簇的熔化行为

采用Monte Carlo方法研究了石墨负载Au团簇的熔化行为,着重考察了载体对负载型团簇结构的影响。建立了一个正二十面体的Au团簇和AB堆积的石墨载体,通过记录每一个状态的结构来研究团簇的熔化过程。模拟结果表明,随着温度的升高,Au原子从外到内,逐层熔化,形成二维岛状结构。Au与石墨载体之间的相互作用使得Au原子单层分散在石墨表面,相互作用越强,金属原子越靠近载体。     

(BN)n(n≤12)团簇的结构及成键性质

利用遗传算法和Gastreich提出的经验势函数研究了(BN)_n(n≤12)团簇的可能稳定结构, 并对能量较低的异构体在HF/6-31G(d)水平进行优化, 得到了(BN)_n(n≤12)团簇的线状、蒲扇形、单环、双环、三环和笼状结构, 讨论了各种结构的特征及相对稳定性. 分析了BN团簇中原子的成键性质, 在单环结构中, N原子以sp2杂化成键, B原子以sp杂化成键, 而在节点处B原子以sp2杂化成键. (BN)6是唯一没有张力的单环结构.     

(Pt_nMn)~(±,0)(n=1～5)掺杂团簇结构与磁性的密度泛函研究

采用密度泛函理论(Demity Function Theory)中的B3LYP方法,在Lanl2dz赝势基组水平上对(PtnMn)±,0(n=1～5)团簇的几何构型进行了全优化,并对基态的能级以及磁性进行了研究.结果表明:PtMn掺杂团簇的自旋多重度比较高,这种性质跟纯Mn团簇相似.并且发现一般情况下Mn原子参与成键数越多,结构越稳定,在成键数相同的情况下,成键的平均键长越短越稳定;其次(PtnMn)±,0团簇的所有稳定结构都表现为铁磁性耦合;掺杂一个Mn原子后的团簇磁性大大增强,磁矩主要来源于未满的d壳层电子,且Mn原子上的局域磁矩远大于Pt原子.随着Pt原子个数的增加,Mn原子的局域磁矩变化不大,但团簇的总磁矩渐渐增大.     

溴化钾团簇相变的分子动力学研究

众所周知，团簇的尺寸介于原子或分子与大块物质之间。由于它有较大的表面／体积比而有独特的物理化学性质。团簇在成核，晶体生长，材料科学和纳米器件加工等领域起了至关重要的作用。近年来，计算机模拟已成为物理，化学，生物，天文，地质等领域的重要研究工具。计算机模拟不仅可以     

(AlN)+n和(AlN)-n(n=1～15)团簇的结构与稳定性

用密度泛函理论(DFT)的B3LYP方法,在6-31G*水平上,对(AlN)+n和(AlN)-n(n=1～15)团簇的几何构型、红外光谱和热力学稳定性进行了研究,并对它们的电离能及电子亲和能进行了讨论.结果表明:(AlN)n团簇的电荷状态对簇合物的结构有较大影响,随着n的增大影响逐渐减小;所有平衡结构中不存在Al-Al和N-N键,Al-N键是惟一键型;(AlN)+n和(AlN)-n结构稳定性幻数与(AlN)n相同,即n=2,4,6,…等偶数结构较为稳定;(A1N)10团簇具有最大的电离能(IP=9.14 eV)和最小的电子亲和能(EA=0.19eV),较其他团簇更稳定.     

(HCaN3)n (n＝1～5)团簇结构与性质的密度泛函理论研究

利用密度泛函理论在B3LYP/6-311G*水平上对叠氮化合物(HCaN₃)_n (n＝1～5)团簇各种可能构型进行了几何优化, 预测了各团簇的最稳定结构. 并对最稳定结构的成键特性、电荷分布、振动特性及稳定性进行理论研究. 结果表明, HCaN₃团簇最稳定结构为折线型, (HCaN₃)_n (n＝2～4)团簇最稳定结构为叠氮基端位N原子与Ca原子相互链接形成平面环状结构, (HCaN₃)₅团簇最稳定结构为立体钟形结构. 团簇最稳定结构中金属Ca原子均显示正电性, H原子均显示负电性, 叠氮基中间的N原子显示正电性, 叠氮基两端的N原子显示负电性, 且和Ca原子直接作用的N原子的负电性更强. Ca-N键和Ca-H键为典型的离子键, 叠氮基内N原子之间是共价键. 最稳定结构的IR光谱主要分为3个部分, 其最强振动峰均位于2193～2302 cm^－1段, 振动模式为叠氮基中N-N键的反对称伸缩振动. 叠氮基在团簇和晶体中结构不变, 始终以直线型存在, 说明金属叠氮化合物团簇可以很好地模拟其晶体的局域成键和局域电荷转移等特性. 稳定性分析显示, (HCaN₃)₃团簇相对于其他团簇更稳定.     

高压下MgO熔化特性的分子动力学模拟

Molecular dynamics simulation was used to study the melting of MgO at high pressures. The melting temperature of MgO was accurately obtained at elevated temperature and high pressure after corrections based on the modern theory of melting. The calculated melting curve was compared with the available experimental data and other theoretical results at the pressure range of 0-135 GPa. The corrected melting temperature of MgO is in good agreement with the results from Lindemann melting equation and the twophase simulated results below 15 GPa.     

Molecular Dynamics Investigation of Differences in Melting Behaviors of Cu57 and Cu58 Clusters

Within the framework of the embedded-atom method, we performed molecular-dynamics calculations to investigate the structural transformation during melting of two copper clus- ters containing 57 and 58 atoms. The simulation results reveal how their different structural changes can strongly influence internal energy and radial distribution functions. The local structural patterns of different regions during the temperature increase, determined by atom density profiles, are identified for the melting of each cluster. The simulations show sensi- tivities of the structural changes for these two small size clusters with different structures.     

金属Ni熔化前后结构变化的分子动力学模拟

使用Tight-binding势函数, 对FCC-Ni升温熔化过程的结构变化进行了分子动力学模拟. 在定压条件下模拟得到的Ni的熔点在1850 K与1900 K之间. 计算得到了体系在各温度下的径向分布函数和配位数分布等静态结构信息以及动力学性质. 计算得出的液体Ni的扩散系数在1900 K时约为5.02×10−9 m2•s−1, 与实验数据相符. 对液态体系中FCC短程有序结构可能发生的畸变以及由此导致的H-A键型变化进行了分析, 结合配位体构型搜索和键对分析方法计算了各温度下不同短程有序结构的分布. 计算表明, Ni在熔化之后仍保留有部分晶态短程结构, 但发生了较大的畸变, 同时液态中有少量的缺陷二十面体结构存在. 而液体Ni中大多数的配位体的几何构型介于FCC与缺陷二十面体之间.     

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

采用密度泛函理论中的广义梯度近似(GGA)对(CoCr)n (n=1-5)团簇的几何结构、电子结构和磁性进行了系统的研究, 确定了团簇的基态和亚稳态. 结果表明, CoCr二元合金团簇的基态几何构型呈对称有序排列, 其磁性均呈反铁磁性耦合; 团簇键长和配位数的大小对原子局域磁性有很明显的影响; 受Cr原子的影响, 在(CoCr)4团簇中, 非相邻的Co原子之间呈现反铁磁性耦合.     

SrF2和BaF2熔化温度的分子动力学模拟

利用壳层模型分子动力学方法,考虑萤石结构分子中的预熔化现象,对SrF2和BaF2的分子动力学模拟熔化温度进行修正,获得了高压下SrF2和BaF2的熔化温度．同时给出了300K、0．1MPa-7GPa和10．1MPa-3GPa时SrF2和BaF2的状态方程,与已有研究结果的最大误差分别为0．3％和2．2％．计算所得SrF2和BaF2常压下的熔点与已有的实验结果符合较好．对于SrF2和BaF2分子体积变化和已有的熔化模拟的差别也做了比较和讨论．     

离子液体[emim]Br熔点的分子动力学模拟

运用分子动力学模拟, 采用直接加热法和微正则(NVE)系综法计算离子液体[emim]Br的熔点, 以期获得较好的熔点预测方法. 直接加热法通过分析体系的非键合能、密度、径向分布函数、扩散系数和平动序参数随温度的变化关系判断熔点; NVE系综法则通过获得固液共存体系判断熔点. 直接加热法中, 体系易出现过热问题; NVE系综法则能有效克服过热问题, 是在模拟研究中应优先选择的离子液体熔点预测方法.     

CaF2熔化温度的分子动力学模拟

The molecular dynamics method has been applied to simulate the melting temperatures of CaF 2 at elevated temperature and high pressure and to calculate the P～V equation of state of CaF 2 up to 100 GPa at 300 K. The interatomic potential was taken to be the sum of pairwise additive Coulomb, van der Waals attractions, and repulsive interactions. In addition, the shell model was used in molecular dynamics simulation. The pressure dependence of the melting temperature of CaF 2 was predicted up to 4 GPa. However, in order to account for the superheating melting of the molecular dynamic simulation, the simulated melting temperatures of CaF 2 were corrected by the modern theory of melting. Consequently, the melting temperatures of CaF 2 were accurately obtained at elevated temperature and high pressure. Therefore, it is shown that shell model molecular dynamics simulation at constant pressure indeed provides a useful tool for studying the melting temperatures of other materials under high pressures.