| 查看: 1547 | 回复: 8 | ||||
| 当前主题已经存档。 | ||||
zhangguangping木虫 (著名写手)
|
[交流]
【求助】siesta计算中关于K点的设置【完结】 已有5人参与
|
|||
|
我想和大家讨论一下siesta以及凡是用siesta做基础的其他软件象如smeagol，transiesta在计算中k点的选取问题！我这个地方一直不是很明白！ 在siesta的计算中，都要构造一个原胞，用周期结构计算的方法来计算体系！其中体系的周期用下面的block设置。 %block LatticeVectors 11.53504 0.00000000 0.00000000 0.00000000 9.98964 0.00000000 0.00000000 0.00000000 30.00000000 %endblock LatticeVectors 但是要想利用siesta计算一个孤立体系，这时候就可以把这个原胞构造的大一点（足够大），使得该体系按照设置的周期平移后没有相互作用！这样计算的结果就是孤立体系的结果！ 在计算周期体系的时候要进行K点的设置，可以用如下的block设置， %block kgrid_Monkhorst_Pack 10 0 0 0.0 0 10 0 0.0 0 0 10 0.0 %endblock kgrid_Monkhorst_Pack 如果计算的体系确确实实是周期性的，例如晶体，这时候这里面每个方向上的K点就要选取的足够多，使得体系计算结果收敛，尤其是计算金属的时候，这个K点可能设置的比较大！ 问题①：我上面的理解对吗？ 问题②：如果我的体系在某一个方向上,如x方向有周期结构，但其他的方向没有，这个时候K点设置里是不是第一行设置的大一些，使得结果收敛，第二行，第三行直接设置1就可以！也就是说没有周期方向对应的K点行上取1就足够！（假设原胞平移矢量是正交的） 问题③：如果问题②是对的，那么上面说到的如果计算一个孤立体系，我把原胞设置的足够大没之后，在设置k点的时候各个行都设置1即可？也就是说只计算自己虚构的周期体系的Gamma点就行了！ 问题④：同一个体系的计算量是不是与k点的多少成正比，也就是说只计算1个K点花10分钟，要是将K点增加到10个，那么是不是计算时间为100分钟？ 附加一个问题，在进行几何结构优化的时候怎么感觉SolutonMethod Diagon TypeOfRun CG 这两个组合比较慢呢？有什么方法能使几何结构优化加快啊！一个简单的分子，G09优化用了5分钟都不到，siesta怎么得花很长时间呢？ 在这儿先谢谢参与该讨论的虫子们！ 谢谢高手指教！ [ Last edited by zhangguangping on 2010-7-20 at 14:30 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta |
» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2010-04-16 12:11:40
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gavinliu7390(金币+1):谢谢交流！ 2010-04-16 13:20
小木虫(金币+0.5):给个红包，谢谢回帖交流
gavinliu7390(金币+1):谢谢交流！ 2010-04-16 13:20
|
第一个问题，我也是这么理解的。 第二个问题，siesta是可以自己控制有几个方向的周期性吗？不知道耶？我每次都是三个方向都做。 别的就不懂了。 我是新手，问你一个小问题可以吗？ BandLines (data block): Speci es the lines along which band energies are calculated (usually along high-symmetry directions). An example for an FCC lattice is: %block BandLines 1 1.000 1.000 1.000 L # Begin at L 20 0.000 0.000 0.000 \Gamma # 20 points from L to gamma 25 2.000 0.000 0.000 X # 25 points from gamma to X 30 2.000 2.000 2.000 \Gamma # 30 points from X to gamma %endblock BandLines 以上出自siesta2.0.2 manual的33页，上面的L，X，Gamma是什么意思？我有一次在算的时候，出现了warning：L应该up to 3.不大懂怎么搞呀？或者你有什么有关这方面的书，推荐我看看呗。（PS：我固体物理还没学到能带那儿。） |
3楼2010-04-16 12:22:29
lz830617
铁杆木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 4856.1
- 散金: 168
- 红花: 3
- 帖子: 312
- 在线: 352.2小时
- 虫号: 584129
- 注册: 2008-08-01
- 性别: GG
- 专业: 光子与光电子器件
4楼2010-04-16 12:53:23
lz830617
铁杆木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 4856.1
- 散金: 168
- 红花: 3
- 帖子: 312
- 在线: 352.2小时
- 虫号: 584129
- 注册: 2008-08-01
- 性别: GG
- 专业: 光子与光电子器件
5楼2010-04-16 13:01:04
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
6楼2010-04-16 20:36:25
7楼2010-04-16 21:27:04
8楼2010-04-19 00:29:44
9楼2010-04-19 10:32:05