版块导航
正在加载中...
客户端APP下载
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(1286)
>
虫友互识
(80)
>
基金申请
(16)
>
博后之家
(5)
>
导师招生
(4)
>
硕博家园
(3)
>
考博
(2)
>
公派出国
(2)
>
考研
(2)
>
休闲灌水
(2)
>
海归之家
(1)
>
计算模拟
(1)
>
晶体
(1)
>
教师之家
(1)
>
有机交流
(1)
>
找工作
(1)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
第一性原理
»
Vasp&MedeA
»
【求助】K-S 方程的迭代解法
11
1/1
返回列表
查看: 1923 | 回复: 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
]
回复此楼
» 猜你喜欢
今年E04面上
已经有23人回复
生命口会评
已经有6人回复
27届辽宁大学应届毕业生申博
已经有3人回复
关于如何从代码看上不上会
已经有23人回复
大龄残疾硕士的一点执念
已经有25人回复
chemdraw
已经有6人回复
祈祷青基必中
已经有10人回复
b口会评
已经有6人回复
信息学部会评时间
已经有3人回复
在职考研
已经有7人回复
高级回复
» 本主题相关价值贴推荐,对您同样有帮助:
有个求解微分方程用到高斯赛德尔迭代的,算不太明白啊,各位帮忙看看
已经有5人回复
指数方程的解法
已经有4人回复
都不知道怎么命名标题:迭代法解方程组?
已经有16人回复
【求助】求助一个二元二次回归方程的解法和答案
已经有11人回复
【求助完毕】为何要计算空带?
已经有5人回复
【原创】第一原理计算的一些心得
已经有160人回复
【求助】一个二阶常微分方程的解法(重奖金币!)
已经有16人回复
非线性方程组的迭代法(数值计算高手请进)
已经有7人回复
【求助】求助:vb编程中用牛顿迭代解三次方程为什么只得到一个根?【已完成】
已经有15人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序
已经有10人回复
【重点讨论】Castep中的Empty Band的用处和设置?(参与讨论就有机会赢取大礼包!)
已经有54人回复
» 抢金币啦!回帖就可以得到:
查看全部散金贴
实验室技术孵化和中试技术转让
+
1
/1970
中山大学医学院(深圳)招聘科研助理(生物信息学方向)
+
1
/100
热分析测试
+
1
/84
北京师范大学(珠海)电催化方向招收2027年硕士研究生化学或材料方向
+
1
/33
上海交通大学-化学化工学院 博士后招聘
+
1
/31
美国圣路易斯大学招聘生物信息学博士后
+
1
/29
北京理工大学长三角研究院(嘉兴)先进材料与智能装备团队 招聘
+
1
/29
北京师范大学(珠海)电催化方向招收2027年硕士研究生化学或材料方向
+
1
/27
干细胞毒理学研究组招收2027年硕博
+
1
/11
澳大利亚昆士兰大学(UQ)】全奖PhD招生:碳捕集与液流电池储能方向
+
1
/10
南方科技大学医学院-神经胶质细胞课题组招收联培(客座)博士生
+
1
/7
悉尼科技大学(QS 第87)数据科学/AI招收2027年秋季入学全奖及CSC博士生
+
1
/6
评职称 / 毕业赶进度,期刊发表选刊、审稿节奏问题可交流答疑
+
1
/5
双一流高校招收博士(现代技术管理/环境技术经济管理)及硕士研究生(长期有效)
+
1
/4
清华大学深圳国际研究生院宋桥课题组招收2027年博士研究生
+
1
/4
浙江大学-化工学院刘平伟课题组-二维材料/功能聚合物开发-博后招聘
+
1
/3
暨南大学开明峰教授团队博士后招聘启事
+
1
/3
诚聘智慧农业方向博士教师
+
1
/2
清华大学韩军功教授课题组博士后招聘启事
+
1
/1
武汉理工大学刘鹏飞院士团队招聘科研助理
+
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云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
信息提示
关闭
请填处理意见
关闭
确定