版块导航
正在加载中...
客户端APP下载
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(3264)
>
虫友互识
(378)
>
文献求助
(353)
>
导师招生
(246)
>
考博
(153)
>
硕博家园
(128)
>
休闲灌水
(88)
>
博后之家
(81)
>
招聘信息布告栏
(54)
>
基金申请
(52)
>
论文投稿
(51)
>
仿真模拟
(46)
>
攻关文献(高奖励)
(38)
>
绿色求助(高悬赏)
(25)
>
教师之家
(22)
>
找工作
(19)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
第一性原理
»
Vasp&MedeA
»
【求助】K-S 方程的迭代解法
11
1/1
返回列表
查看: 1762 | 回复: 10
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
ss_kampf
铜虫
(初入文坛)
应助: 0
(幼儿园)
金币: 78
帖子: 18
在线: 22.2小时
虫号: 1119678
[交流]
【求助】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
]
回复此楼
» 猜你喜欢
求碳排放博导;方向是LCA、生命周期可持续发展以及碳排放
已经有5人回复
26年博士申请自荐-电催化
已经有7人回复
2026年博士申请求捞
已经有3人回复
2026博士或科研助理转27年博士
已经有6人回复
国自科送审了吗
已经有11人回复
博士招生
已经有5人回复
夜,静悄悄的
已经有8人回复
研究生做的很差,你们会让毕业吗?
已经有6人回复
26年申博自荐-计算机视觉
已经有4人回复
高级回复
» 本主题相关价值贴推荐,对您同样有帮助:
有个求解微分方程用到高斯赛德尔迭代的,算不太明白啊,各位帮忙看看
已经有5人回复
指数方程的解法
已经有4人回复
都不知道怎么命名标题:迭代法解方程组?
已经有16人回复
【求助】求助一个二元二次回归方程的解法和答案
已经有11人回复
【求助完毕】为何要计算空带?
已经有5人回复
【原创】第一原理计算的一些心得
已经有160人回复
【求助】一个二阶常微分方程的解法(重奖金币!)
已经有16人回复
非线性方程组的迭代法(数值计算高手请进)
已经有7人回复
【求助】求助:vb编程中用牛顿迭代解三次方程为什么只得到一个根?【已完成】
已经有15人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
【重点讨论】Castep中的Empty Band的用处和设置?(参与讨论就有机会赢取大礼包!)
已经有54人回复
» 抢金币啦!回帖就可以得到:
查看全部散金贴
香港城市大学 数据科学系+能源与环境学院 李堂华 课题组诚招全奖博士/博士后/研究助理
+
3
/180
坐标广州,征女友
+
2
/124
河南师范大学水产学院博士招生
+
1
/82
太行山徒步之凌水河(一)
+
1
/79
山东征女友,坐标济南
+
1
/70
昆明理工大学交通载运车辆方向招收2026年秋季入学博士生2名
+
4
/44
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招2026级博士(5月15号截止)
+
1
/31
坐标南京
+
1
/20
南京邮电大学李巍教授招收2026博士生(5月12日前有效)
+
1
/18
2026年西安交通大学博士招生!
+
1
/16
港科大广州 - 量子计算全奖PhD招生
+
1
/11
申请账号解封
+
1
/10
论坛管理,垃圾帖子充满整个论坛
+
1
/9
注册小木虫周年大礼包
+
1
/5
易度质量流量计:燃料电池电堆气密性检测的高精度之选
+
1
/5
招收2026年秋季入学博士生1名(河北工业大学/北京科技大学联合 增材制造/生物材料)
+
1
/4
UIO-66-NH2/RGO如何碳化,求助!!!!!
+
1
/3
【急招】北京航空航天大学 | 超声脑机接口方向 | 2026年博士研究生
+
1
/3
李老师课题组急招收生物信息学方向博士研究生!
+
1
/1
知名外资仪器厂家急招Application Scientist(2名
+
1
/1
1楼
2011-03-02 14:02:14
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
zzy870720z
荣誉版主
(文坛精英)
1ST强帖: 3
应助: 47
(小学生)
贵宾: 9.05
金币: 30914.3
帖子: 12592
在线: 23567.6小时
虫号: 745488
★ ★ ★
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
铜虫
(初入文坛)
应助: 0
(幼儿园)
金币: 78
帖子: 18
在线: 22.2小时
虫号: 1119678
我才发现问了一个很白痴的问题,看来晚上吃饭时得补补智商了~
赞
一下
回复此楼
3楼
2011-03-02 18:03:21
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
franch
版主
(著名写手)
1ST强帖: 2
应助: 62
(初中生)
贵宾: 0.185
金币: 7766.8
帖子: 1964
在线: 731.1小时
虫号: 587928
关注此贴
回复此楼
4楼
2011-03-02 18:51:18
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
cenwanglai
荣誉版主
(知名作家)
1ST强帖: 5
应助: 46
(小学生)
贵宾: 8.842
金币: 7441.4
帖子: 5306
在线: 1961.4小时
虫号: 537452
★ ★
小木虫(金币
+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
荣誉版主
(文坛精英)
1ST强帖: 3
应助: 47
(小学生)
贵宾: 9.05
金币: 30914.3
帖子: 12592
在线: 23567.6小时
虫号: 745488
★ ★
小木虫(金币
+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
铜虫
(初入文坛)
应助: 0
(幼儿园)
金币: 78
帖子: 18
在线: 22.2小时
虫号: 1119678
引用回帖:
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
铜虫
(初入文坛)
应助: 0
(幼儿园)
金币: 78
帖子: 18
在线: 22.2小时
虫号: 1119678
引用回帖:
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
金虫
(著名写手)
1ST强帖: 5
应助: 2
(幼儿园)
贵宾: 0.446
金币: 431.9
帖子: 1207
在线: 114.8小时
虫号: 487243
★
小木虫(金币
+0.5
):给个红包,谢谢回帖交流
何为“只有迭代的方法”,迭代的目的不是为了能够对角化吗?
赞
一下
(1人)
回复此楼
9楼
2011-03-02 22:29:30
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
ss_kampf
铜虫
(初入文坛)
应助: 0
(幼儿园)
金币: 78
帖子: 18
在线: 22.2小时
虫号: 1119678
★
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
金虫
(著名写手)
1ST强帖: 5
应助: 2
(幼儿园)
贵宾: 0.446
金币: 431.9
帖子: 1207
在线: 114.8小时
虫号: 487243
★
小木虫(金币
+0.5
):给个红包,谢谢回帖交流
不是还有迭代对角化吗?
赞
一下
(1人)
回复此楼
11楼
2011-03-03 09:00:19
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
相关版块跳转
第一性原理
量子化学
计算模拟
分子模拟
仿真模拟
程序语言
我要订阅楼主
ss_kampf
的主题更新
11
1/1
返回列表
如果回帖内容含有宣传信息,请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
高级回复
(可上传附件)
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
信息提示
关闭
请填处理意见
关闭
确定