24小时热门版块排行榜    

查看: 1231  |  回复: 5

hutudiangong

新虫 (初入文坛)

[求助] 硅的光学计算结果在低于带隙部分有光吸收 已有2人参与

1.对硅进行光学计算,使用LOPTICS = T 得到结果。可以看到即使考虑了DFT对于带隙的低估,仍然可以看到在低于带隙的地方出现介电常数虚部不为零的情况,硅的光学计算结果在低于带隙部分有光吸收,虽然很小。从VASP手册上的计算公式来看硅的光学计算结果在低于带隙部分有光吸收-1,这里面的delta函数应该是限制了低于带隙的介电常数虚部的出现的,所以对于结算结果很是不解。

2.此公式中将能带分成导带价带两部分,那么程序里边具体是如何判断是价带还是导带的呢?是通过电子占有数吗?这个光学计算模块又是如何对待带隙中的杂质带的,是通过电子占有数将他们分到价带导带中去?

谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。 2014-05-05 23:11:24
注意那个delta函数,想想实际程序怎么处理就明白了。你需要减小smearing并检查结果对k点的收敛性。
2楼2014-05-02 15:40:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。 2014-05-05 23:11:30
ALGO = Exact
NBANDS  = 尽量大点,原始的3倍左右
LOPTICS = .TRUE.
CSHIFT = 0.1    # this is the default
NEDOS = 2000

ISMEAR =  0
SIGMA  =  0.01
3楼2014-05-02 22:08:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hutudiangong

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 锐利的碎片 at 2014-05-02 15:40:16
注意那个delta函数,想想实际程序怎么处理就明白了。你需要减小smearing并检查结果对k点的收敛性。

很有启发,非常感谢。
您所指的减小smearing应该是指减小SIGMA吧?
另外能否明示程序中具体用什么函数实现delta函数的
4楼2014-05-04 15:12:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hutudiangong

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by mywai520 at 2014-05-02 22:08:10
ALGO = Exact
NBANDS  = 尽量大点,原始的3倍左右
LOPTICS = .TRUE.
CSHIFT = 0.1    # this is the default
NEDOS = 2000

ISMEAR =  0
SIGMA  =  0.01

多谢。
关于ISMEAR 取值为0 和 -5,会带来哪些方面的差别,还请不吝赐教。
5楼2014-05-04 15:13:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hutudiangong

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by mywai520 at 2014-05-02 22:08:10
ALGO = Exact
NBANDS  = 尽量大点,原始的3倍左右
LOPTICS = .TRUE.
CSHIFT = 0.1    # this is the default
NEDOS = 2000

ISMEAR =  0
SIGMA  =  0.01

我发现我的算例中,光学计算时如果采用ISMEAR=0,那么NEDOS的设置就没有用了,只能够给出200个点的矩阵元
6楼2014-05-04 15:23:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hutudiangong 的主题更新
信息提示
请填处理意见