计算化学公社

标题: VASP中如何固定电子态? [打印本页]

作者
Author:
zhongyuabc    时间: 2019-1-9 21:24
标题: VASP中如何固定电子态?
今天看到一篇用vasp计算单线态激子和三线激子的文献,具体计算细节如下图所示。对于其中提到的电子态固定是如何操作的?希望各位前辈能指点迷津,十分感谢!


作者
Author:
卡开发发    时间: 2019-1-10 16:20
应该是通过ISMEAR=-2来进行控制,需要设置FERWE、FERDO,手册上面说的比较少,你可以参考https://wiki.kfki.hu/nano/Easy_m ... ith_FERWE_and_FERDO
作者
Author:
zhongyuabc    时间: 2019-1-10 19:00
卡开发发 发表于 2019-1-10 16:20
应该是通过ISMEAR=-2来进行控制,需要设置FERWE、FERDO,手册上面说的比较少,你可以参考https://wiki.kfki ...

救星呀,非常感谢前辈,我去看看勒
作者
Author:
Adair    时间: 2021-10-18 11:41
卡开发发 发表于 2019-1-10 16:20
应该是通过ISMEAR=-2来进行控制,需要设置FERWE、FERDO,手册上面说的比较少,你可以参考https://wiki.kfki ...

老师 您好,我看了您分享的这个介绍的网页(https://wiki.kfki.hu/nano/Easy_m ... ith_FERWE_and_FERDO),有两点不太明白,想向您请教一下。
1. 在考虑自旋极化(ISPIN=2)的时候,将一个电子由VBM拿到CBM(以下简称为单电子激发吧)的时候,需要设置ISTART=1,这地方读取的WAVECAR是在不考虑单电子激发的情况下优化体系结构(ISPIN=2)时所产生的WAVECAR吗?
2.关于文中提到的参数设置:FERWE = 1020*1.0 1*0.0 1*1.0 1000*0.0 #The syntax is BANDS*OCC,其中“1000*0.0”的意思是CBM对应轨道后剩余的轨道数乘以轨道占据数0吗?也就是(总的轨道数-CBM轨道数)*0
作者
Author:
卡开发发    时间: 2021-10-18 16:58
Adair 发表于 2021-10-18 11:41
老师 您好,我看了您分享的这个介绍的网页(https://wiki.kfki.hu/nano/Easy_manual_occupancy_of_Kohn-S ...

1、建议拿到优化结构后最好自洽一下得到WAVECAR做。
2、1000*0.0表示后面1000个轨道是空的。前面1020*1.0是占据态,之后两个轨道1*0.0和1*1.0位于第1021和1022,这两个的占据被调换了。
作者
Author:
Adair    时间: 2021-10-19 22:42
卡开发发 发表于 2021-10-18 16:58
1、建议拿到优化结构后最好自洽一下得到WAVECAR做。
2、1000*0.0表示后面1000个轨道是空的。前面1020*1. ...

谢谢老师,也就是说读取的WAVECAR是基于不考虑单电子激发的情况下优化出的结构做自洽计算后产生的WAVECAR。
作者
Author:
卡开发发    时间: 2021-10-23 18:30
Adair 发表于 2021-10-19 22:42
谢谢老师,也就是说读取的WAVECAR是基于不考虑单电子激发的情况下优化出的结构做自洽计算后产生的WAVECAR ...

对,建议用这个。
作者
Author:
Adair    时间: 2021-10-27 23:33
卡开发发 发表于 2021-10-23 18:30
对,建议用这个。

谢谢老师
作者
Author:
finefun    时间: 2024-11-8 11:36
老师,想问一下,对于多个k点的情况,是只需要对价带顶(如伽马点)的k点设置为激发态,还是所有k点都要设置为激发态呢




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3