分享一个使用杂化泛函计算稀磁半导体的ppt - 第一性原理 - Vasp&MedeA

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
散金: 114
红花: 4
帖子: 390
在线: 236.9小时
虫号: 582657
注册: 2008-07-25
性别: GG
专业: 凝聚态物性 II ：电子结构

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

10楼: Originally posted by huangll99 at 2012-04-24 08:48:52:
谢指教，我的第一个问题其实是做静态计算时报错，但我把ENCUTFOCK=0和NKRED=2都删掉，就可以算了，但是不知为什么啊

ENCUTFOCK=0

这个设置貌似更改了我遇到了一次不记得了
你提交任务的时候有相关提示

12楼2012-04-24 09:34:48

m2125

13楼2012-04-24 09:41:56

tfl03

木虫 (职业作家)

应助: 24 (小学生)
金币: 3168.5
散金: 1126
红花: 42
帖子: 4452
在线: 1812小时
虫号: 772313
注册: 2009-05-17
性别: GG
专业: 凝聚态物性 II ：电子结构

送鲜花一朵

good ddd

此处回帖者，声名水上流。

14楼2012-04-24 10:17:24

huangll99

木虫 (职业作家)

引用回帖:

14楼: Originally posted by tfl03 at 2012-04-24 10:17:24:
good ddd

15楼2012-04-24 10:17:58

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

1ST强帖: 5
应助: 46 (小学生)
贵宾: 8.842
金币: 7440.4
散金: 9654
红花: 118
沙发: 23
帖子: 5306
在线: 1961.4小时
虫号: 537452
注册: 2008-04-01
性别: GG
专业: 理论和计算化学
管辖: 仿真模拟

先看看~

16楼2012-04-24 10:20:46

techfan

顶！

17楼2012-04-26 10:19:24

buct2010

金虫 (正式写手)

1ST强帖: 10
应助: 27 (小学生)
金币: 2328.6
散金: 306
红花: 11
帖子: 351
在线: 231.2小时
虫号: 1129027
注册: 2010-10-22
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+2, 专家考核, 谢谢回帖！ 2012-04-26 17:28:36

我觉得HSE不能算静态的吧，因为静态计算得到的是电荷密度文件CHGCAR，而HSE的本质是混合计算，即用到电荷密度（DFT），又用到波函数（HF），可以理解为用波函数去修正纯DFT得到的结果。所以我觉得算DOS的时候直接用DFT得到的CHGCAR，先用标准DFT计算DOS，再用HSE重新算一下。这是我的理解，你可以试试。
我用HSE算的是能带，CHGCAR用的是DFT算的，而且发现HSE计算并不会改变CHGCAR文件，因为能带计算和DOS计算有点类似，你可以看一下能带计算的例子。
我参考了版主 WDD880227 的帖子
http://muchong.com/bbs/viewthread.php?tid=4232787&fpage=1
还有版主 liliangfang 的帖子
http://muchong.com/bbs/viewthread.php?tid=3990683&fpage=1

[ Last edited by buct2010 on 2012-4-26 at 15:58 ]

18楼2012-04-26 15:56:20

huangll99

木虫 (职业作家)

应助: 43 (小学生)
金币: 775.9
散金: 2040
红花: 20
沙发: 8
帖子: 4051
在线: 1447.4小时
虫号: 1521828
注册: 2011-12-03
性别: GG
专业: 计算机软件

引用回帖:

18楼: Originally posted by buct2010 at 2012-04-26 15:56:20:
我觉得HSE不能算静态的吧，因为静态计算得到的是电荷密度文件CHGCAR，而HSE的本质是混合计算，即用到电荷密度（DFT），又用到波函数（HF），可以理解为用波函数去修正纯DFT得到的结果。所以我觉得算DOS的时候直接 ...

谢谢您的回帖，我看了那两个帖子，我认为您说的计算DOS的方法是有道理的，我会尝试。

其实我提第一个关于静态计算的问题，主要是想用HSE来计算体系的能量。
我的具体的计算内容是计算缺陷晶体的优化，然后做自洽计算（使用的GGA)
现在我的思路是:首先用GGA做结构优化，其次用GGA作自洽计算，最后用HSE作静态计算（这一步使用上一步由GGA计算得到的WAVECAR)

我已经使用了上面这个思路计算了一个小的体系，并用公式：能隙Egap=E(N+1)+E(N-1)-2E(N) 来计算了能隙，结果是比实验值大3%左右，我想是因为默认HF的成分太多的原因

麻烦您帮忙鉴定一下我的思路是否有问题

[ Last edited by huangll99 on 2012-4-26 at 16:48 ]

19楼2012-04-26 16:38:02

buct2010

金虫 (正式写手)

1ST强帖: 10
应助: 27 (小学生)
金币: 2328.6
散金: 306
红花: 11
帖子: 351
在线: 231.2小时
虫号: 1129027
注册: 2010-10-22
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

19楼: Originally posted by huangll99 at 2012-04-26 16:38:02:
谢谢您的回帖，我看了那两个帖子，我认为您说的计算DOS的方法是有道理的，我会尝试。

其实我提第一个关于静态计算的问题，主要是想用HSE来计算体系的能量。
我的具体的计算内容是计算缺陷晶体的优化，然 ...

能量的话我觉得用标准DFT算没什么问题，因为HSE主要是用来修正电子结构的，而且很多DFT算出来的能量和实验很吻合。我觉得你算能量没有必要用HSE。
如果你非要算的话就算一步弛豫吧，把所有原子固定，用弛豫后的波函数。
带隙还是用能带算出来的准确，个人建议，呵呵~

[ Last edited by buct2010 on 2012-4-26 at 17:11 ]

20楼2012-04-26 17:08:45