24小时热门版块排行榜    

查看: 2749  |  回复: 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的回帖

redskywei

木虫 (著名写手)

收敛速度是否和体系的结构有管呢,比如说体系本来就是个亚稳态的结构,是否收敛会很慢~~~~
坠入无边的网~~
5楼2010-06-09 23:09:32
已阅   回复此楼   关注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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见