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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » Dmol 计算出的能带结构和态密度图不一致的问题
121/2返回列表12下一页
查看: 3691  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

earthwater

铜虫 (初入文坛)

[求助] Dmol 计算出的能带结构和态密度图不一致的问题 已有1人参与

band structure 和 Dos  见附件图中：
能带图显示有带隙，按理说态密度图中处于禁带的那一段能量对应的态密度应该为零，但是DOS图中不为零，不知道怎么回事。


[ Last edited by earthwater on 2011-9-22 at 11:03 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-22 11:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

【答案】应助回帖

★
earthwater(金币+3): 谢谢关注 2011-09-24 15:03:54
zzy870720z(金币+1): 谢谢指教 2011-09-24 17:25:57
有可能是band计算时选择的高对称点不够，没有涵盖所有的高对称点，所以出现了gap
赞一下(1人) 回复此楼
2楼2011-09-24 13:07:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

earthwater

铜虫 (初入文坛)
引用回帖:
2楼: Originally posted by uuv2010 at 2011-09-24 13:07:05:
有可能是band计算时选择的高对称点不够，没有涵盖所有的高对称点，所以出现了gap

你好，
我还换用castep计算了下，相同结构文件，k-points一样。
结果如下：
castep:

Dmol:


Dmol 的band 和 Dos 不一致
Castep的band和Dos 一致

Dmol和Castep的DOS基本相同，
貌似说明了Dmol对于该体系计算得到的Band不可信，有问题？
赞一下 回复此楼
3楼2011-09-24 15:08:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
★
youzhizhe(金币+1): 谢谢交流。 2011-09-26 16:44:37
计算的什么体系？
我怀疑是计算参数设置的问题。
赞一下(1人) 回复此楼
4楼2011-09-24 16:00:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

earthwater

铜虫 (初入文坛)
uuv2010: 可以同时将做dos和band的输入文件贴上来，已经贴上来的这个文件是用来做单点能地计算的 2011-09-24 23:04:44
引用回帖:
4楼: Originally posted by uuv2010 at 2011-09-24 16:00:16:
计算的什么体系？
我怀疑是计算参数设置的问题。

我的体系是二维多层BN，上下面吸氢。
结构文件见附件。
参数设置如下：
===========================
# Task parameters
Calculate                     energy
Symmetry                      on
# Electronic parameters
Spin_polarization             restricted
Charge                        0
Basis                         dnd
Pseudopotential               none
Functional                    pbe
Aux_density                   octupole
Integration_grid              medium
Occupation                    thermal 0.0050
Cutoff_Global                 3.4000 angstrom
Scf_density_convergence       1.0000e-005
Scf_charge_mixing             0.2000
Scf_iterations                50
Scf_diis                      6 pulay
# Kpoint definition file (intervals/offset):
Kpoints                       file     11 11 1 0.0000 0.0000 0.0000
4BN_AB__H1.kpoints
# Print options
Print                         eigval_last_it
# Calculated properties
==============================
赞一下(1人) 回复此楼
5楼2011-09-24 17:51:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

earthwater

铜虫 (初入文坛)
引用回帖:
4楼: Originally posted by uuv2010 at 2011-09-24 16:00:16:
计算的什么体系？
我怀疑是计算参数设置的问题。

我把三个 input 文件都传上来
赞一下 回复此楼
6楼2011-09-24 23:43:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
★
youzhizhe(金币+1): 谢谢交流。 2011-09-26 16:45:19
引用回帖:
6楼: Originally posted by earthwater at 2011-09-24 23:43:07:
我把三个 input 文件都传上来

从输入文件没有发现问题。您在做band图的时候没有加入剪刀因子吧？
赞一下(1人) 回复此楼
7楼2011-09-25 09:57:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

earthwater

铜虫 (初入文坛)
引用回帖:
7楼: Originally posted by uuv2010 at 2011-09-25 09:57:45:
从输入文件没有发现问题。您在做band图的时候没有加入剪刀因子吧？

嗯，没有。
剪刀因子应该怎么设置呢
赞一下 回复此楼
8楼2011-09-25 12:42:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李佳佳佳

金虫 (正式写手)

Alice

【答案】应助回帖

★
youzhizhe(金币+1): 谢谢交流。 2011-09-26 16:45:28
引用回帖:
2楼: Originally posted by uuv2010 at 2011-09-24 13:07:05:
有可能是band计算时选择的高对称点不够，没有涵盖所有的高对称点，所以出现了gap

真有这种可能吗，我用CASTEP算得时候也是，和DOS图对不上，怎么样能改变这个啊？
赞一下(1人) 回复此楼
停不下的马蹄，滴答滴答~
9楼2011-09-26 13:32:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huti88

铁虫 (初入文坛)

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示! 2011-10-01 09:46:28
看起来是Dmol的band算错了。如果Dmol算的dos是对的，说明电子结构收敛到了基态，按理不会得不到正确的能带结构。没用过ms，我猜要么这个不是你这个体系的能带，要么就是能带被投影分解了（没见过这个）。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

earthwater
10楼2011-10-01 01:23:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 earthwater 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭