| 查看: 2332 | 回复: 13 | ||||
redskywei木虫 (著名写手)
|
[交流]
【求助】siesta中transiesta模块的收敛问题 已有3人参与
|
|
最近在做transieta的计算，发现收敛困难，看maillist里主要考虑增加屏蔽层和设置TS.UpdateDMCROnly，我在设置TS.UpdateDMCROnly时，一直不明白这个东西到底是作用在整个散射区还是仅仅与电极结构相重复的散射区中左右单元的过度层？这个coupling terms该怎么理解？ 还有期中有个人提到TS.UpdateDMCROnly + extra elec layers，又该怎么理解，是说对外电极设置这个参数么？ 我用3层Al电极做左右屏蔽层，TS.UpdateDMCROnly设置为 T了，还是收敛非差差。近300步才收敛。也不知道哪里问题哦。我的输入参数如下，请大家指教： NumberOfSpecies 3 NumberOfAtoms 287 %block ChemicalSpeciesLabel 1 6 C 2 7 N 3 13 Al %endblock ChemicalSpeciesLabel PAO.EnergyShift 200 meV PAO.BasisType split PAO.BasisType split PAO.BasisSize DZP LatticeConstant 1.0 Ang %block LatticeVectors 16.19832 0.00000000 0.00000000 0.00000000 16.19832 0.00000000 0.00000000 0.00000000 26.80535368 %endblock LatticeVectors #kgrid_cutoff 20. Ang %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 1 0.5 %endblock kgrid_Monkhorst_Pack # SPIN options xc.functional GGA # Exchange-correlation functional xc.authors PBE # Exchange-correlation version SpinPolarized F # Logical parameters are: yes or no FixSpin F TotalSpin 0.0 NonCollinearSpin F # 'T', 'F' MeshCutoff 150. Ry # Equivalent plane wave cutoff for the grid # SCF options DM.Tolerance 1.d-4 # Tolerance in maximum difference between input and output DM MaxSCFIterations 300 # Maximum number of SCF iter DM.UseSaveDM F # to use continuation files DM.MixingWeight 0.03 # New DM amount for next SCF cycle DM.NumberPulay 5 DM.MixSCF1 F DM.PulayOnFile F # Store in memory ('F') or in files ('T') # NeglNonOverlapInt T # 'F'=do not neglect SolutionMethod transiesta # OrderN or Diagon or transiesta ElectronicTemperature 300 K # Temp. for Fermi smearing # MD options MD.TypeOfRun CG # Type of dynamics: MD.VariableCell F MD.NumCGsteps 000 # Number of CG steps for coordinate optimization ZM.UnitsLength Ang #the units of length used during Z-matrix input ZM.UnitsAngle deg #the units of angles used during Z-matrix input ZM.ForceTolLength 0.02 eV/Ang #controls the convergence with respect to forces on Z-matrix lengths ZM.ForceTolAngle 0.00356549 Ry/rad #controls the convergence with respect to forces on Z-matrix angles ZM.MaxDisplLength 0.1 Ang # controls the maximum change in a Zmatrix length during an optimisation step ZM.MaxDisplAngle 0.003 rad # controls the maximum change in a Z-matrix angle during an optimisation step ZM.CalcAllForces T # Default value [ Last edited by redskywei on 2010-5-29 at 18:21 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 SIESTA |
» 猜你喜欢
315求调剂
已经有15人回复
-
070300化学学硕311分求调剂
已经有9人回复
-
一志愿南昌大学，085600，344分求调剂
已经有7人回复
-
一志愿 江南大学 085602 化工专硕 338分求调剂
已经有13人回复
-
化学0703-一志愿211-338分求调剂
已经有8人回复
-
0703调剂，一志愿天津大学319分
已经有8人回复
-
085602调剂 初试总分335
已经有9人回复
-
336材料与化工085600求调剂
已经有8人回复
-
085600材料与化工301分求调剂院校
已经有6人回复
-
071000生物学调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手请进，帮忙看下我翻译的一段话（中译英），有没有问题 谢谢
已经有2人回复
- MS 各个模块的功能
已经有3人回复
- siesta-3.1编译成功，但有几个问题求助
已经有7人回复
- siesta 计算电子占据数的问题
已经有3人回复
- 用siesta画声子谱的问题
已经有12人回复
- 利用siesta做verlet MD时出现问题siesta: WARNING: Atoms
已经有3人回复
- siesta 表面不收敛
已经有3人回复
- 【求助】siesta多个/种分子优化的问题
已经有7人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【求助】关于siesta的几何结构优化中断后继续的问题【完结】
已经有6人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
2楼2010-05-29 18:23:30
163.com
金虫 (著名写手)
实习版主
- 应助: 2 (幼儿园)
- 贵宾: 0.016
- 金币: 1489.8
- 散金: 1347
- 红花: 32
- 帖子: 1215
- 在线: 87小时
- 虫号: 923817
- 注册: 2009-12-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★
redskywei(金币+2):谢谢建议，我感觉这个参数应该不是很多影响的，不过希望可以把这个参数的意义说的更清楚些～ 2010-05-29 21:22:18
redskywei(金币+3):我想同样的第二周期的原子，应该没这么大的影响吧，我把原子替换看看～～ 2010-05-29 21:24:34
lvjian8596(金币+1):感谢交流！鼓励发有价值的话题！ 2010-06-22 22:23:49
redskywei(金币+2):谢谢建议，我感觉这个参数应该不是很多影响的，不过希望可以把这个参数的意义说的更清楚些～ 2010-05-29 21:22:18
redskywei(金币+3):我想同样的第二周期的原子，应该没这么大的影响吧，我把原子替换看看～～ 2010-05-29 21:24:34
lvjian8596(金币+1):感谢交流！鼓励发有价值的话题！ 2010-06-22 22:23:49
|
这个收敛参数没有问题，200meV。 MeshCutoff=150.Ry 你在看一下其他参数。 还有可能就是体系的问题，同样结构的BN比Graphene就会要难收敛好多，体系完全一样，把原子替换成BN之后 计算量是Graphene的好几倍 |
3楼2010-05-29 18:45:14
redskywei
木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 2935.1
- 散金: 57
- 红花: 4
- 帖子: 1413
- 在线: 300小时
- 虫号: 222757
- 注册: 2006-03-21
- 专业: 凝聚态物性 II ：电子结构
4楼2010-06-08 21:27:52
redskywei
木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 2935.1
- 散金: 57
- 红花: 4
- 帖子: 1413
- 在线: 300小时
- 虫号: 222757
- 注册: 2006-03-21
- 专业: 凝聚态物性 II ：电子结构
5楼2010-06-08 21:30:32
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★ ★
redskywei(金币+3):感谢交流，我用al做屏蔽层的，3层不知道够不够，呵呵～～ 2010-06-20 10:18:20
qasd(金币+2):谢谢 2010-06-20 11:32:58
redskywei(金币+3):感谢交流，我用al做屏蔽层的，3层不知道够不够，呵呵～～ 2010-06-20 10:18:20
qasd(金币+2):谢谢 2010-06-20 11:32:58
|
我看到mailist上面确实是有些人讨论过碳纳米管的收敛性问题。我没有涉及这个，所以我基本上没有遇到收敛性的问题。当我遇到收敛问题的时候，一般通过增加屏蔽层的层数就可以解决问题！或者是通过TS.UpdateDMCROnly这个选项设置为F就可以了。要是这个选项设置为F解决问题的话，一般说明是屏蔽曾数目太少了。我个人认为你的屏蔽层是不是有点少啊？C原子属于轻元素，所以相应的屏蔽效应比较差，所以你的屏蔽曾要足够多。我用的都是重金属原子，所以屏蔽层一般在4~5层就解决了！ |
6楼2010-06-19 21:43:23
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
7楼2010-06-20 11:04:12
redskywei
木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 2935.1
- 散金: 57
- 红花: 4
- 帖子: 1413
- 在线: 300小时
- 虫号: 222757
- 注册: 2006-03-21
- 专业: 凝聚态物性 II ：电子结构
8楼2010-06-21 12:05:39
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
9楼2010-06-21 15:34:39
redskywei
木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 2935.1
- 散金: 57
- 红花: 4
- 帖子: 1413
- 在线: 300小时
- 虫号: 222757
- 注册: 2006-03-21
- 专业: 凝聚态物性 II ：电子结构
10楼2010-06-22 21:24:58
