| 查看: 1441 | 回复: 9 | |||
| 当前主题已经存档。 | |||
leoaqu金虫 (正式写手)
|
[交流]
【求助】请指教VASP 总能，DOS，band-structure计算步骤
|
||
|
各位虫友，各位大虾： 我对VASP计算体系总能，计算体系DOS以及band-structure的步骤有些疑问。以ZnO为例，如果我要计算2个Co掺杂的ZnO体系（332的supercell）的磁学性子，那么我就要比较FM和AFM结构的体系总能，所以，我的步骤是这样： 一、对Co-ZnO体系进行结构优化（FM和AFM结构分别优化），K点用小点，比如332。得到CONTCAR 和CHGCAR。 二、增加K点到443，用第一步得到的CONTCAR做位这步的POSCAR，NSW＝0， ICHARG=1。这样得到DOSCAR，CHGCAR和OUTCAR用来做结果分析，即DOSCAR－－DOS，CHGCAR－－spin density， OUTCAR－－系统总能E。 三、还是用第一步得到的CONTCAR 和CHGCAR来做位其实文件，并设置NSW＝0， ICHARG=11， KPIONTS沿某些特殊方向取。这样得到的EIGENVAL用来做band分析。 我一直是以为总能的计算与DOS的计算是一样的，我的方法对吗？有什么需要改进的请不吝赐教！ 我的问题还有，就是我用我的方法，有的时候会遇到： 1.在结构优化的时候，本来是FM很稳定（能量低），而我做完第二步的时候，发现变成很AFM稳定了。 2.同一结构（e.g. FM 排列的Co－ZnO）在优化的时候，即KPT＝332的能量反而比KPT＝443的DOS计算时候得到的能量稳定很多。 谢谢! |
回复此楼» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有13人回复
-
285求调剂
已经有8人回复
-
22408 266求调剂
已经有6人回复
-
264求调剂
已经有8人回复
-
材料工程专业日语生求调剂
已经有9人回复
-
297求调剂
已经有15人回复
-
277、学硕，求调剂 数一104，
已经有12人回复
-
机械专硕273请求调剂
已经有3人回复
-
080100力学316求调剂
已经有3人回复
-
288环境专硕,求调材料方向
已经有39人回复
化学小工
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 贵宾: 0.379
- 金币: 5805.4
- 散金: 200
- 红花: 2
- 帖子: 808
- 在线: 62.1小时
- 虫号: 722185
- 注册: 2009-03-14
- 性别: GG
- 专业: 无机纳米化学
2楼2009-10-08 10:04:13
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
3楼2009-10-08 11:08:19
leoaqu
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1221.8
- 散金: 2
- 红花: 1
- 帖子: 336
- 在线: 47.3小时
- 虫号: 592997
- 注册: 2008-09-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
我有很仔细的看vasp的manu，但是还是不太清楚。我的INCAR是这样的： 一，结构优化 INCAR file： ENCUT=520 ISTART=0 ICHARG=2 INIWAV=1 GGA=PE LWAVE= .FALSE. LCHARG=.TRUE. NSW=100 NELM=100 IBRION=2 ISIF=3 EDIFF=1.0E-04 EDIFFG = -0.01 ISMEAR= 0 SIGMA = 0.1 LORBIT=11 LREAL=Auto KPONITS file： Monkhorst Pack 0 Monkhorst Pack 3 3 2 0 0 0 二 DOS和系统能量 INCAR file： ENCUT=520 ISTART=0 ICHARG=1 INIWAV=1 GGA=PE LWAVE= .FALSE. LCHARG=.TRUE. NSW=0 NELM=100 IBRION=－1 ISIF=0 EDIFF=1.0E-06 ＃EDIFFG = -0.01 ISMEAR= 0 SIGMA = 0.1 LORBIT=11 LREAL=Auto EMAX=20 EMIN=-20 NEDOS=2000 KPONITS file： Monkhorst Pack 0 Monkhorst Pack 4 4 3 0 0 0 请问我的INCAR需要做什么样的改动？再就是能量计算跟这个DOS计算是一样的吗？ 关于KPT的问题，我是觉得对332的supercell而言，332的应该已经可以收敛（文献里面的也是用这个kpt），所以，在DOS计算的时候就增加到443（也是别人跟我说的）。 请帮我看看，谢谢。 |
4楼2009-10-08 12:04:48
leoaqu
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1221.8
- 散金: 2
- 红花: 1
- 帖子: 336
- 在线: 47.3小时
- 虫号: 592997
- 注册: 2008-09-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2009-10-08 12:12:58
化学小工
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 贵宾: 0.379
- 金币: 5805.4
- 散金: 200
- 红花: 2
- 帖子: 808
- 在线: 62.1小时
- 虫号: 722185
- 注册: 2009-03-14
- 性别: GG
- 专业: 无机纳米化学
6楼2009-10-08 22:12:32
化学小工
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 贵宾: 0.379
- 金币: 5805.4
- 散金: 200
- 红花: 2
- 帖子: 808
- 在线: 62.1小时
- 虫号: 722185
- 注册: 2009-03-14
- 性别: GG
- 专业: 无机纳米化学
7楼2009-10-08 22:22:25
leoaqu
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1221.8
- 散金: 2
- 红花: 1
- 帖子: 336
- 在线: 47.3小时
- 虫号: 592997
- 注册: 2008-09-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2009-10-08 22:29:59
leoaqu
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1221.8
- 散金: 2
- 红花: 1
- 帖子: 336
- 在线: 47.3小时
- 虫号: 592997
- 注册: 2008-09-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2009-10-08 22:33:56
化学小工
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 贵宾: 0.379
- 金币: 5805.4
- 散金: 200
- 红花: 2
- 帖子: 808
- 在线: 62.1小时
- 虫号: 722185
- 注册: 2009-03-14
- 性别: GG
- 专业: 无机纳米化学
10楼2009-10-08 22:46:17
