| 查看: 3761 | 回复: 12 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
zhangguangping木虫 (著名写手)
|
[求助]
K点在TranSIESTA计算中的作用
|
|||
|
在计算TranSIESTA输运的时候，如果横向(x,y方向)使用了周期性边界条件(这儿指的是真正的电极周期性，而不是程序包默认的周期性边界条件)，那么就需要在对应的二维布里渊区要取k点。 首先，根据固体物理我们知道，对于周期体系，所有的物理量都是k点的函数。所以要想计算一个物理量，需要对k点求平均。原则上k点需要取无穷多个才能计算准确。但是，考虑到数值计算时间和精度的容忍度，我们可以选择有限的k点来近似计算这些物理量。因此，k点越多，计算结果越接近准确值。 在TranSIESTA中，体系的哈密顿对k点收敛很快，而透射谱的计算对k点收敛却很慢。如图FIG1是分别采用4*4 k点和20*20 k点的结果。从计算结果可以看出，在费米能级上方附近，4*4的k点有很多尖锐的透射峰。当取20*20之后，这个情况得到缓解。但是如果对20*20的再继续放大，如FIG2,就会发现，这种锯齿状的波动还是存在。 问题1：这种锯齿状的情况以前从来没有碰到过，并且对于以前的体系，采用较少的4*4的k点就很平滑了。想请教大家，这种锯齿状的透射谱，能说明什么呢？为什么随着k点的收敛这么慢？这个能给出一些什么信息呢？ 问题2：对于计算透射谱的时候，这个k点什么含义呢？为什么要对k点求平均呢？ 对于某一个能量E，要对所有取的k点求平均。这个k点对应电子的能量吗？如果对应的话，是否对应一个k点（k_x,k_y）有E=(k_x*k_x+k_y*k_y)/2m吗？那k_z怎么处理呢？ 真心请大牛出来指教.   |
回复此楼» 收录本帖的淘帖专辑推荐
 Siesta（Transiesta、ATK） | 计算-siesta |
» 猜你喜欢
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有13人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 讨论：在ATK和TranSIESTA输运中，大家对k-sampling的选取问题
已经有16人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
9楼2013-06-11 09:28:33
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
sunyang1988: 金币+2, 谢谢交流 2013-06-09 12:42:25
zhangguangping: 金币+10, ★有帮助 2013-06-09 17:26:44
感谢参与，应助指数 +1
sunyang1988: 金币+2, 谢谢交流 2013-06-09 12:42:25
zhangguangping: 金币+10, ★有帮助 2013-06-09 17:26:44
|
对于输运图，计算的时候应该是单点计算，也就是跟K点没有关系，所以一般我计算输运的时候都是取单个的K点，，具体步骤是先用很多K点优化，然后进行单点的DFT计算，最后进行输运计算，关于为什么取单个K点，我问过我导师，他说跟K没关系，不过好像Datta的第一本关于输运的里面涉及了 ，最后你这个输运图好奇怪,不知道你算得是什么的输运？ |
2楼2013-06-09 12:35:14
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
3楼2013-06-09 16:15:30
4楼2013-06-09 17:00:21
20