| 查看: 2063 | 回复: 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 |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有10人回复
-
EST投稿状态问题
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
垃圾破二本职称评审标准
已经有10人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+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
凌云志123
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1203.9
- 散金: 454
- 红花: 6
- 帖子: 452
- 在线: 98.7小时
- 虫号: 2073907
- 注册: 2012-10-20
- 性别: GG
- 专业: 金属结构材料
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
凌云志123
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1203.9
- 散金: 454
- 红花: 6
- 帖子: 452
- 在线: 98.7小时
- 虫号: 2073907
- 注册: 2012-10-20
- 性别: GG
- 专业: 金属结构材料
匿名
用户注销 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2932.1
- 散金: 317
- 红花: 3
- 帖子: 539
- 在线: 421.3小时
- 虫号: 0
- 注册: 2007-12-05
- 性别: GG
- 专业: 材料学
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
凌云志123: 金币+5 2016-01-11 20:07:21
小木虫: 金币+0.5, 给个红包，谢谢回帖
凌云志123: 金币+5 2016-01-11 20:07:21
凌云志123
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1203.9
- 散金: 454
- 红花: 6
- 帖子: 452
- 在线: 98.7小时
- 虫号: 2073907
- 注册: 2012-10-20
- 性别: GG
- 专业: 金属结构材料