24小时热门版块排行榜    

查看: 3050  |  回复: 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-10 12:09:32
redskywei(金币+5):多谢指点,我尝试过你提到的DM.MixingWeight参数,不过效果不好,但是DM.NumberPulay我最大也就设置10 2010-06-12 15:54:07
引用回帖:
Originally posted by redskywei at 2010-06-09 23:09:32:
收敛速度是否和体系的结构有管呢,比如说体系本来就是个亚稳态的结构,是否收敛会很慢~~~~

咱们都只是相当于作单点计算,收敛快慢,我想应该和你设置的参数有关。
你设置的DM.MixingWeight 是多少?可以设置小点的比如 0.01。而DM.NumberPulay设置大些,8,10,20,40这些可以试试。一般情况20应该是足够了。我一般先设置在8左右。
6楼2010-06-10 09:36:40
已阅   回复此楼   关注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的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 3/150 2026-11-10 01:40 by FChKX8meu1ds
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见