| 查看: 2324 | 回复: 7 | ||
xwu10铜虫 (小有名气)
|
[求助]
SIESTA 计算石墨烯，模型增大体系不收敛 已有1人参与
|
|
各位好 我在用SIESTA进行石墨烯的简单计算时，遇到了一个问题。我用MS进行石墨烯体系的建模，当石墨烯体系较小时（小于30个原子），SIESTA计算可以很好地收敛；然而当我增大石墨烯体系时（其他参数不变），去发现体系始终无法收敛，这个时候不管我修改何种参数，体系都无法收敛。我想问的是，SIESTA针对大体系计算难道就会出现收敛性问题吗？这是软件本身的bug还是什么其他原因？各位是否有可靠建议？（我试过修改MixingWeight和NumberPulay参数，以及Kgrid，结果依旧不能收敛） 两种体系的fdf文件如下，也供初学者学习参考 #graphene-small supercell SystemName graphene SystemLabel graphene NumberOfAtoms 30 NumberOfSpecies 1 %block ChemicalSpeciesLabel 1 6 C %endblock ChemicalSpeciesLabel %block PAO.Basis # Define Basis set C 2 # Species label, number of l-shells n=2 0 1 # n, l, Nzeta 4.088 1.000 n=2 1 1 # n, l, Nzeta, Polarization, NzetaPol 4.870 1.000 %endblock PAO.Basis XC.functional GGA LatticeConstant 1.0 Ang %block LatticeVectors 20.00000000 0.0000000000 0.0000000000 0.000000000 15.0000000000 0.0000000000 0.000000000 0.0000000000 20.0000000000 %endblock LatticeVectors AtomicCoordinatesFormat Ang %block AtomicCoordinatesAndAtomicSpecies 1.7 2.13 1.23 1 1.7 2.841 0 1 1.7 4.261 0 1 1.7 4.971 1.23 1 1.7 6.391 1.23 1 1.7 7.101 0 1 1.7 8.522 0 1 1.7 9.232 1.23 1 1.7 10.652 1.23 1 1.7 11.362 0 1 1.7 2.13 3.69 1 1.7 2.841 2.46 1 1.7 4.261 2.46 1 1.7 4.971 3.69 1 1.7 6.391 3.69 1 1.7 7.101 2.46 1 1.7 8.522 2.46 1 1.7 9.232 3.69 1 1.7 10.652 3.69 1 1.7 11.362 2.46 1 1.7 2.13 6.15 1 1.7 2.841 4.92 1 1.7 4.261 4.92 1 1.7 4.971 6.15 1 1.7 6.391 6.15 1 1.7 7.101 4.92 1 1.7 8.522 4.92 1 1.7 9.232 6.15 1 1.7 10.652 6.15 1 1.7 11.362 4.92 1 %endblock AtomicCoordinatesAndAtomicSpecies %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 10 0 0.0 0 0 10 0.5 %endblock Kgrid_Monkhorst_Pack MeshCutoff 500.0 Ry MaxSCFIterations 1000 DM.MixingWeight 0.05 DM.NumberPulay 10 DM.Tolerance 0.0001 WriteCoorXmol .true. BandLinesScale ReciprocalLatticeVectors %block BandLines 1 0.0 0.0 0.0 \Gamma 100 0.0 0.0 0.5 X %endblock BandLines SolutionMethod diagon TS.WriteHS .true. #graphene-large supercell SystemName graphene SystemLabel graphene NumberOfAtoms 50 NumberOfSpecies 1 %block ChemicalSpeciesLabel 1 6 C %endblock ChemicalSpeciesLabel %block PAO.Basis # Define Basis set C 2 # Species label, number of l-shells n=2 0 1 # n, l, Nzeta 4.088 1.000 n=2 1 1 # n, l, Nzeta, Polarization, NzetaPol 4.870 1.000 %endblock PAO.Basis XC.functional GGA XC.authors PBE LatticeConstant 1.0 Ang %block LatticeVectors 20.00000000 0.0000000000 0.0000000000 0.000000000 15.0000000000 0.0000000000 0.000000000 0.0000000000 20 %endblock LatticeVectors AtomicCoordinatesFormat Ang %block AtomicCoordinatesAndAtomicSpecies 1.7 2.13 1.23 1 1.7 2.841 0 1 1.7 4.261 0 1 1.7 4.971 1.23 1 1.7 6.391 1.23 1 1.7 7.101 0 1 1.7 8.522 0 1 1.7 9.232 1.23 1 1.7 10.652 1.23 1 1.7 11.362 0 1 1.7 2.13 3.69 1 1.7 2.841 2.46 1 1.7 4.261 2.46 1 1.7 4.971 3.69 1 1.7 6.391 3.69 1 1.7 7.101 2.46 1 1.7 8.522 2.46 1 1.7 9.232 3.69 1 1.7 10.652 3.69 1 1.7 11.362 2.46 1 1.7 2.13 6.15 1 1.7 2.841 4.92 1 1.7 4.261 4.92 1 1.7 4.971 6.15 1 1.7 6.391 6.15 1 1.7 7.101 4.92 1 1.7 8.522 4.92 1 1.7 9.232 6.15 1 1.7 10.652 6.15 1 1.7 11.362 4.92 1 1.7 2.13 8.61 1 1.7 2.841 7.38 1 1.7 4.261 7.38 1 1.7 4.971 8.61 1 1.7 6.391 8.61 1 1.7 7.101 7.38 1 1.7 8.522 7.38 1 1.7 9.232 8.61 1 1.7 10.652 8.61 1 1.7 11.362 7.38 1 1.7 2.13 11.07 1 1.7 2.841 9.84 1 1.7 4.261 9.84 1 1.7 4.971 11.07 1 1.7 6.391 11.07 1 1.7 7.101 9.84 1 1.7 8.522 9.84 1 1.7 9.232 11.07 1 1.7 10.652 11.07 1 1.7 11.362 9.84 1 %endblock AtomicCoordinatesAndAtomicSpecies %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 10 0 0.0 0 0 10 0.5 %endblock Kgrid_Monkhorst_Pack MeshCutoff 500.0 Ry MaxSCFIterations 1000 DM.MixingWeight 0.1 DM.NumberPulay 10 DM.Tolerance 0.0001 WriteCoorXmol .true. BandLinesScale ReciprocalLatticeVectors %block BandLines 1 0.0 0.0 0.0 \Gamma 100 0.0 0.0 0.5 X %endblock BandLines SolutionMethod diagon TS.WriteHS .true. |
回复此楼» 猜你喜欢
青基如何破局
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有7人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 静态计算，加U计算不收敛！当U>=7时，就不收敛了，当U小于7时很快收敛！
已经有8人回复
- 计算机程序设计艺术（第三版）1~3卷
已经有46人回复
- 石墨烯转移膜层不连续
已经有6人回复
- 请教如何加快siesta的收敛速度
已经有9人回复
- 双层石墨烯的间距多大啊 优化后 谢谢
已经有7人回复
- X-crysden安装各种问题
已经有6人回复
- VASP NEB 不收敛问题
已经有25人回复
- VASP 计算石墨烯与石墨烯带的区别
已经有11人回复
- Siesta从官网下载赝势计算出错
已经有7人回复
- VASP计算表面吸附 但是能量总是不收敛/。。。
已经有9人回复
- 求助 高斯计算过程中一直不收敛 scf(maxcyc=500,xqc)
已经有9人回复
- 算例第一个时间步就不收敛？
已经有14人回复
- 用Dmol3结构优化时不收敛
已经有5人回复
- 被积函数很快收敛到0，为什么积分不收敛。
已经有3人回复
- siesta结构优化收敛标准
已经有5人回复
- 结构优化 siesta 和 vasp
已经有7人回复
- siesta-3.1 计算石墨烯条带不收敛，不知是不是我的输入有问题。
已经有24人回复
- SIESTA计算能级不准怎么解决
已经有10人回复
- 求助关于用vasp计算二维石墨烯纳米带所得到的能带结果与报道对不上
已经有7人回复
- vasp计算石墨烯能带时，dirac点是打开的，带隙不为零，与静态的DOSCAR也对不上
已经有10人回复
- 【求助】请问用siesta算输运怎么建模啊 我是搞石墨烯
已经有7人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
- 【求助】siesta中transiesta模块的收敛问题
已经有13人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
xwu10
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 544.6
- 散金: 283
- 帖子: 240
- 在线: 273.8小时
- 虫号: 1366976
- 注册: 2011-08-12
- 专业: 金属材料表面科学与工程
2楼2015-01-26 14:55:53
锐利的碎片
木虫 (正式写手)
star watcher
- 1ST强帖: 1
- 应助: 136 (高中生)
- 金币: 3637.1
- 散金: 252
- 红花: 22
- 帖子: 988
- 在线: 1224.9小时
- 虫号: 961933
- 注册: 2010-03-05
- 专业: 凝聚态物性 II ：电子结构
3楼2015-01-26 21:09:01
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
4楼2015-01-26 21:42:12
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
5楼2015-01-26 21:45:10
guohuazhong
至尊木虫 (职业作家)
- 1ST强帖: 1
- 应助: 52 (初中生)
- 贵宾: 0.032
- 金币: 35131.5
- 散金: 5
- 红花: 19
- 沙发: 1
- 帖子: 3924
- 在线: 885.6小时
- 虫号: 550098
- 注册: 2008-04-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2015-01-27 20:44:34
xwu10
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 544.6
- 散金: 283
- 帖子: 240
- 在线: 273.8小时
- 虫号: 1366976
- 注册: 2011-08-12
- 专业: 金属材料表面科学与工程
7楼2015-01-28 09:44:28
8楼2016-05-12 22:58:07
10