| 查看: 1454 | 回复: 1 | |||
[求助]
杂化泛函计算出能带需要的费米能级应该取哪一个？
|
|
在杂化泛函计算能带结构的时候，通常有以下顺序 针对杂化泛函特殊K点的PBE计算——〉 针对杂化泛函特殊K点的HSEscf计算，利用前一步的波函数（ISTART=1）——〉 针对杂化泛函特殊K点的HSEnoscf计算（ISTART=1，ICHARG=11） 经过测试发现，利用第二步或者第三步出的能带结构均可。。没啥太大改变。。 但是。出现与PBE计算不同的是，比如纯半导体的体系。在PBE计算出能带结构后，取PBE静态自洽的或者DOS的费米能级，那么价带顶正好是在能量零点处。。。 而用HSE方法计算能带发现无论取HSEscf还是HSEnoscf步的费米能级都达不到这样。。 后来却发现一个很奇怪的现象，就是利用 普通K点（自动生成）PBE的电荷和波函数，直接进行HSE的非自洽计算 （即ISTART=1，ICHARG=11），在前N步电子步都会出现不收敛的情况 entering main loop N E dE d eps ncg rms ort N E dE rms(c) ROT: 1 -0.507082286949E+02 -0.50708E+02 ROT: 2 -0.506938482460E+02 0.14380E-01 。。。。。。。 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 ROT: ** -0.769335096046E+02 -0.88443E+01 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037469E+02 0.88443E+01 ROT: ** -0.769335096047E+02 -0.88443E+01 ROT: ** -0.680892037470E+02 0.88443E+01 gam= 0.000 g(H,U,f)= 0.949E+00-0.182E+01 0.157E-14 ort(H,U,f) = 0.000E+00 0.000E+00 0.000E+00 SDA: 1 -0.237059317507E+02 -0.23706E+02 0.34776E+00 864 -0.869E+00 0.000E+00 。。。。。。。 ROT: 1 -0.107404641718E+03 0.47400E-01 ROT: 2 -0.107404640662E+03 0.10556E-05 ROT: 3 -0.107404641256E+03 -0.59404E-06 gam= 0.382 g(H,U,f)= 0.108E-05 0.878E-05 0.111E-16 ort(H,U,f) = 0.141E-05 0.598E-05 0.827E-16 DMP: 19 -0.243777224012E+02 -0.40814E-05 -0.50766E-05 864 0.987E-05 0.739E-05 1 F= -.24377722E+02 E0= -.24377722E+02 d E =-.153305E-14 mag= 0.0000 writing wavefunctions 直到最后收敛，这样得到的费米能级却是合理的。。。但是觉得这样的计算不合理啊 HSE应该先用波函数进行自洽的啊。。不知道大家有什么看法？ |
回复此楼» 收录本帖的淘帖专辑推荐
 杂化泛函方面的研究 |
» 猜你喜欢
材料与化工调剂
已经有26人回复
-
一志愿211，化学310分，本科重点双非，求调剂
已经有19人回复
-
求调剂288
已经有7人回复
-
298求调剂
已经有13人回复
-
281求调剂
已经有3人回复
-
265求调剂
已经有11人回复
-
301求调剂
已经有5人回复
-
343求调剂
已经有3人回复
-
272分材料子求调剂
已经有18人回复
-
273求调剂
已经有57人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于用vasp计算铜锌锡硫或铜锌锡硒的能带结构和态密度
已经有3人回复
- CASTEP计算半导体能带
已经有4人回复
- PBE0方法在高斯里怎么输入来确定含有的杂化泛函的比例的？
已经有3人回复
- 大神看这里，VASP能带中费米能级问题
已经有4人回复
- 穿越费米能级的杂质带的性质讨论
已经有13人回复
- CASTEP算出GaAs的band gap为0.082eV，该怎么处理？
已经有7人回复
- 下图中的怎么计算？
已经有4人回复
- 【文献学习】Li,Na掺杂ZnO空穴的杂化密度泛函计算(PRB，2009)
已经有21人回复
- 【专家在线答疑贴之三】杂化泛函计算FAQ
已经有28人回复
- vasp 杂化泛函HSE 它怎么就不run涅？
已经有12人回复
- 计算一种氧化物，光吸收系数计算的光吸收边远大于能带带隙，求解释
已经有24人回复
- 请问态密度是这个样子，有磁性吗？谢谢
已经有14人回复
- 怎样判断掺入的原子起空穴还是载流子作用？
已经有4人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】关于杂化泛函HSE计算
已经有20人回复
xueht987
木虫 (正式写手)
- 应助: 112 (高中生)
- 贵宾: 0.081
- 金币: 3514.2
- 散金: 106
- 红花: 6
- 帖子: 593
- 在线: 273.4小时
- 虫号: 1059317
- 注册: 2010-07-17
- 性别: GG
- 专业: 金属材料表面科学与工程
2楼2012-04-01 09:27:30
