| 查看: 2988 | 回复: 17 | ||||
xirainbow木虫 (正式写手)
|
[交流]
【原创】The meaning of nband 已有1人参与
|
|
MS版上发了个贴，讨论nband的意义 PWscf也牵涉到这个问题。 以下为我自己的理解，其中最后一段有待商榷 原帖网址：http://blog.sina.com.cn/s/blog_5f15ead20100e1l5.html 由于上传图片比较麻烦 所以原帖中的很多公式我没有上传 有兴趣的可以去博客上看；） This is my second try. I spend almost two days to get some progress in understanding of Davidson diagonalization. I am not interested in Davidson method itself. But without knowing it, you can not understand the real meaning of parameter "nband". 在计算中指定了截止动能，对于一个固定的k点，也就意味着指定了有多少个G（=2pi/a的整数倍）参与傅立叶变换（假设为M个）。而G的个数就是相应哈密顿量的维度，M*M维。通常M是非常大的，而我们只需要知道较低的几个被占据的能级，即可求出总能等一些列刚兴趣的物理量。因此直接对M*M维的哈密顿量进行对角化是不合适的。所以需要采用迭代的方法来求解几个低能的本征值和本征能量即可。 其实迭代方法在矩阵对角化中是非常常用的数值计算方法，他可以通过迭代的方式只求出某一个或者几个本征值。Davidson是其中一个非常出名的，因为他在量子力学中取得了很大的成功，这个成功来源于量子力学的哈密顿矩阵通常都是稀疏的，就是非对角元很多都是0，或者接近0。 当nband被指定后，本征方程的就被给定。。本征方程的维度N应该是等于nband的。此时N远小于有ecutwfc决定的G的个数，即H的维度M*M。而波函数psi则是由M个平面波的叠加。 然后求得对初始猜测波函数的修正。这样的修正波函数作为独立的基矢，加入到初始猜测的N个波函数上，构成N+1维的子空间。如果对N个初始猜测波函数均采用这种办法，则得到2*N维的子空间。然后，在2*N维的子空间中求得H矩阵，并对其对角化，比较前N个本征值是否有变化，判定是否收敛。 (本征方程的维度是否在不断扩大，有待商榷) 纳米盘下载： http://d.namipan.com/sd/1558396 [ Last edited by wuli8 on 2009-12-22 at 19:13 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理计算 |
» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有13人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有11人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有9人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【分享】《The Practice of Programming》(应该是编程必备吧，程序设计实践)
已经有147人回复
- what's the meaning of this sentence?
已经有8人回复
- what's the meaning of this sentence?
已经有3人回复
- 【Pure water】The meaning of life
已经有19人回复
- 【Discussion】What's the meaning of this sentence ?
已经有14人回复
- 【Original】The meaning of life
已经有14人回复
- The Old Woman and the Wine-Jar, what is the meaning of this fable?
已经有4人回复
- the meaning of the word polymerosomes?
已经有3人回复
- 【Pray】the meaning of the word polymerosomes？
已经有13人回复
- 【Activity】Do you know what's the meaning of O T L!
已经有36人回复
- 【交流】what is the meaning of 'target exposure'
已经有3人回复
- what's the meaning of ‘running water’
已经有7人回复
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2009-12-22 17:17:38
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
3楼2009-12-22 17:21:06
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
4楼2009-12-22 17:43:52
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
5楼2009-12-22 17:45:05
6楼2009-12-22 18:58:05
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
7楼2009-12-23 14:01:05
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
8楼2010-02-09 19:19:25
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
9楼2010-02-09 22:44:31
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
10楼2010-02-10 09:12:15
