| 查看: 2340 | 回复: 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 那么到底他判断收敛是什么标准呢~~~ 另外，俺们的收敛特别慢，除了跟参数有关外，还跟什么东西有关？ 这里面的这些量都是些啥？ |
回复此楼» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有13人回复
-
285求调剂
已经有8人回复
-
22408 266求调剂
已经有6人回复
-
264求调剂
已经有8人回复
-
材料工程专业日语生求调剂
已经有9人回复
-
297求调剂
已经有15人回复
-
277、学硕，求调剂 数一104，
已经有12人回复
-
机械专硕273请求调剂
已经有3人回复
-
080100力学316求调剂
已经有3人回复
-
288环境专硕,求调材料方向
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 梳理心情，你有啥高招？！出来晒晒吧，O(∩_∩)O哈哈~
已经有92人回复
- 求助Siesta计算熟手，想学习此软件
已经有8人回复
- 用Dmol优化收敛的结构，用vasp不收敛？
已经有9人回复
- ppt模板常用框图模块【转载】
已经有69人回复
- Lyapunov函数的选取和控制系统收敛速度快慢是怎么一种关系?
已经有10人回复
- 计算表面时：Vacuum真空层和slab层厚度收敛性测试，保持超晶胞体积不变吗？？
已经有10人回复
- siesta 表面不收敛
已经有3人回复
- 【求助】MS中dmol3模块——两种分子的共吸附优化
已经有9人回复
- 【求助】siesta中transiesta模块的收敛问题
已经有13人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
2楼2010-06-09 14:20:39
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★
zzy870720z(金币+1):谢谢补充，呵呵 2010-06-09 15:14:31
redskywei(金币+5):多谢交流，让我对参数理解更清楚了~~ 2010-06-09 23:08:06
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
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
4楼2010-06-09 16:02:51
redskywei
木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 2935.1
- 散金: 57
- 红花: 4
- 帖子: 1413
- 在线: 300小时
- 虫号: 222757
- 注册: 2006-03-21
- 专业: 凝聚态物性 II ：电子结构
5楼2010-06-09 23:09:32
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
6楼2010-06-10 09:36:40
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
★
redskywei(金币+4):多谢交流，对TS.ComplexContour.Emin我倒是没注意过～～ 2010-06-17 13:33:11
zzy870720z(金币+1):谢谢交流 2010-06-20 13:22:15
redskywei(金币+4):多谢交流，对TS.ComplexContour.Emin我倒是没注意过～～ 2010-06-17 13:33:11
zzy870720z(金币+1):谢谢交流 2010-06-20 13:22:15
|
改变DM.MixingWeight参数，对于siesta计算收敛很好，值越小收敛的速度越慢。对于金属性的一般值要小，特别是能隙比较小的体系。但对于transiesta我也不是很清楚要调什么参数。现在也在摸索阶段。如果在有电压下，我想TS.BiasContour.NumPoints这个参数很重要。可以试试这个，并且可以把TS.ComplexContour.Emin降点。虽然你设的值已经是低于最小值了，有时在有电压下，情况会有变化，再降一点有时候有作用。 |
7楼2010-06-12 17:36:44
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
8楼2010-06-19 21:49:32
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
9楼2010-06-20 21:00:18
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
10楼2010-12-04 11:17:02
