24小时热门版块排行榜    

查看: 2086  |  回复: 5

xmc212

木虫 (正式写手)

[交流] 【求助】transiesta和tbtrans中的电压有什么不同 已有1人参与

siesta程序包中的transiesta和tbtrans都有电压选项:TS.Voltage

transiesta运行的目的是产生散射区域的.TSHS文件,用以tbtrans计算

以下是计算的电流值:
(1)transiesta TS.Voltage 0.0 eV
tbtrans TS.Voltage 0.0 eV, I= 0.0 A
tbtrans TS.Voltage 0.1 eV, I= 4.57078560E-14 A
tbtrans TS.Voltage 1.0 eV, I= 1.67146015E-09 A
(2)transiesta TS.Voltage 0.1 eV
tbtrans TS.Voltage 0.0 eV, I= 0.0 A
tbtrans TS.Voltage 0.1 eV, I= 4.56622327E-14 A
(3)transiesta TS.Voltage 1.0 eV
tbtrans TS.Voltage 0.0 eV, I= 0.0 A
tbtrans TS.Voltage 1.0 eV, I= 1.61023936E-09 A

无论transiesta加不加电压,用tbtrans计算出来的电流相差都很小,0.1 eV的时候小于0.1%,1.0 eV的时候约4%

tbtrans是实际加电压算电流的,transiesta加电压的所用是不是得到在加电压情况下的Hamiltonian?

另外,实际计算显示在误差不大的情况下,transiesta是不是就不用加电压了,算一个transiesta然后算多个tbtrans?

谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vv_china

铜虫 (小有名气)


cenwanglai(金币+1):鼓励新虫发贴。木虫期待你的精彩表现! 2010-08-19 13:25:07
TS.Voltage 是给TBTrans使用的参数,它的值应该和transiesta的计算没有关系。楼主可以仔细读读软件手册。
2楼2010-08-19 12:12:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmc212

木虫 (正式写手)

谢谢你的回答。从结果看来,电压对transiesta计算的影响确实很小,但不是0,不知道是什么原因?

(1)说明书上TS.Voltage是transiesta的选项
8.7 TranSiesta Options
8.7.1 General options
SolutionMethod
TS.SaveHS
TS.Voltage
(2)输出文件显示TS.Voltage也是transiesta的选项
ts_read_options: TranSIESTA Voltage           =    0.1000 Volts

                     ************************  
                     *   TRANSIESTA BEGIN   *  
                     ************************  

InitMesh: MESH =    64 x   128 x    80 =      655360
InitMesh: Mesh cutoff (required, used) =   100.000   108.298 Ry

TSvoltage: Bias  0.100V
引用回帖:
Originally posted by vv_china at 2010-08-19 12:12:13:
TS.Voltage 是给TBTrans使用的参数,它的值应该和transiesta的计算没有关系。楼主可以仔细读读软件手册。

3楼2010-08-19 15:36:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

★ ★
cenwanglai(金币+2):谢谢回复! 2010-08-19 17:16:24
xmc212(金币+10): 2010-08-19 19:51:00
引用回帖:
Originally posted by xmc212 at 2010-08-18 05:17:21:
siesta程序包中的transiesta和tbtrans都有电压选项:TS.Voltage

transiesta运行的目的是产生散射区域的.TSHS文件,用以tbtrans计算

以下是计算的电流值:
(1)transiesta TS.Voltage 0.0 eV
tbtrans TS.Vo ...

不知道你计算的是一个什么体系,但是可以肯定地说,transiesta的偏压和tbtrans的偏压要一致。在transiesta计算中是用到TS.Voltage的,至于用来做什么,我不好说,应为我不是很清楚,但是并不象楼上说的没有关系。起码HS文件是与偏压有关的。并且左右的GF文件也是偏压相关的。只能说,你测试的那个体系加上偏压后,对transiesta的计算结果没有什么大的影响。如果真的是没有关系的话,应该出来的结果完全一致才是对的。所以还是有关系的。这个结果是你的体系造成的。
弘德明志博学笃行
4楼2010-08-19 16:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vv_china

铜虫 (小有名气)

xmc212(金币+10): 2010-08-19 19:51:30
楼上比我对transiesta更熟悉。我只是对你的结果的猜测,也许偏压(对体系费米能级)的影响到了TBTrans才会正式考虑。因为你的模型加偏压后的电流极小,建议你换个导电性好的模型再做测试。

我在手册上找到一些相关的描述,你可以参考:
TS.TBT.HSFile (string): Scattering region .TSHS file.
Default value: SystemLabel.TSHS
respectively for the left and right electrodes and the scattering region .TSHS file.
The energy scale in tbtrans is shifted so that the Fermi level of the system, if no voltage were
applied, is zero. When computing the transmission function of a zero bias calculation, the
transmission at the Fermi level is then given by T(E=0). When there is a finite bias, the Fermi
energy of the left electrode is placed at V/2, and that of the right electrode at -V/2.
The voltage is specified by TS.Voltage.

To summarize, here we give a list of the parameters read by tbtrans from the input file (the
fdf flags):
• TS.Voltage
• kgrid Monkhorst Pack (block)
• TS.HSFileLeft。。。。
5楼2010-08-19 18:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖
引用回帖:
Originally posted by xmc212 at 2010-08-18 05:17:21:
siesta程序包中的transiesta和tbtrans都有电压选项:TS.Voltage

transiesta运行的目的是产生散射区域的.TSHS文件,用以tbtrans计算

以下是计算的电流值:
(1)transiesta TS.Voltage 0.0 eV
tbtrans TS.Vo ...

TS.Voltage这个选项对TranSIESATA和SIESTA的计算中都是有用的,具体如下:
TranSIESTA:用来提供电势自洽的边界条件,并且加上偏压后,左右电极的自能和表面格林函数都要平移±V/2,等自洽完成后,所有的输出量都是在偏压下收敛的结果,也就是说,包含了偏压的影响了。
TBtrans:TBtrans计算的时候直接读取TranSIESTA计算生成的电极和散射区的TSHS文件,这个时候虽然包含了偏压的影响。但是在计算电流的时候,公式

中的左右电极的化学势分别是电极的Fermi energy±V/2所以这个地方需要用到V,仅此而已。

并且需要补充一点的是:在TranSIESTA计算有一个温度的参数,这个参数只在TranSIESTA的自洽中用到,起到一个加速收敛的目的。而在TBtrans中,也就是上面的公式中,Fermi-Dirac分布函数中的T=0K。
弘德明志博学笃行
6楼2011-05-16 15:22:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xmc212 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天务委会开完了,明天出结果吗 +13 angus9576 2026-08-25 13/650 2026-08-25 18:24 by miguan3344
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 有没有大神帮我看看基金代码 28+3 1234567wang 2026-08-24 9/450 2026-08-25 16:51 by 北京莱茵编辑
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 +6 ziyangfang 2026-08-25 8/400 2026-08-25 16:28 by 宝贝的猪
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 什么时候开奖? +12 CrisMessi 2026-08-18 13/650 2026-08-25 14:40 by ZJTJZ
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见