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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】K-S 方程的迭代解法
111/1返回列表
查看: 1992  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 让我中一个面上吧！ +7 大萍1987 2026-08-20 8/400 2026-08-20 21:19 by 大萍1987
[基金申请] 科研孤儿太难了 +12 我4大白菜 2026-08-20 13/650 2026-08-20 20:55 by gltch
[基金申请] 言之凿凿，明天官网会通知放榜？ +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗？20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭