应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】K-S 方程的迭代解法
北京石油化工学院2026年研究生招生接收调剂公告
111/1返回列表
查看: 1732  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ss_kampf

铜虫 (初入文坛)

[交流] 【求助】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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-02 14:02:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)
★ ★ ★
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一般都是用平面波基组或是原子轨道基组做展开
赞一下(1人) 回复此楼
2楼2011-03-02 16:46:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ss_kampf

铜虫 (初入文坛)
我才发现问了一个很白痴的问题，看来晚上吃饭时得补补智商了～
赞一下 回复此楼
3楼2011-03-02 18:03:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

franch

版主 (著名写手)
关注此贴
回复此楼
4楼2011-03-02 18:51:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bingmou(金币+1): 这个在electronic structure里面都有讲到 2011-03-02 22:21:01
呵呵，这个过程我也没懂的。呵呵。比如为什么是子空间(subspace)对角化？AX=B，A对角化后，X就可以解出来了，但是为什么是子空间对角化呢？

H应该是考虑赝势和交换关联泛函的“人造”算符。第二步构造出新的波函数后，直接得到新的电荷密度，与上一步电荷密度混合后就得到下一步输入电荷密度。对角化看起来是可以跳过的？
赞一下(2人) 回复此楼
5楼2011-03-02 18:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bingmou(金币+1): 感谢交流 2011-03-02 22:21:15
引用回帖:
Originally posted by cenwanglai at 2011-03-02 18:55:33:
呵呵，这个过程我也没懂的。呵呵。比如为什么是子空间(subspace)对角化？AX=B，A对角化后，X就可以解出来了，但是为什么是子空间对角化呢？

H应该是考虑赝势和交换关联泛函的“人造”算符。第二步构造出新的波 ...

密度不是唯一的判别标准
scf也有收敛标准，比如能量，对角化就可以得到能量等其他性质的数据，还是需要对角化的
赞一下(2人) 回复此楼
6楼2011-03-02 19:36:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ss_kampf

铜虫 (初入文坛)
引用回帖:
Originally posted by zzy870720z at 2011-03-02 19:36:42:
密度不是唯一的判别标准
scf也有收敛标准，比如能量，对角化就可以得到能量等其他性质的数据，还是需要对角化的

彻底糊涂了。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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ss_kampf

铜虫 (初入文坛)
引用回帖:
Originally posted by cenwanglai at 2011-03-02 18:55:33:
呵呵，这个过程我也没懂的。呵呵。比如为什么是子空间(subspace)对角化？AX=B，A对角化后，X就可以解出来了，但是为什么是子空间对角化呢？

H应该是考虑赝势和交换关联泛函的“人造”算符。第二步构造出新的波 ...

对啊，用不用对角化呢？
赞一下 回复此楼
8楼2011-03-02 22:22:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
何为“只有迭代的方法”，迭代的目的不是为了能够对角化吗？
赞一下(1人) 回复此楼
9楼2011-03-02 22:29:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ss_kampf

铜虫 (初入文坛)
★
youzhizhe(金币+1): 感谢回复虫友问题，多多交流。 2011-03-03 00:26:02
引用回帖:
Originally posted by bingmou at 2011-03-02 22:29:30:
何为“只有迭代的方法”，迭代的目的不是为了能够对角化吗？

看样子，只迭代，不对角化也可以的。
赞一下(1人) 回复此楼
10楼2011-03-02 23:10:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
不是还有迭代对角化吗？
赞一下(1人) 回复此楼
11楼2011-03-03 09:00:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ss_kampf 的主题更新
111/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿安徽大学计算机科学与技术学硕，331分求调剂 +3 蒋昌鹏qtj 2026-04-01 3/150 2026-04-01 23:56 by 义文wang
[考研] 材料科学与工程调剂 +16 深V宿舍吧 2026-03-30 17/850 2026-04-01 21:00 by 赖春艳
[考研] 一志愿郑大材料工程290求调剂 +18 Youth_ 2026-03-30 18/900 2026-04-01 19:23 by 无际的草原
[考研] 省双一流重点一本大学招收调剂 +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[考研] 一志愿同济大学323分（080500）求调剂 +4 yikeniu 2026-04-01 4/200 2026-04-01 14:06 by asdfzly
[考研] 349求调剂 +6 zwjjjjjj 2026-03-31 6/300 2026-04-01 09:16 by JourneyLucky
[考研] 085601英二数二求调剂 总分325 +4 余航航 2026-03-31 4/200 2026-03-31 17:38 by 唐沐儿
[考研] 材料工程专硕求调剂 +10 hyl3153942 2026-03-29 10/500 2026-03-31 16:31 by hypershenger
[基金申请] 面上5B能上会吗？ +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 266分，求材料相关专业调剂 +10 哇呼哼呼哼 2026-03-30 12/600 2026-03-31 11:00 by 熊一刀
[考研] 22408 359分调剂 +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[考研] 08工科求调剂286 +5 tgs_001 2026-03-28 5/250 2026-03-31 08:18 by 一只好果子?
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
[考研] 抱歉 +4 田洪有 2026-03-30 4/200 2026-03-30 21:26 by mumin1990
[考研] 一志愿华中师范化学332分求调剂 +3 Lyy930824@ 2026-03-29 3/150 2026-03-30 20:15 by DHUSHUAI
[考研] 105500药学求调剂，一志愿山东大学药学，348分 +3 gr哈哈哈 2026-03-28 3/150 2026-03-30 18:56 by 源_2020
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 一志愿吉大071010，316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
信息提示
关闭
请填处理意见
关闭