| 查看: 2703 | 回复: 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 |
» 猜你喜欢
什么时候开奖？
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手请进，帮忙看下我翻译的一段话（中译英），有没有问题 谢谢
已经有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
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
9楼2010-06-21 15:34:39
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.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
