计算化学公社

 找回密码 Forget password
 注册 Register

如何求取晶体中沿不同方向空穴和电子的有效质量?

查看数: 39333 | 评论数: 25 | 收藏 Add to favorites 11
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2017-9-20 16:41

正文摘要:

我在用origin处理能带计算的结果,用二次拟合法拟合在导带底和价带顶这两点的有效质量。比如我在G点附近选取几组数据进行拟合得到E对K的二阶导数,从而得到一个有效质量。问题是:1.我看到文献中可以得到在某一点沿 ...

回复 Reply

珊珊来迟 发表于 Post on 2022-4-14 09:35:09
QFF 发表于 2022-3-17 11:26
请问一下,input第二行可以直接取0.01吗?第一行直接用KPOINTS文件里面的reciprocal格式吗?

1.可以。这个值相当于stepsize,数值越小 考虑进去的k空间的范围越小。因为程序一共计算61个点,kcut=0.01*60
2. 对的
QFF 发表于 Post on 2022-3-17 11:26:05
请问一下,input第二行可以直接取0.01吗?第一行直接用KPOINTS文件里面的reciprocal格式吗?
Marineking 发表于 Post on 2021-9-22 11:56:52
ggdh 发表于 2017-9-20 18:37
理论上有效质量是个张量,你对角化这个矩阵,得到三个本征值和三个本征方向。应该就是文献中给的三个值。
...

我与楼主也有类似的问题,

“怎么做出特定某个方向的能带呢?我们平常计算出的不直接是晶体的沿特定K点路径的能带吗?”

怎么把高对称点的路径联系到具体晶体的方向呢?
珊珊来迟 发表于 Post on 2021-2-28 21:44:11
谢谢您的回复!
您第一行讲到的“方向是在笛卡尔坐标系中的”,是不是和POSCAR给出的晶格矢量在一个坐标系下的,也就是以a轴为x方向,y在ab平面内,垂直于x,然后z垂直于xy平面?
ggdh 发表于 Post on 2021-2-28 20:35:09
珊珊来迟 发表于 2021-2-27 22:48
老师您好,请问emc程序输出的结果文件中给出的三个本征方向,对应的坐标系是什么呀?是晶轴的abc方向吗? ...

方向是在笛卡尔坐标系中的。
如果要投影到abc方向
首先写出abc在笛卡尔方向上的归一化矢量。写成列向量,然后排成一行如下
[a b c] (这是一个3*3的矩阵)
假设你的第一个本征方向是0.3 0.2 0.1,用这个行向量乘以刚才那个3*3矩阵就行
[0.3 0.2 0.1]*[ a b c]
代静静 发表于 Post on 2019-6-25 11:11:56
ggdh 发表于 2017-9-21 20:35
1, 设价带顶,或者导带底部的K点的位置。
2. 是的。
3. 在input_file 里面的band number 那里可以改具 ...

您好,请问回答一中,K点的位置时高对称点的坐标还是什么意思,这个位置具体指什么?
binganfe 发表于 Post on 2018-7-19 23:05:06
楼主你好,我也是算出来有效质量有正有负,你找到问题在哪了吗?
珊珊来迟 发表于 Post on 2018-3-26 22:14:53
ggdh 发表于 2017-9-20 18:37
理论上有效质量是个张量,你对角化这个矩阵,得到三个本征值和三个本征方向。应该就是文献中给的三个值。
...

老师请问这三个本征方向的物理意义是什么呀?
这三个方向与晶格矢量有什么关系吗?
珊珊来迟 发表于 Post on 2018-3-8 20:08:05
maoxinxina 发表于 2018-3-8 09:18
我想问下第二个参数step size怎么设置?还有为什么我的不管价带还是导带得到的有效质量都是负值。谢谢楼 ...

我看文献中一般都设的0.01
是不是导带选错了?导带是NELECT/2再加1
maoxinxina 发表于 Post on 2018-3-8 09:18:29
ggdh 发表于 2017-9-21 20:35
1, 设价带顶,或者导带底部的K点的位置。
2. 是的。
3. 在input_file 里面的band number 那里可以改具 ...

我想问下第二个参数step size怎么设置?还有为什么我的不管价带还是导带得到的有效质量都是负值。谢谢楼主。
ggdh 发表于 Post on 2017-10-2 19:43:32
珊珊来迟 发表于 2017-9-30 22:58
老师,这是我的能带图

附件是我的EIGENVAL文件,input及output文件。

94  应该是HOMO把,HOMO的话有效质量应该是负的,结果你算出来是正的。。。
这就肯定有问题
你取的k点是0 0.5 0.5 这对应于你能带图中的哪个点呢?

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
卡开发发 + 1 LZ说是U点:-]

查看全部评分 View all ratings

ggdh 发表于 Post on 2017-9-29 23:16:12
1.最好把能带图贴上来。看看emc的给出的信息是否和能带图一致。然后再结合具体的结构分析。
2 你可以用ms,在相应的分数坐标上放置原子,比如你这里的第三个矢量,就在000和 0.3313,0.3802,1 这两个位置分别放两个原子,然后连起来,就可以看到方向了。
3.在outcar里面找到类似下面的话
NELECT =      78.0000    total number of electrons
然后这里的78/3=39,就是最上面的价带,TVB; 40,就是最下面的导带了,BCB
另外 EIGENVAL 文件中每个能带后面应该是有占据数的,我不知道你怎么没有贴出来。占据为0的都是导带。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 17:30 , Processed in 1.650517 second(s), 32 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list