24小时热门版块排行榜    

查看: 2293  |  回复: 20
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

redskywei

木虫 (著名写手)

[交流] 【求助】siesta中transiesta模块的收敛标准是啥~~ 已有4人参与

如题,我在用transesta计算输运时发现,最后收敛的的时候,这些值都比较大,他就停止迭代了:
TranSiesta: Qsol,Qtot:  515.152  512.000    1.000
transiesta: 290   -19912.8627   -19914.3989   -19914.3989  0.4114 -5.3182
TranSiesta: Qsol,Qtot:  514.839  512.000    1.000
transiesta: 291   -19908.1417   -19913.3302   -19913.3302  0.4050 -5.3182

那么到底他判断收敛是什么标准呢~~~
另外,俺们的收敛特别慢,除了跟参数有关外,还跟什么东西有关?
这里面的这些量都是些啥?
回复此楼
坠入无边的网~~
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ifmc1234

木虫 (著名写手)


zzy870720z(金币+1):感谢交流,谢谢 2010-06-09 16:20:27
redskywei(金币+5):感谢指点,qsol意思是解方程得到的体系电荷值哦 2010-06-09 22:52:21
引用回帖:
Originally posted by redskywei at 2010-06-08 21:23:46:
如题,我在用transesta计算输运时发现,最后收敛的的时候,这些值都比较大,他就停止迭代了:
TranSiesta: Qsol,Qtot:  515.152  512.000    1.000
transiesta: 290   -19912.8627   -19914.3989   -19914.3989 ...

Qsol这行指的是电荷,Qtot是体现原来的电荷值,前一个是计算得到的电荷值,看他们的值好像是加电压下的值。要不就是有电荷累积,不合理了。
transiesta这行-19912等等是能量值,最好一个是fermi值,倒数第二个是dDmax,它要小于DM.Tolerance才行。
4楼2010-06-09 16:02:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

zhangguangping

木虫 (著名写手)


redskywei(金币+5):多谢指点,可是siesta这么快收敛,transiesta这么这么难的收敛~~~ 2010-06-09 14:23:35
zzy870720z(金币+1):谢谢指点 2010-06-09 15:14:21
引用回帖:
Originally posted by redskywei at 2010-06-08 14:23:46:
如题,我在用transesta计算输运时发现,最后收敛的的时候,这些值都比较大,他就停止迭代了:
TranSiesta: Qsol,Qtot:  515.152  512.000    1.000
transiesta: 290   -19912.8627   -19914.3989   -19914.3989 ...

你这个都迭代291次了,是不是你设置的最大迭代次数是300,其中在transiesta开始的siesta计算中用了9个迭代。所以说你现在超出了迭代的设置。如果不超的话,它的收敛标准是DM.Tolerance所设置的,显然你现在是0.4050,估计没有达到你设置的标准。一般都设置为1.0d-4或者是1.0d-3
弘德明志博学笃行
2楼2010-06-09 14:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)


zzy870720z(金币+1):谢谢补充,呵呵 2010-06-09 15:14:31
redskywei(金币+5):多谢交流,让我对参数理解更清楚了~~ 2010-06-09 23:08:06
transiesta难收敛的情况也不好说!我也遇到一次。那次是我选取的缓冲层有点少。缓冲层不但起到屏蔽作用,还有助于加快收敛。所以你可以试试把扩展分子两边的电极部分取得多一点。或者有个选项你设置一下也可以起作用的。TS.UpdateDMCROnly  F.我自己认为如果这个选项设置为T就表示在进行Transiesta迭代计算的时候密度矩阵只是更新扩展分子中间散射区的,不更新属于电极的那些原子。要是设置这个为F,那么整个扩展分子的密度矩阵在迭代过程中都要更新。要是这个设置为F之后收敛好了,我认为缓冲区选取的有点少!所以应该加大这个缓冲区。这个紧紧是我自己的理解!一般自由分子距离电极中间要有4-5层的金属原子,这个缓冲一般情况足够了。但是为了计算的需要(减少计算时间),不是选的这么多!所以这个问题遇到的就多一些!
弘德明志博学笃行
3楼2010-06-09 14:35:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)

收敛速度是否和体系的结构有管呢,比如说体系本来就是个亚稳态的结构,是否收敛会很慢~~~~
坠入无边的网~~
5楼2010-06-09 23:09:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +7 rare12345 2026-03-18 7/350 2026-03-18 14:31 by laoshidan
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
请填处理意见