| 查看: 7554 | 回复: 12 | ||
spring5908金虫 (小有名气)
|
[求助]
求助：带电缺陷HSE06总能计算结果匪夷所思，请赐教！ 已有1人参与
|
|
我最近在尝试用HSE06泛函做一些缺陷及掺杂结构的理论计算，结果得到了一些匪夷所思的数据，至今未找到原因，着急得不行，所以把计算结果和计算参数贴出来，请高手指点！ 这个工作是关于带有不同电荷数的缺陷态的理论计算，由于超胞80个原子，用HSE06构型优化太耗时了，所以我用一般的DFT方法进行了全优化，然后用HSE06进行静态计算以得到HSE06方法的总能量。结果，用HSE06泛函得到的不同电荷态的缺陷构型的总能较一般DFT方法差别很大，主要体现在两个方面：一是带电缺陷总能比中性缺陷总能小100多ev，二是不同电荷态之间的总能差变的非常小，几乎可以忽略了！ 在做超胞计算之前，我测试了用HSE06泛函计算了单胞的能带结构，结果是合理的，可见泛函本身的使用并没有问题。我将不同电荷缺陷的总能计算结果比较，以及HSE06泛函进行静态计算的INCAR文件贴在下面，请高手帮忙分析下到底是什么错误导致的，如能赐教，不胜感激啊！！！ q=0: -625.5eV (DFT) and -766.2eV (HSE06) q=-1: -624.2eV (DFT) and -885.1eV (HSE06) q=-2: -619.6eV (DFT) and -885.2eV (HSE06) q=-3: -615.0eV (DFT) and -885.3eV (HSE06) INCAR (HSE06, vasp5.2): PREC=Accurate NWRITE=1 ICHARG=2 ISTART=0 INIWAV=1 LREAL=.TRUE. NELECT=489 EDIFF=1.E-7 NELM=1000 NELMIN=4 NELMDL=-5 ISMEAR=0 SIGMA=0.1 ISPIN=1 LHFCALC=.TRUE. HFSCREEN=0.2 ALGO=Damped TIME=0.4 PRECFOCK=Fast SYMPREC=1E-6 NGXF=144 NGYF=144 NGZF=144 NSW=0 IBRION=-1 ADDGRID=.TRUE. EDIFFG=-0.05 ISIF=0 POTIM=0.2 GGA=PE LSCALAPACK=.TRUE. LPLANE=.TRUE. |
回复此楼» 猜你喜欢
关于filecode
已经有6人回复
-
关于代码变化问题，想知道的进来
已经有8人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】vasp 5.2中用HSE06计算不能进入主循环
已经有7人回复
- 【求助】Hse06计算时KPOINTS出错
已经有30人回复
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
2楼2012-10-23 21:19:55
spring5908
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1339.3
- 散金: 1
- 红花: 3
- 帖子: 130
- 在线: 79小时
- 虫号: 715360
- 注册: 2009-03-05
- 性别: MM
- 专业: 理论和计算化学
3楼2012-10-23 21:30:23
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
4楼2012-10-23 21:40:07
spring5908
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1339.3
- 散金: 1
- 红花: 3
- 帖子: 130
- 在线: 79小时
- 虫号: 715360
- 注册: 2009-03-05
- 性别: MM
- 专业: 理论和计算化学
|
你所说的用力收敛具体是怎么做，可否说的详细一点？ 另外，我还有两个问题。之前的DFT必须用pbe么，还是用pw91也可以？另外，http://muchong.com/html/201205/4421729.html这个帖子中提到HSE不能进行静态计算，并且提到总能用DFT结果即可，你认为合理么？ 谢谢！ |
5楼2012-10-23 21:54:47
liliangfang
荣誉版主 (著名写手)
- 1ST强帖: 7
- 应助: 138 (高中生)
- 贵宾: 0.952
- 金币: 10648.7
- 散金: 4575
- 红花: 42
- 沙发: 3
- 帖子: 1856
- 在线: 825.2小时
- 虫号: 1275010
- 注册: 2011-04-23
- 性别: GG
- 专业: 微/纳机械系统
- 管辖: 第一性原理
6楼2012-10-24 08:17:40
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
7楼2012-10-24 22:15:16
spring5908
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1339.3
- 散金: 1
- 红花: 3
- 帖子: 130
- 在线: 79小时
- 虫号: 715360
- 注册: 2009-03-05
- 性别: MM
- 专业: 理论和计算化学
8楼2012-10-26 01:35:43
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
9楼2012-10-26 10:23:19
spring5908
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1339.3
- 散金: 1
- 红花: 3
- 帖子: 130
- 在线: 79小时
- 虫号: 715360
- 注册: 2009-03-05
- 性别: MM
- 专业: 理论和计算化学
10楼2012-10-26 22:05:28
100