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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】K-S 方程的迭代解法
51/1返回列表
查看: 1996  |  回复: 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的回帖

ss_kampf

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

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的回帖

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的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 让我中一个面上吧！ +8 大萍1987 2026-08-20 10/500 2026-08-20 21:48 by 大萍1987
[基金申请] 言之凿凿，明天官网会通知放榜？ +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
[基金申请] 时间戳变了，能看出什么问题？ +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
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +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
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
关闭
请填处理意见
关闭