24小时热门版块排行榜    

查看: 3340  |  回复: 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怎么处理呢?

真心请大牛出来指教.

K点在TranSIESTA计算中的作用
K点在TranSIESTA计算中的作用-1
回复此楼
弘德明志博学笃行
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lu_wei07

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

引用回帖:
2楼: Originally posted by lu_wei07 at 2013-06-09 12:35:14
对于输运图,计算的时候应该是单点计算,也就是跟K点没有关系,所以一般我计算输运的时候都是取单个的K点,,具体步骤是先用很多K点优化,然后进行单点的DFT计算,最后进行输运计算,关于为什么取单个K点,我问过我 ...

谢谢你的回答。

计算输运的时候,是需要对K点取样的。我不明白你说的跟K点没有关系是什么意思,能详细指教一下吗?
弘德明志博学笃行
3楼2013-06-09 16:15:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lu_wei07

银虫 (小有名气)

引用回帖:
3楼: Originally posted by zhangguangping at 2013-06-09 16:15:30
谢谢你的回答。

计算输运的时候,是需要对K点取样的。我不明白你说的跟K点没有关系是什么意思,能详细指教一下吗?...

我用的是我导师自己开发的软件(PRB 70.085410),类似于ATK的那种,我确实问过我导师,为什么计算输运前是进行单点的DFT计算(不是优化),而不是取很多的K点,老师的回答只是说跟K点没有关系了,可是我明明发现Qup-Qdown 不一样。
如果你有兴趣读我导师的博士后,当然你可以直接问我导师的(同济大学柯三黄),我相信他在电子输运这块应该是国内顶尖的,
这只是一场游戏
4楼2013-06-09 17:00:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

内容已删除
弘德明志博学笃行
5楼2013-06-09 17:17:23
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
建议读读有关波函数匹配的方法来计算电子输运的文献,就可以比较清楚地理解transimission spectra等计算中的波矢k的含义。
比如:PHYSICAL REVIEW B 70, 045417 (2004)

有关k点取样对transimission spectra计算的影响:
M. Stilling , K. Stokbro & K. Flensberg (2007): Crystalline Magnetotunnel Junctions: Fe-mgo-fe, Fe-
feomgo-fe And Fe-aumgoau-fe, Molecular Simulation, 33:7, 557-561
http://www.tandfonline.com/doi/a ... 930508#.UbZstpzRL0c
6楼2013-06-11 08:18:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

可能是因为你垂直于传输方向上的横截面积取的比较大,中间的导体取的是分子体系吧。

但是对于像Fe/MgO/Fe这些晶体组成的磁隧道结,用单个k点就不够了。
7楼2013-06-11 08:20:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

引用回帖:
6楼: Originally posted by valenhou001 at 2013-06-11 08:18:13
建议读读有关波函数匹配的方法来计算电子输运的文献,就可以比较清楚地理解transimission spectra等计算中的波矢k的含义。
比如:PHYSICAL REVIEW B 70, 045417 (2004)

有关k点取样对transimission spectra计算 ...

谢谢。我去读读。
弘德明志博学笃行
8楼2013-06-11 09:27:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

引用回帖:
7楼: Originally posted by valenhou001 at 2013-06-11 08:20:36
可能是因为你垂直于传输方向上的横截面积取的比较大,中间的导体取的是分子体系吧。

但是对于像Fe/MgO/Fe这些晶体组成的磁隧道结,用单个k点就不够了。

我电极采用的是3*3的Au,中间是分子。
你说的应该是如果是“Fe/MgO/Fe这些晶体组成的磁隧道结,用单个k点就够了”是吧?
弘德明志博学笃行
9楼2013-06-11 09:28:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
9楼: Originally posted by zhangguangping at 2013-06-11 10:28:33
我电极采用的是3*3的Au,中间是分子。
你说的应该是如果是“Fe/MgO/Fe这些晶体组成的磁隧道结,用单个k点就够了”是吧?...

Fe/MgO/Fe这些晶体组成的磁隧道结,用单个k点不够。
10楼2013-06-11 11:02:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangguangping 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 欢迎采矿、地质、岩土、计算机、人工智能等专业的同学报考 +6 pin8023 2026-02-28 8/400 2026-03-02 17:13 by 0854蹲调剂
[考研] 303求调剂 +5 今夏不夏 2026-03-01 5/250 2026-03-02 15:01 by 向上的胖东
[考研] 282求调剂 +4 2103240126 2026-03-02 5/250 2026-03-02 13:45 by littlehu66
[考研] 材料学硕318求调剂 +9 February_Feb 2026-03-01 9/450 2026-03-02 13:31 by njzyff
[考研] 338求调剂 +3 18162027187 2026-03-02 3/150 2026-03-02 13:12 by houyaoxu
[考研] 一志愿山东大学材料与化工325求调剂 +3 半截的诗0927 2026-03-02 3/150 2026-03-02 12:58 by houyaoxu
[基金申请] 面上模板改不了页边距吧? +6 ieewxg 2026-02-25 7/350 2026-03-02 12:44 by stidwellNK
[基金申请] 成果系统访问量大,请15分钟后再尝试。由此给您造成的不便,敬请谅解。 +5 xhuama 2026-03-02 5/250 2026-03-02 12:34 by stidwellNK
[考研] 材料工程专硕283求调剂 5+4 ,!? 2026-03-02 6/300 2026-03-02 11:07 by 黑!在干嘛
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪!?!
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +6 kyf化工 2026-02-28 7/350 2026-03-02 10:56 by 无际的草原
[考研] 0856材料调剂 +4 沿岸有贝壳OUC 2026-03-02 4/200 2026-03-02 10:19 by 公瑾逍遥
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师!
[考研] 材料工程269求调剂 +3 白刺玫 2026-03-02 3/150 2026-03-02 09:25 by 一休哥FU
[考研] 材料复试调剂 +4 学材料的点 2026-03-01 5/250 2026-03-02 08:26 by houyaoxu
[考研] 0857调剂 +4 一ll半 2026-02-28 5/250 2026-03-02 02:33 by 908055542
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
信息提示
请填处理意见