24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2353  |  回复: 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

木虫 (著名写手)

★ ★
hedaors(金币+1):Thanks 2010-06-21 07:51:37
zzy870720z(金币+1):谢谢提示 2010-06-21 08:12:17
引用回帖:
Originally posted by zhangguangping at 2010-06-19 21:49:32:

你是怎么看出TS.ComplexContour.Emin值的设置已经低于最小值了?

首先估计一下,一般的-40.0eV就可以了。
然后就是用siesta单独算一下。
9楼2010-06-20 21:00:18
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 267求调剂 +5 再忙也要吃饭啊 2026-04-09 5/250 2026-04-09 18:47 by stone_128
[考研] 085600材料与化工,求调剂 +11 won_qii 2026-04-07 11/550 2026-04-09 17:03 by luoyongfeng
[考研] 314求调剂 +20 wakeluofu 2026-04-09 21/1050 2026-04-09 16:14 by zhuimr
[考研] 085404 293求调剂 +7 勇远库爱314 2026-04-08 7/350 2026-04-09 16:02 by 猪会飞
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 287求调剂 +10 Fnhc 2026-04-07 16/800 2026-04-08 10:07 by xingguangj
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
[考研] 0854电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-05 6/300 2026-04-07 22:16 by hemengdong
[考研] 一志愿北京化工085600 310分求调剂 +20 0856材料与化工3 2026-04-04 22/1100 2026-04-07 15:14 by 上岸快快
[考研] 0854求调剂 +4 assdll 2026-04-05 4/200 2026-04-06 12:29 by 中飞院空管学院
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 328分调剂 +6 门men 2026-04-04 6/300 2026-04-05 13:40 by imissbao
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 387求调剂 +4 爱吃片豆土 2026-04-03 5/250 2026-04-04 08:10 by 岸上的一条鱼
信息提示
请填处理意见