| 查看: 1559 | 回复: 3 | ||
[交流]
如何实现液态纯钼的凝固过程的模拟？
|
|
如题，想要模拟不同冷速下液态纯钼的凝固过程，目前只知钼熔点２８８３Ｋ，知道学长用ＬＡＭＭＰＳ做过液态锌的凝固（可是也没有具体程序），小虫又对这块不太了解，所以求助一下大家，大家尽量知无不言，言无不尽，小虫必当倾囊而尽，答谢大家！！ 下面贴出学长“液态锌凝固模拟”的相关内容，希望大家就“液态钼凝固”给点建议，例如，ＬＡＭＭＰＳ程序如何写等等，不胜感激！！ “将1000 个Zn 原子随机置于一个立方体盒子中，并按照Maxwell-Boltzmann 速率分布随机赋予每个原子初始速度值。为了减少尺寸效应，采用周期性边界条件，选取NTP 系综。选用Gear 预测矫正算法进行积分，运行的时间步长为5fs。模拟计算均从770K 开始(Zn 的熔点为692K)，首先让系统等温运行20000 步使之处 于平衡状态；然后以不同冷却速度让系统从770K 冷却到70K。在每个温度点等温运行2000步，以便计算测量各个温度点的构型。采用20 个不同构型的平均结果，计算每个温度点的径向分布函数和Voronoi 多面体结构。”   [ Last edited by ytudou on 2013-12-22 at 09:55 ] |
回复此楼
★ ★ ★ ★ ★ ★ ★
ytudou(金币+1): 谢谢参与
月只蓝: 金币+3, 感谢指导！ 2014-03-11 18:37:01
ytudou(月只蓝代发): 金币+3, 鼓励交流！ 2014-05-31 21:05:16
ytudou(金币+1): 谢谢参与
月只蓝: 金币+3, 感谢指导！ 2014-03-11 18:37:01
ytudou(月只蓝代发): 金币+3, 鼓励交流！ 2014-05-31 21:05:16
|
据我说知，降温或者升温速率是由三个参数来控制，fix 1 all npt T1 T2 tdamp iso P1 P2 pdamp，另外一个就是“run 步数”，冷却速度的话基本上可以简单的用 （T2-T1）/步数 来计算，你可以试试看，另外我自己做下来tdamp和pdamp这两个参数对凝固的最终构型（bcc，fcc or hcp）影响比较大，最近也正在研究~~~ |
2楼2014-03-11 15:23:31
3楼2014-03-12 18:11:37
4楼2017-02-28 19:12:16