| 查看: 2946 | 回复: 24 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
小木虫于金虫 (正式写手)
|
[求助]
siesta-3.1 计算石墨烯条带不收敛，不知是不是我的输入有问题。
|
|||
|
大家好，我最近用siesta计算一个zigzag型石墨烯条带，宽度为16条zigzag链（16-zgnr）。可是自洽不收敛啊，我的输入结构使用vasp优化后的结构。想用siesta算电子结构。 一下是我的fdf文件： SystemName zigzag graphene nanoribbon-16 lains SystemLabel ZGNR16 NumberOfAtoms 34 NumberOfSpecies 2 %block ChemicalSpeciesLabel 1 6 C 2 1 H %endblock ChemicalSpeciesLabel PAO.BasisSize DZP LatticeConstant 1.0 Ang %block LatticeVectors 50.000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 2.4676929835305224 0.0000000000000000 0.0000000000000000 0.0000000000000000 10.0000000000000000 %endblock LatticeVectors AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies 0.4005863261073461 0.7635675166241384 0.5000000000000000 1 0.4861138744030586 0.7631620370088510 0.5000000000000000 1 0.5715306303683311 0.7652203390752147 0.5000000000000000 1 0.6569497014874649 0.7649995785970080 0.5000000000000000 1 0.7423983813159606 0.7636392092441970 0.5000000000000000 1 0.8279538544680551 0.7645027036095580 0.5000000000000000 1 0.9137570545789509 0.7656507928653636 0.5000000000000000 1 0.3151078859472453 0.7656007593586642 0.5000000000000000 1 0.9273615966545191 0.2658042298488465 0.5000000000000000 1 0.3292355303324607 0.2653823675374359 0.5000000000000000 1 0.4148299038026551 0.2630882289567006 0.5000000000000000 1 0.5003518177665096 0.2635412543771478 0.5000000000000000 1 0.5857451241441004 0.2652320086258215 0.5000000000000000 1 0.6711672235607071 0.2645919901159354 0.5000000000000000 1 0.7566318781910900 0.2636219003256479 0.5000000000000000 1 0.8421527606483915 0.2648575848498046 0.5000000000000000 1 0.8707644696675394 0.2653823675374359 0.5000000000000000 1 0.2726384033454806 0.2658042298488465 0.5000000000000000 1 0.3578472393516086 0.2648575848498046 0.5000000000000000 1 0.4433681218089099 0.2636219003256479 0.5000000000000000 1 0.5288327764392927 0.2645919901159354 0.5000000000000000 1 0.6142548758558992 0.2652320086258215 0.5000000000000000 1 0.6996481822334901 0.2635412543771478 0.5000000000000000 1 0.7851700961973451 0.2630882289567006 0.5000000000000000 1 0.5430502985125357 0.7650005785970054 0.5000000000000000 1 0.7994136738926537 0.7635665166241409 0.5000000000000000 1 0.2862429454210487 0.7656517928653611 0.5000000000000000 1 0.3720461455319449 0.7645037036095554 0.5000000000000000 1 0.4576016186840393 0.7636402092441945 0.5000000000000000 1 0.8848921140527546 0.7655997593586668 0.5000000000000000 1 0.6284693696316690 0.7652193390752172 0.5000000000000000 1 0.7138861255969410 0.7631610370088535 0.5000000000000000 1 0.2507856796782075 0.2675344989796367 0.5000000000000000 2 0.9492143203217924 0.2675344989796367 0.5000000000000000 2 %endblock AtomicCoordinatesAndAtomicSpecies %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 36 0 0.0 0 0 1 0.0 %endblock Kgrid_Monkhorst_Pack SolutionMethod diagon ElectronicTemperature 50 meV xc.functional GGA xc.authors PBE #SpinPolarized true MeshCutoff 400. Ry MaxSCFIterations 200 DM.MixingWeight 0.5 DM.NumberPulay 1 DM.Tolerance 1.d-4 #SpinPolarized .true. #%DM.InitSpin #9 +2. #18 -2. #%endblock DM.InitSpin WriteMullikenPop 1 哪位高手帮我看一下，不胜感激！ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | good resources |
» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
国自然面上复盘~欢迎讨论
已经有15人回复
-
中青基了要发朋友圈吗？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 计算石墨烯不收敛的问题
已经有10人回复
redskywei
木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 2935.1
- 散金: 57
- 红花: 4
- 帖子: 1413
- 在线: 300小时
- 虫号: 222757
- 注册: 2006-03-21
- 专业: 凝聚态物性 II ：电子结构
6楼2012-09-17 18:57:39
redskywei
木虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 2935.1
- 散金: 57
- 红花: 4
- 帖子: 1413
- 在线: 300小时
- 虫号: 222757
- 注册: 2006-03-21
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
小木虫于: 金币+5, ★有帮助 2012-09-17 15:15:47
liliangfang: 金币+1, 谢谢交流 2012-09-30 18:34:01
感谢参与，应助指数 +1
小木虫于: 金币+5, ★有帮助 2012-09-17 15:15:47
liliangfang: 金币+1, 谢谢交流 2012-09-30 18:34:01
|
DM.NumberPulay 1 设置的低了点吧。 %block kgrid_Monkhorst_Pack 1 0 0 0.0 ####这个感觉也有点小。 0 36 0 0.0 ####这个感觉也有点大。 0 0 1 0.0 %endblock Kgrid_Monkhorst_Pack MeshCutoff 400. Ry ####感觉这个有点太大了。 还有的你晶格数据有点奇怪，真空层是否小了点，特别是y方向上。 |
2楼2012-09-17 15:04:09
小木虫于
金虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1055.4
- 散金: 160
- 红花: 1
- 帖子: 419
- 在线: 293.1小时
- 虫号: 1027530
- 注册: 2010-05-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2012-09-17 15:11:58
小木虫于
金虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1055.4
- 散金: 160
- 红花: 1
- 帖子: 419
- 在线: 293.1小时
- 虫号: 1027530
- 注册: 2010-05-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2012-09-17 18:33:46
