| 查看: 1951 | 回复: 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一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有6人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
小木虫看见有人已经知道结果了
已经有19人回复
-
2026年国自然面上资助率
已经有8人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
微信指数没变化，科研之友没阅读
已经有17人回复
-
申请2027年材料/化学类博士
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有个求解微分方程用到高斯赛德尔迭代的，算不太明白啊，各位帮忙看看
已经有5人回复
- 指数方程的解法
已经有4人回复
- 都不知道怎么命名标题：迭代法解方程组？
已经有16人回复
- 【求助】求助一个二元二次回归方程的解法和答案
已经有11人回复
- 【求助完毕】为何要计算空带？
已经有5人回复
- 【原创】第一原理计算的一些心得
已经有160人回复
- 【求助】一个二阶常微分方程的解法（重奖金币！）
已经有16人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】
已经有15人回复
- 【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
- 【重点讨论】Castep中的Empty Band的用处和设置？（参与讨论就有机会赢取大礼包！）
已经有54人回复
» 抢金币啦！回帖就可以得到:
-
助理研究员——吉林大学-苏州国家实验室光电半导体智能设计与应用团队
+1/479
-
诚征女友
+1/259
-
华东理工大学 生物反应器工程全国重点实验室 博士后招聘
+1/79
-
Postdoctoral fellow in Medicinal Chemistry, McGill University
+1/71
-
中科院东莞材料所依托中科大招收2027届推免生
+1/50
-
没基金，自费做nature课题值得吗?
+1/44
-
西安交通大学补亚忠课题组招聘2027年博士生
+1/27
-
直博夏令营--北京理工大学/国家人工智能学院AI交叉学科联合招生
+1/21
-
坐标北京，诚征女友
+2/12
-
招生宣讲会倒计时！7月29日20：00，不见不散！
+1/10
-
浙江大学王宇课题组诚聘博士后/博士研究生/科研助理
+1/9
-
四川大学周加境课题组招博士研究生/研究助理（自组装材料/生物医学方向）
+1/6
-
浙江大学医学院附属妇产科医院张霈婧课题组招科研助理——生物信息方向
+2/6
-
新西兰坎特伯雷大学招博士后
+1/6
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘（长期有效）
+5/5
-
多本Scopus收录国际英文期刊征稿，全免 APC！含医学期刊
+1/5
-
给合作老师打个广告，找个做 NMR/PHIP 超极化的 CSC 博士生
+1/5
-
（Scopus期刊）Biofunctional Materials《生物功能性材料》招募青年编委及征稿
+1/4
-
中山大学孙逸仙纪念医院潘越教授团队招聘博士后(40万年薪)和科研助理
+1/3
-
iThenticate查SCI论文，一般查几次？有没有技巧一次过？
+1/1
★ ★ ★
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
