| 查看: 2177 | 回复: 6 | |||
凌云志123金虫 (正式写手)
|
[交流]
VASP MD 模拟Fe的相变
|
|
VASP可以做MD运算，现在打算模拟Fe的相变做练习，我的思路是首先按照实验值建一个5x5x5的supercell，然后直接SMASS=3，TEBEG=TEEND=2500K，跑20000步，看径向分布函数显示为液体结构，开始降温冷却，估计在1800k附近出现能量突变，说明凝固出现。 这种思路对吗？ 我在2500k的模拟结果显示是液体结果，但是能量一直在上升，没有在某一个能量值附近波动，大家看是什么原因？ 降温到1800K附近没有出现能量突变，能量还是像2500K时状态一样，在上升？ 我的问题在哪？ 我的INCAR如下：NPAR=4 ISMEAR = 0 ; SIGMA = 0.1 # insulating NCORE = 4 ISYM = 0 # relax NSW = 10000 ; NWRITE = 0 ALGO = V PREC = L MAXMIX = 50 # electronic degrees ENMAX = 250 ; LREAL = A IBRION = 0 ; SMASS = 3 ; NBLOCK = 50 ; POTIM = 2.0 ; NELMIN = 4 EDIFF = 1E-3 TEBEG = 2500 ; TEEND = 2500 # mixing AMIX = 0.8 BMIX = 2.5001 |
回复此楼» 猜你喜欢
275求调剂
已经有9人回复
-
083000学硕274求调剂
已经有5人回复
-
【求调剂】085601材料工程专硕 | 总分272 |
已经有5人回复
-
275求调剂
已经有13人回复
-
0856，材料与化工321分求调剂
已经有8人回复
-
085701环境工程，267求调剂
已经有12人回复
-
320分，材料与化工专业，求调剂
已经有7人回复
-
085701环境工程求调剂
已经有8人回复
-
295材料工程专硕求调剂
已经有5人回复
-
265求调剂
已经有8人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
凌云志123: 金币+5, 谢谢你的答复 2015-12-01 21:03:28
凌云志123: 金币+10, 谢谢 2016-01-10 16:29:26
小木虫: 金币+0.5, 给个红包，谢谢回帖
凌云志123: 金币+5, 谢谢你的答复 2015-12-01 21:03:28
凌云志123: 金币+10, 谢谢 2016-01-10 16:29:26
|
并没有仔细看,我说说我觉得有问题的地方 1.思路对的，但是问题很多，你算相变有没有体积变化?为何用NVT系综,是否需要用NPT 2.降温过程,MD直接模拟像变过程，非常大的问题是会有一个迟滞过程，就是说液态降温到固态，相变点偏低，固态升温到液态，相变点偏高，所以一般相变用thermodynamic integration算自由能差来看，具体看文献： Sugino, O., & Car, R. (1995). Ab initio molecular dynamics study of first-order phase transitions: melting of silicon. Physical review letters, 74(10), 1823. 3.在2500k的模拟结果显示是液体结果，但是能量一直在上升，没有在某一个能量值附近波动，大家看是什么原因？你有没有看看温度是不是稳定了，我怀疑你温度还没平衡。一般如果你是perfect crystal开始跑，也就是能量最低点开始跑得化，初始的温度很快会先降低到你设置值的一半，记得有个定理说过的，大概意思是体系平衡后平均动能等于平均势能。 |
2楼2015-12-01 19:28:20
凌云志123
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1203.9
- 散金: 454
- 红花: 6
- 帖子: 452
- 在线: 98.7小时
- 虫号: 2073907
- 注册: 2012-10-20
- 性别: GG
- 专业: 金属结构材料
|
针对你给我提的问题，我的考虑如下： 1、选择NVT系综是因为我刚开始学习，有人给我个例子，他说这样可以，但我觉得心里没底，我是在用VASP跑MD，不知道VASP的NPT准确性怎么样，manual上说不确保。 2、我模拟液体Fe的温度是2500K，确保高于熔点1800K，降温迟滞就是过冷度吧，这个就目前的降温过程来看，到1700K（NVT）还是液态结果，没有出现能量突变。关于你提的文献，我要仔细学习下，十分感谢！！！ 3、我看了下温度是在2500K上下波动，没有想能量一样一直上升。我开始的结构就是5x5x5的perfect crystal 直接加温到2500K MD。 4、我最近在测试Si的相同过程，液体结构在MD的时候对POTIM很敏感，我降低了POTIM后，能量就基本稳定了。继续降低能量下降了。 5、在MD时，我们有时候并不知道我们的结构到底体积是多少，尤其是在由温度变化时，像在600K模拟MD C在Fe中的扩散用NVT呢？还是NPT？ 希望与你交流，想你多学习！ 谢谢 |
3楼2015-12-01 21:16:42
4楼2016-01-10 21:29:01
凌云志123
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1203.9
- 散金: 454
- 红花: 6
- 帖子: 452
- 在线: 98.7小时
- 虫号: 2073907
- 注册: 2012-10-20
- 性别: GG
- 专业: 金属结构材料
|
我自给也持怀疑态度，现在我已经基本否定了当初的想法，因为金属凝固还有滞后效应，这样模拟是不可能的应该，但是Phase selection during solidification of liquid magnesium via ab initio molecular dynamics simulations Tekalign T. Debela, X. D. Wang, Q. P. Cao, D. X. Zhang, Jian-Jun Zhu, and J. Z. Jiang 这篇文章介绍了他们的计算，不知道他们怎么搞到，我的压力无法出现他们所说的情况 ps：你是用VASP-MD 还是其他的？ |
5楼2016-01-11 16:38:00
匿名
用户注销 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2932.1
- 散金: 317
- 红花: 3
- 帖子: 539
- 在线: 421.3小时
- 虫号: 0
- 注册: 2007-12-05
- 性别: GG
- 专业: 材料学
6楼2016-01-11 17:08:42
凌云志123
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1203.9
- 散金: 454
- 红花: 6
- 帖子: 452
- 在线: 98.7小时
- 虫号: 2073907
- 注册: 2012-10-20
- 性别: GG
- 专业: 金属结构材料
7楼2016-01-11 20:07:24
5