24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2147  |  回复: 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的回帖

zhangguangping

木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
mazuju028(金币+2):谢谢交流 2010-12-05 12:46:30
引用回帖:
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 ...

假如你设置的最大收敛步数是300,如果在生成*DM的时候就用了100步,那么剩下的transiesta部分就只有200步的容量了。只要这个部分超过200步就不再自洽了。为了避免这个问题,可以将最大收敛步数设置的大一点。至于你说的收敛慢的问题,你可以通过调节参数加快收敛。比如减小mixing的大小,到0.01甚至更小0.001。另外升高温度也是一个办法。还有就是你可以减小电压的间隔等等。
弘德明志博学笃行
11楼2010-12-05 09:03:16
已阅   回复此楼   关注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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见