应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】Dmol中带隙的问题
141/2返回列表12下一页
查看: 1316  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

宏hong

木虫 (正式写手)

[交流] 【求助】Dmol中带隙的问题 已有3人参与

我在计算完带隙后，在.outmol文件中找到了3个energy gap值，如何确认这3个值哪一个是准确的？
另说一下我的发现，原来在看带密度分析资料时，以为0（费米能级）值，如果有曲线穿过，那么该物体为导体，现在发现不是这样的了，如下图：

该图中，显然有曲线穿过0值，但是仍有带隙。带隙为导带低与价带顶的差值，对应于改图，应该是0值两侧的峰值，左侧峰值为价带顶，右侧峰值为导带底，因此，才会有了.outmol中的带隙，不知我这样理解对不对？
谢谢
回复此楼

» 猜你喜欢
1楼 2010-04-10 10:23:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宏hong

木虫 (正式写手)
木人回答？
回复此楼
2楼2010-04-13 11:06:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)
★ ★
xiaohunhun(金币+2):xiexie 2010-04-13 21:42
宏hong(金币+5):谢谢，学习了 2010-04-14 09:06
从态密度图读gap确实有一定的误差，因为绘制图的时候有一定展宽，好处是可以是曲线更平滑。但此时就不能单纯的根据曲线是否为0来计算gap。从峰值判断是个很好的想法，但实行起来还是有困难的，比如你现在的dos图中Ef附近，展宽导致峰发生重叠，所以很难说清这里的峰位在哪，尤其是Ef附近的小峰，几乎被两侧高峰的展宽所淹没了。虽然画DOS图的时候适当减低smearing可以使峰更加明显，但直接从Dos图取值还是要格外小心。
能带图不存在展宽的影响，所以读值还是比较准确的，但注意的是能带图只是绘制了特殊k点，因此必须保证体系的VB_max和CB_min确实出现在这些特殊点上
另外你说的3个energy gap都是怎么得到的呢？DMOL会直接给出gap的值？
赞一下(1人) 回复此楼
3楼2010-04-13 21:01:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宏hong

木虫 (正式写手)
引用回帖:
Originally posted by acridine at 2010-04-13 21:01:30:
从态密度图读gap确实有一定的误差，因为绘制图的时候有一定展宽，好处是可以是曲线更平滑。但此时就不能单纯的根据曲线是否为0来计算gap。从峰值判断是个很好的想法，但实行起来还是有困难的，比如你现在的dos图中 ...

DMOL中的.outmol文件确实给出了带隙值，而且还是3个，如下：
within the k-points mesh accuracy this solid is an insulator or semiconductor
DFT energy gap:      0.01 eV   indirect   ! upper limit from sampling with scf k-mesh
valence band edge at     0.444  0.444  0.444  in fractional units of recip
conduction band edge at  0.000  0.000  0.000

df              ATOMIC  COORDINATES (au)                  DERIVATIVES (au)
df            x          y          z            x          y          z
df   Ni     0.000000   0.000000   0.000000     0.000000   0.000000   0.000000
df   Ni     0.000000   4.014405   4.014405     0.000000   0.000000   0.000000
df   Ni     4.014405   0.000000   4.014405     0.000000   0.000000   0.000000
df   Ni     4.014405   4.014405   0.000000     0.000000   0.000000   0.000000
df    O     4.014405   4.014405   4.014405     0.000000   0.000000   0.000000
df    O     4.014405   0.000000   0.000000     0.000000   0.000000   0.000000
df    O     0.000000   4.014405   0.000000     0.000000   0.000000   0.000000
df    O     0.000000   0.000000   4.014405     0.000000   0.000000   0.000000
df  binding energy      -1.6347149Ha       -44.48287eV       -1025.820kcal/mol


            Total Energy           Binding E                   Time   Iter
Ef        -6364.025407Ha        -1.6347149Ha                  15.5m      9

Crystal symmetry:  full cubic group                                 
  with inversion symmetry  

Crystal symmetry:  full cubic group                                 
  with inversion symmetry  

     type    Max Gradient   Displacement    Energy change    Time Cycle Hessmin
Opt Cartes    0.000000au    0.000000au      -0.0147963Ha    15.5m   1  0.96881

Optimisation of atomic coordinates is already converged
Lat_sc 1dim   1.01926                        -1.6347149
ck_Out NiO.incoor                                                                                                                                                                                                                                                      

cell volume     517.552
              8.10909836653810    0.00000000000000    0.00000000000000
              0.00000000000000    8.10909836653810    0.00000000000000
              0.00000000000000    0.00000000000000    8.10909836653810
ck_vol    533.234

Integration points, checksum  0      4266     34032     144.000066

            Total Energy           Binding E       Cnvgnce     Time   Iter
Ef        -6364.023295Ha        -1.6326036Ha      1.41E-02    15.8m      1
Ef        -6364.026465Ha        -1.6357731Ha      1.80E-02    16.0m      2
Ef        -6364.024407Ha        -1.6337150Ha      9.93E-03    16.2m      3
Ef        -6364.024686Ha        -1.6339942Ha      9.70E-03    16.5m      4
Ef        -6364.024428Ha        -1.6337365Ha      6.42E-04    16.7m      5
Ef        -6364.024426Ha        -1.6337339Ha      8.18E-04    16.9m      6
Ef        -6364.024443Ha        -1.6337515Ha      5.14E-05    17.1m      7
Ef        -6364.024444Ha        -1.6337518Ha      1.04E-05    17.3m      8
Ef        -6364.024444Ha        -1.6337521Ha      4.18E-06    17.5m      9
Message: SCF converged


Fermi Energy   -0.187190 Ha    -5.094 eV

within the k-points mesh accuracy this solid is an insulator or semiconductor
DFT energy gap:      0.11 eV   indirect   ! upper limit from sampling with scf k-mesh
valence band edge at     0.444  0.444  0.444  in fractional units of recip
conduction band edge at  0.000  0.000  0.000

df              ATOMIC  COORDINATES (au)                  DERIVATIVES (au)
df            x          y          z            x          y          z
df   Ni     0.000000   0.000000   0.000000     0.000000   0.000000   0.000000
df   Ni     0.000000   4.054549   4.054549     0.000000   0.000000   0.000000
df   Ni     4.054549   0.000000   4.054549     0.000000   0.000000   0.000000
df   Ni     4.054549   4.054549   0.000000     0.000000   0.000000   0.000000
df    O     4.054549   4.054549   4.054549     0.000000   0.000000   0.000000
df    O     4.054549   0.000000   0.000000     0.000000   0.000000   0.000000
df    O     0.000000   4.054549   0.000000     0.000000   0.000000   0.000000
df    O     0.000000   0.000000   4.054549     0.000000   0.000000   0.000000
df  binding energy      -1.6337519Ha       -44.45667eV       -1025.215kcal/mol


            Total Energy           Binding E                   Time   Iter
Ef        -6364.024444Ha        -1.6337519Ha                  17.8m     10

Crystal symmetry:  full cubic group                                 
  with inversion symmetry  

Crystal symmetry:  full cubic group                                 
  with inversion symmetry  

     type    Max Gradient   Displacement    Energy change    Time Cycle Hessmin
Opt Cartes    0.000000au    0.000000au       0.0009629Ha    17.8m   1  0.96880

Optimisation of atomic coordinates is already converged
Lat_sc 1dim   1.02945                        -1.6337519
ck_Out NiO.incoor                                                                                                                                                                                                                                                      

cell volume     533.234
              7.94852216126012    0.00000000000000    0.00000000000000
              0.00000000000000    7.94852216126012    0.00000000000000
              0.00000000000000    0.00000000000000    7.94852216126012
ck_vol    502.180

Integration points, checksum  0      4266     34032     144.000033

            Total Energy           Binding E       Cnvgnce     Time   Iter
Ef        -6364.022222Ha        -1.6315304Ha      1.40E-02    18.1m      1
Ef        -6364.023435Ha        -1.6327430Ha      1.83E-02    18.3m      2
Ef        -6364.023740Ha        -1.6330477Ha      9.12E-03    18.5m      3
Ef        -6364.023588Ha        -1.6328966Ha      9.89E-03    18.8m      4
Ef        -6364.023397Ha        -1.6327053Ha      1.06E-03    19.1m      5
Ef        -6364.023402Ha        -1.6327099Ha      9.91E-04    19.3m      6
Ef        -6364.023379Ha        -1.6326875Ha      2.25E-04    19.6m      7
Ef        -6364.023375Ha        -1.6326830Ha      2.48E-05    19.8m      8
Ef        -6364.023376Ha        -1.6326837Ha      9.95E-06    20.1m      9
Message: SCF converged


Fermi Energy   -0.182643 Ha    -4.970 eV

this solid is a metal
with density of states at Fermi energy:       2.99 states/per unit cell and eV

df              ATOMIC  COORDINATES (au)                  DERIVATIVES (au)
df            x          y          z            x          y          z
df   Ni     0.000000   0.000000   0.000000     0.000000   0.000000   0.000000
df   Ni     0.000000   3.974261   3.974261     0.000000   0.000000   0.000000
df   Ni     3.974261   0.000000   3.974261     0.000000   0.000000   0.000000
df   Ni     3.974261   3.974261   0.000000     0.000000   0.000000   0.000000
df    O     3.974261   3.974261   3.974261     0.000000   0.000000   0.000000
df    O     3.974261   0.000000   0.000000     0.000000   0.000000   0.000000
df    O     0.000000   3.974261   0.000000     0.000000   0.000000   0.000000
df    O     0.000000   0.000000   3.974261     0.000000   0.000000   0.000000
df  binding energy      -1.6326837Ha       -44.42760eV       -1024.545kcal/mol


            Total Energy           Binding E                   Time   Iter
Ef        -6364.023375Ha        -1.6326837Ha                  20.4m     10

Crystal symmetry:  full cubic group                                 
  with inversion symmetry  

Crystal symmetry:  full cubic group                                 
  with inversion symmetry  

     type    Max Gradient   Displacement    Energy change    Time Cycle Hessmin
Opt Cartes    0.000000au    0.000000au       0.0010683Ha    20.4m   1  0.96876

Optimisation of atomic coordinates is already converged
Lat_sc 1dim   1.00906                        -1.6326837
Lat_op    -1.6326837    -1.6347149    -1.6337519    -0.0005341     0.0029941     0.1783957
lattice optim (2): stress [GPa], redu c-matrix [GPa]: (f0 =  -1.6347149 Ha  step= 0.03000 )
    -33.74       189.11
eigenvalues  2.99E-03
strain0 =         0.178
predicted minimum energy, [Ha]     -1.6347625     -1.6347149
ck_Out NiO.incoor                                                                                                                                                                                                                                                      

cell volume     502.180
              8.04313331824503    0.00000000000000    0.00000000000000
              0.00000000000000    8.04313331824503    0.00000000000000
              0.00000000000000    0.00000000000000    8.04313331824503
cumul-defo   1.0211   0.0000   0.0000   0.0000   1.0211   0.0000   0.0000   0.0000   1.0211  (rel-vol   106.46 %)
deformation: NiO.os35                                                                                                                                             
ck_vol    520.326

Integration points, checksum  0      4266     34032     144.000056

            Total Energy           Binding E       Cnvgnce     Time   Iter
Ef        -6364.025414Ha        -1.6347226Ha      2.52E-03    20.7m      1
Ef        -6364.025661Ha        -1.6349690Ha      3.19E-03    21.0m      2
Ef        -6364.025413Ha        -1.6347216Ha      1.65E-03    21.2m      3
Ef        -6364.025468Ha        -1.6347764Ha      1.82E-03    21.5m      4
Ef        -6364.025452Ha        -1.6347604Ha      9.25E-05    21.7m      5
Ef        -6364.025451Ha        -1.6347595Ha      3.67E-05    22.0m      6
Ef        -6364.025452Ha        -1.6347605Ha      1.10E-05    22.2m      7
Ef        -6364.025452Ha        -1.6347602Ha      2.17E-06    22.5m      8
Message: SCF converged


Fermi Energy   -0.185977 Ha    -5.061 eV

within the k-points mesh accuracy this solid is an insulator or semiconductor
DFT energy gap:      0.03 eV   indirect   ! upper limit from sampling with scf k-mesh
valence band edge at     0.444  0.444  0.444  in fractional units of recip
conduction band edge at  0.000  0.000  0.000
这个该怎么判断哪个是正确值？
赞一下 回复此楼
4楼2010-04-14 09:05:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)
★ ★
hometownlove(金币+2):谢谢交流 2010-04-14 19:59
你这个结果是在优化晶格吧，三个gap数值分别对应晶格常数为
8.02881（这个值你没有贴出来，应该在贴出来部分的之前）
8.0431333182450
8.10909836653810的情况。
而晶格常数为7.94852216126012 的时候，得到的是金属性
从能量上看应该是最后那个结构能量最低吧，这个是最后得到的结构？
另外这个gap值还是很小的，0.11也许还勉强可以说是半导体，0.01eV这个值也太小了。常温下原子热运动的能量就能使电子越过这个gap吧。
顺便说一下感觉dmol优化晶格的功能还是比较弱的，所以没怎么用过。不如自己拟合状态方程舒服
赞一下(1人) 回复此楼
5楼2010-04-14 19:05:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宏hong

木虫 (正式写手)
引用回帖:
Originally posted by acridine at 2010-04-14 19:05:05:
你这个结果是在优化晶格吧，三个gap数值分别对应晶格常数为
8.02881（这个值你没有贴出来，应该在贴出来部分的之前）
8.0431333182450
8.10909836653810的情况。
而晶格常数为7.94852216126012 的时候，得到的 ...

确实是在优化晶格，得到的结果是计算成功后的结果。我也感觉计算的带隙太小了，实验值是4.0左右，而我做的这个才0.03，感觉很无语，在不加自旋设置时，.outmol中根本没有带隙值，设置后，计算的结果就是这样的了。
可不可以在CASTEP里优化晶格，然后利用DMOL算表面的吸附，这样合理吗？
赞一下 回复此楼
6楼2010-04-14 19:57:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)
★ ★
宏hong(金币+5):好，谢谢了。 2010-04-14 21:32
qli0218(金币+2):谢谢 2010-04-15 09:30
dmol优化立方结构还是可以的，结果应该是很准确的，所以就没有必要再引入用CASTEP了。
至于能隙的问题，这个应该是方法的限制，计算一些过渡金属氧化物的时候需要LDA+U才能正确得到能隙的信心，castep有这个功能，dmol好像还没有吧。所以是否讨论能隙、是否+U，才是楼主现在最需要考虑的
赞一下(1人) 回复此楼
7楼2010-04-14 21:22:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

沉默是金8138

铁杆木虫 (正式写手)

副研究员
★
宏hong(金币+5):谢谢 2010-04-15 09:12
qli0218(金币+1):谢谢 2010-04-15 09:30
宏 你的DOS示意图是导体 没有能隙
我的建议 首先你的DOS精确值 不要迷信所谓的Fine ulatra Fine精度之类的设置
查阅相关文献 看看人家DOS的K点设置是多少
还有计算band structure 我不知道你的结构是不是自己构建的
如果是自己构建的 你的band structure的对称点K值是否设置的正确
如果没有设置 直接用默认值
我想你的问题也就应该在此处
K点设置可以查阅相关资料获得
祝楼主好运 早日解决
赞一下(1人) 回复此楼
Whilethereislife,thereishope!
8楼2010-04-15 08:27:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宏hong

木虫 (正式写手)
引用回帖:
Originally posted by acridine at 2010-04-14 21:22:11:
dmol优化立方结构还是可以的，结果应该是很准确的，所以就没有必要再引入用CASTEP了。
至于能隙的问题，这个应该是方法的限制，计算一些过渡金属氧化物的时候需要LDA+U才能正确得到能隙的信心，castep有这个功能 ...

我在小木虫上也看了许多+U的相关问题，而我计算的也是过渡金属。我也试着利用CASTEP的LDA+U计算，但是LDA+U有个问题，它只能计算Setup下的Energy项，并且我在Properties里勾选了计算带隙的density of state，结果输出文件中根本没有带隙值，为了防止是带隙为0这种情况，我利用Analysis绘制带隙，发现View按钮为灰色，意味着根本不能绘制带隙图。另外一个问题就是，利用+U也不能优化晶格了，所以对于+U很矛盾。
赞一下 回复此楼
9楼2010-04-15 09:11:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宏hong

木虫 (正式写手)
引用回帖:
Originally posted by 沉默是金8138 at 2010-04-15 08:27:50:
宏 你的DOS示意图是导体 没有能隙
我的建议 首先你的DOS精确值 不要迷信所谓的Fine ulatra Fine精度之类的设置
查阅相关文献 看看人家DOS的K点设置是多少
还有计算band structure 我不知道你的结构是不是自己构 ...

K点是我查了文献的，大部分都这么设置，应该不会有问题的。structure是从库中直接导入的，我现在算了键能，Ni磁矩，晶格参数，精度还都可以，唯独这个带隙，很无语了
赞一下 回复此楼
10楼2010-04-15 09:14:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 宏hong 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭