我想和大家讨论一下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 ] 返回小木虫查看更多