| 查看: 2164 | 回复: 10 | |||
[交流]
【求助】K-S 方程的迭代解法
|
麻烦您详细解释一下图片中的解法。我刚接触vasp，还在看书阶段，还没有操作过。这张图是vasp官方的指导幻灯片里的。vasp手册里有一些说明，我看不懂 。这里的能量是怎么求的？我怎么看来看去，残余矢量是0？？？ 我现在金币较少，如果能帮我弄明白的话，金币还会补加的。具体的细节可能会很复杂，我现在只求能理解整个解的过程。  非常感谢 原文件：http://cms.mpi.univie.ac.at/vasp-workshop/slides/comput_mat.pdf [ Last edited by ss_kampf on 2011-3-2 at 14:14 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
27博士申请（市政/环境）
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有个求解微分方程用到高斯赛德尔迭代的，算不太明白啊，各位帮忙看看
已经有5人回复
- 指数方程的解法
已经有4人回复
- 都不知道怎么命名标题：迭代法解方程组？
已经有16人回复
- 【求助】求助一个二元二次回归方程的解法和答案
已经有11人回复
- 【求助完毕】为何要计算空带？
已经有5人回复
- 【原创】第一原理计算的一些心得
已经有160人回复
- 【求助】一个二阶常微分方程的解法（重奖金币！）
已经有16人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】
已经有15人回复
- 【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
- 【重点讨论】Castep中的Empty Band的用处和设置？（参与讨论就有机会赢取大礼包！）
已经有54人回复
★ ★ ★
ss_kampf(金币+25): 谢谢版主应助！ 2011-03-02 18:02:47
youzhizhe(金币+3): 谢谢指导交流。 2011-03-03 00:24:52
ss_kampf(金币+25): 谢谢版主应助！ 2011-03-02 18:02:47
youzhizhe(金币+3): 谢谢指导交流。 2011-03-03 00:24:52
|
开始的trial vectors应该不是H的本征矢量 如果是的话就没有必要修正trail vector，因为trail vector对应的能量就是所要求的能量 所以所说的residual vector不会为零，为零的时候表示已经收敛了 修正之后通过对比前后的电荷密度，当密度之间差距小于指定值时就收敛了 这就是K-S方程的迭代求解吧 另外trail vector一般都是用平面波基组或是原子轨道基组做展开 |
2楼2011-03-02 16:46:26
3楼2011-03-02 18:03:21
4楼2011-03-02 18:51:18
5楼2011-03-02 18:55:33
6楼2011-03-02 19:36:42
|
彻底糊涂了。vasp手册7.1节里面讲了好多方法，看不懂…… 111页的图弄上来了，是不是说： 第一步先从尝试波函数构建Hamitonian，然后对角化得出新的波函数，和能量E。 第二步弄迭代残余矢量的时候，是不是哈密顿要用第一步得出新的波函数来重新构建，E是上一步得出来的，这样就会产生残余矢量。进而得出有一个新的波函数。 第三步用上步得出的波函数，再构建哈密顿，再对角化，得出更新的波函数和能量。 不知道这样对不对？ 我老师怎么说VASP里面不用对角化的方法，因为计算量大。只有迭代的方法。   [ Last edited by ss_kampf on 2011-3-2 at 22:06 ] |
7楼2011-03-02 22:02:35
8楼2011-03-02 22:22:55
9楼2011-03-02 22:29:30
10楼2011-03-02 23:10:18
11楼2011-03-03 09:00:19
