24小时热门版块排行榜    

查看: 2336  |  回复: 16

mazuju028

木虫 (正式写手)

[求助] GW计算流程问题

GW计算在vasp网站上有简单说明,还是有些不明确的地方想和大家交流一下。
先简单罗列一下步骤,
1,较少空带的自洽计算。
System  = Si groundstate occupied orbitals
ISMEAR = 0 ; SIGMA = 0.05  ! small sigma is required to avoid partial occupancies
EDIFF = 1E-8               ! required tight tolerance for groundstate orbitals
2,较多能带的自洽计算。
System  = Si unoccupied orbitals
ALGO = Exact               ! use exact diagonalization of the Hamiltonian
NELM = 1                   ! since we are already converged stop after one step
NBANDS = 96               
ISMEAR = 0 ; SIGMA = 0.05  ! small sigma is required to avoid partial occupancies
LOPTICS = .TRUE.
3,G0W0计算。
System  = Si
NBANDS = 96
ISMEAR = 0 ; SIGMA = 0.05
LOPTICS = .TRUE.
ALGO = GW0 ; NOMEGA = 50
---疑问如下---------------------------------------------------------------------------
1,这三步都应该是在同一目录下进行,所以第二步应该是默认读取第一步的波函数,第三步默认读取第二步的波函数?
2,有关K点设置,是否可以从第一步自洽计算开始就用特殊K点?
3,第二步中NELM=1如果不设置,对最后GW计算的结果会有多大影响?
回复此楼
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ApolloYang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
fzx2008: 金币+2, 谢谢指教 2012-05-25 15:29:34
mazuju028: 金币+4, ★★★很有帮助, 谢谢。 2012-05-25 17:03:55
NELM = 1                   ! since we are already converged stop after one step
看注释,是因为step1进行了自洽运算,故这里设置NELM=1是为了减少计算量

2,有关K点设置,是否可以从第一步自洽计算开始就用特殊K点?
楼主是为了算能带吗? 我曾经试过,特殊的k点一开始可以使用的,但是最后GW的计算仍然得不到能带图。
2楼2012-05-25 10:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

引用回帖:
2楼: Originally posted by ApolloYang at 2012-05-25 10:06:12
NELM = 1                   ! since we are already converged stop after one step
看注释,是因为step1进行了自洽运算,故这里设置NELM=1是为了减少计算量

2,有关K点设置,是否可以从第一步自洽计算开始就 ...

你说第二个,在GW算能带这步直接用手动k点算不了吧,请问你是如何算的
3楼2012-05-25 15:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

引用回帖:
2楼: Originally posted by ApolloYang at 2012-05-25 10:06:12
NELM = 1                   ! since we are already converged stop after one step
看注释,是因为step1进行了自洽运算,故这里设置NELM=1是为了减少计算量

2,有关K点设置,是否可以从第一步自洽计算开始就 ...

确实不能用特殊K点算GW能带啊,即便前面自洽计算加了特殊K点,后面计算GW一样出错。说是:
internal error in KPOINT_IN_FULL_GRID: can not find  ***
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
4楼2012-05-25 17:03:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

引用回帖:
3楼: Originally posted by liliangfang at 2012-05-25 15:57:26
你说第二个,在GW算能带这步直接用手动k点算不了吧,请问你是如何算的...

呵呵,我也没能算出来,设置了特殊K点后,GW计算中总是出错。
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
5楼2012-05-25 17:04:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ApolloYang

木虫 (正式写手)

引用回帖:
3楼: Originally posted by liliangfang at 2012-05-25 15:57:26
你说第二个,在GW算能带这步直接用手动k点算不了吧,请问你是如何算的...

我的意思是无法算,等VASP

[ 发自手机版 http://muchong.com/3g ]
6楼2012-05-26 08:04:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ApolloYang

木虫 (正式写手)

引用回帖:
3楼: Originally posted by liliangfang at 2012-05-25 15:57:26
你说第二个,在GW算能带这步直接用手动k点算不了吧,请问你是如何算的...

我的意思是无法算,等VASP

[ 发自手机版 http://muchong.com/3g ]
7楼2012-05-26 08:05:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
fzx2008: 金币+1, 鼓励交流 2012-05-26 12:41:09
现有的版本是算不了特殊k点的GW的,而且其实第一步和第二步没什么太多的区别,完全可以从第二步直接开始。另外vasp新出了一个GW的赝势库,不知道谁有,估计那个赝势库和GW的特殊k点计算有关。
8楼2012-05-26 11:55:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

引用回帖:
8楼: Originally posted by dxcharlary at 2012-05-26 11:55:50
现有的版本是算不了特殊k点的GW的,而且其实第一步和第二步没什么太多的区别,完全可以从第二步直接开始。另外vasp新出了一个GW的赝势库,不知道谁有,估计那个赝势库和GW的特殊k点计算有关。

依然不可以算的
9楼2012-05-26 16:16:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖


fzx2008: 金币+1, 谢谢回帖! 2012-05-26 19:58:25
引用回帖:
9楼: Originally posted by liliangfang at 2012-05-26 16:16:13
依然不可以算的...

我也发现了,不过其实有个办法,前提是内存足够大的话,可以算一个很密的k点然后在里面挑。但是估计内存需求太大了。
10楼2012-05-26 16:40:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mazuju028 的主题更新
信息提示
请填处理意见