资源描述：
《研究钯团簇的熔点》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、MonteCarlo实验报告一、项目名称：研究耙I才I簇的熔点二、项目内容概要本实验中使用的是软件MonteLab是使用模特卡罗模拟方法來让算大块物体或团簇的热力学性质的软件。金属I才I簇是由儿个到儿百个原子的聚合体组成的，金属团簇有不同于大块金屈或金屈原子的性质。木实验中的研究对象为Pdl3团簇，有前面对于团簇性质的讨论可以知道，我们并不能拿到固态的团簇，也并没冇讨论液体团簇的意义，因为液体、固体这种常'规的定义对于易于聚合的团簇来说并没有意义。但是，“融化”对于团簇可以有其它的类似意义。Pdl3团簇的稳定构型为二十面体构型，这意味着一个稳定的团簇分子，但是当温
2、度升高时，团簇的构型可能会很大程度地偏离二十而体结构，我们把这样的一个温度范围称为熔点，也就是使团簇结构发生巨变的温度。虽然与一般的固体物质融化从微观角度来看是十分相近的，但是对于团簇分子，其结构巨变的温度一般是一个温度范围而不是一个确定的值。本实验的口的就是确定这一温度范围。在金属I才I簇中，不同于一般的分子化合物，两个团簇相撞将会聚合成更大的团簇而不是弹开，这导致一般对丁团簇的研究是在真空屮或者使用其它的隔离团簇的手段。MonteLab所模拟的是处在特定温度下的时候，采用理想气体模型的单一大小的团簇，但是这种模型其实是有缺陷，因为实际的气体和固体小并不存在这样
3、的形态。但这并不妨碍我们从中研究有益的物理性质。本程序中认为因簇形成了低压，这是为了避免团簇的聚合，此时团簇Z间独立存在，也就意味着它们没冇相撞的可能，所以我只需要模拟一个团簇分了即可。但是此时由于一个团簇存在于空间中而没有特定的压强，因此团簇将在全温度范围内存在自动分裂问题，解决方法是将团簇分子限定在一个特定的区域内，拒绝任何离开该区域的操作，但该方法将会使I才I簇结构与实际略有不符；另外的解决方案是通过定义键长限制来确定是否发生分裂。MonteLab通过以上限定就基本完成了对丁•体系的定义和存在问题的解决，可以用来计算体系的物理性质。当温度上升至熔点时，许多参
4、量，如势能或与质心的平均距离曲线会演化出一个弯曲，我们把它当做熔点。与块状材料不同的是，团簇的熔化发生在一个温度区间而非某个特定的温度。运行一系列的模拟，我们能够获得Epot二f(T),定位曲线的弯曲段，从而确定熔点。融化前的原子团簇，可以从半长轴定义一个原子团簇随着温度升高，原子间开始相互位移，团簇的形状不再固定，可以认为团簇开始融化原子团箴Q碎片OpQoo两个原子间最人距离，如果一个原子违法通过B跨到另一个原子，说明产生了碎片任选一个原子，给它一个位移，计算其转移概率之比，并考虑是否产住碎片，以进行接受和拒绝，大量行走后确定整个原子团簇的结构。三、项目实施方法
5、进入MBACLUSTER文件，编译MonteLab.6.f文件。用到的命令为：f95MonteLab.6.f-oMonteLab.6.exe之后运行编译出来的可执行文件，用到的命令为；./MonteLab・6.exe然后按照inputClustcr文件里面给出的顺序输入命令。Z后可以得到三个文件，分别是Data.pdl050k,Conf.pdl3_1050k,s.r.pdl3_1050ko用gedit命令打开Data,pd1050k文件，得到每走1000步以后得到的数据。查看最终的数据可以得到如下图所示的数据：其中Averagepotentialenergy为最后
6、原子的平均是能,Averagesurfacepotentialenergy为平均表面能，Heatcapacity为热容。将所给出的六个数据记录下来。之后改变输入的温度，从50k开始，一直到2200k,每50k为一个间隔，重复上述步骤，得到一系列温度下的能量和热容的数据。先做岀一条曲线，Z后在热容的最高点附近，即1200k左右，每10k取一个点，将温度分别改为1100k,1110k,1120k,1130k,1140k,1150k,1160k,1170k,1180,1190k,1200k,1210k,1220k,1230k,1240k,1250k,1260k,1270
7、k,1280k,1290k,1300k。最后分别以六个数据为纵坐标，以温度为横坐标做出曲线。得到这些数拯与温度之间的关系。四、项目实施结果得到的六个数据和相应的温度温度(K)EE(finalconfigrationE(surface)E(surfacepotentialCvAvereigemetindistanccto)energyoffinalconfigration)thecentreofmass50-34.548773-34.55351998-31.58447222-31.5452171116.865928872.335853366100-34.476059
8、46-34
