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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】K-S 方程的迭代解法
51/1返回列表
查看: 1674  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭