24小时热门版块排行榜    

查看: 985  |  回复: 7
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

limaojlu

金虫 (正式写手)

[交流] 【求助】请教关于做分子动力学晶胞大小选择的要求。

请教各位大侠了,
我用VASP做一个材料在室温下300K(远低于熔点)的NVT MD,我选择用两种晶胞,一个是87个原子的大晶胞,另外一个是57个原子的小晶胞,我发现做完MD后,大晶胞的原子只是在平衡位置微动而完全保持晶体原子分布特点,而小晶胞原子却不再保持晶体的原子分布,这一点从他们各自的对关联函数也可以看出来,这时小的晶胞已经有非晶化。请教各位大侠这个是不是由于我的晶胞选择太小所致呢? 在我的体统中材料均含有空位,小的晶胞空位比较集中在边界,大的晶胞空位分布在里面,是不是也跟我的空位分布有关?。是不是做分子动力学之前需要先做一个测试(比如像我上面提到的在室温做一个NVT MD 然后来跟真实情况相比较)从而选择合适大小的晶胞?

请大家多多指教

[ Last edited by aylayl08 on 2010-3-22 at 19:36 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)

引用回帖:
Originally posted by limaojlu at 2009-2-27 18:08:
非常感谢您的宝贵意见。
两个包的含量(空位占所有格点比列):大的9.4%,小的10.93%,我原来的目的是模拟10%的空位均匀分配到整个材料上。按照一片PRL上的文献,他们就是做NVT的MD,从对关联函数的角度,我的大 ...

看你做什么了,比如要是模拟一些性质随温度的变化的时候肯定就有用了
5楼2009-02-27 21:14:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

ustbmars

铜虫 (小有名气)

★ ★ ★ ★
limaojlu(金币+1,VIP+0):非常感谢您的解答,对不起啊我总共只有两个金币。 2-27 12:58
limaojlu(金币+1,VIP+0):还有一个 2-27 18:11
wuchenwf(金币+2,VIP+0):谢谢 O(∩_∩)O 2-27 20:36
不知道你所用的两个晶胞空位含量是否一样? 依你文中所述, 你两种晶胞中空位的分布是不一样的。在300K时,不排除原子通过空位进行长程扩散,如果有较多原子占据空位的话则容易形成非晶?!再者,如果较小的那个单胞如果原子数和空位数目相差不大的话(也与空位分布有关), 则很难说这个单胞时晶态还是非晶态,且径向分布函数的特征与晶态的差别较大。

建议你试试两种晶胞都不含空位,看看较小的那个是否还会非晶化?依我的经验,应该是不会的。
2楼2009-02-26 22:14:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limaojlu

金虫 (正式写手)

★ ★
fegg7502(金币+2,VIP+0):鼓励交流! 2-27 21:10
非常感谢您的宝贵意见。
两个包的含量(空位占所有格点比列):大的9.4%,小的10.93%,我原来的目的是模拟10%的空位均匀分配到整个材料上。按照一片PRL上的文献,他们就是做NVT的MD,从对关联函数的角度,我的大胞的模拟结果更符合这个文献。我是不是能这样理解学长的话呢,不管胞的大小,这个10%空位的分布情况如果不同很有可能导致材料热力学性质差别很大,正如我MD出现的情况一样。所以我会好好分析我的两个体系空位的分部差别。
再多请教师兄一些问题,单纯分析分子动力学某一个瞬态的结构,电子结构,光学性质是不是没有意义?因为一般实验的测量值是反映平衡时状态的平均效果,所以是不是一般需要先将结构求平均然后再分析的她的结构和电子结构呢? 是不是我们理解分子动力学一个时刻的结构,需要结合其运动状态去理解呢?
3楼2009-02-27 18:08:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)


fegg7502(金币+1,VIP+0):鼓励交流!wxDragon 2-27 21:11
你这里的空位是相对于节点来说的还是说类似于那种自由体积的

[ Last edited by wuchenwf on 2009-2-27 at 21:15 ]
4楼2009-02-27 20:35:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见