24小时热门版块排行榜    

查看: 3730  |  回复: 17

xirainbow

木虫 (正式写手)

★ ★
aylayl08(金币+2):谢谢交流讨论 2010-02-18 09:35
引用回帖:
Originally posted by goldenfisher at 2010-02-10 09:12:15:
nbnd的个数必须远远小于平面波个数,通常来说满足,但是很多时候还是不满足的,比如高温的时候。

你可以看看pwscf的mail-list,很多人也在讨论这个问题

多谢您的指正:)
本人孤陋寡闻,没用过分子动力学
所以只知道通常的DFT只能计算0K的性质:)

davidson使用的子空间迭代法
如果[哈密顿量的维度/子空间的维度]过小,
那么少数几次迭代后,子空间将搜索完整个空间
而此时还没有收敛
因此会出问题:(
但其中更细节的数学知识我不懂:(

[ Last edited by xirainbow on 2010-2-10 at 11:20 ]
11楼2010-02-10 11:02:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1):谢谢讨论,欢迎常来 2010-02-18 09:34
引用回帖:
Originally posted by xirainbow at 2010-02-10 11:02:33:

多谢您的指正:)
本人孤陋寡闻,没用过分子动力学
所以只知道通常的DFT只能计算0K的性质:)

davidson使用的子空间迭代法
如果[哈密顿量的维度/子空间的维度]过小,
那么少数几次迭代后,子空间将搜索完 ...

客气客气,我经常在pwscf的mailist里面发现你的身影,哈哈
DFT有含温的DFT。其中很多东西需要详细讨论。
12楼2010-02-10 22:49:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awmc2008

至尊木虫 (文坛精英)

学习者

xuexizhe
13楼2010-02-11 11:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

引用回帖:
Originally posted by goldenfisher at 2010-02-10 22:49:38:


客气客气,我经常在pwscf的mailist里面发现你的身影,哈哈
DFT有含温的DFT。其中很多东西需要详细讨论。

不是身影,是呻吟:)
你说的含温的DFT是什么呢?:)
14楼2010-02-11 21:35:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guoguo3138

银虫 (小有名气)

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
一般为了计算的准确性应多设置几个空带数,但是会加大计算量
15楼2010-02-12 14:09:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+2):谢谢解答讨论 2010-02-18 09:37
含温的DFT可以参考文献Mermin, Phys. Rev. 137, A1441.这是当今温度较高的分子动力学的基础。增加空带数并不只是为了增加准确性,而是在高温的时候必须设置足够多的带数,才能考虑到足够高的能级分布。计算量是与带数成3次方的关系的,这也是当今第一原理的一个难点之一。
16楼2010-02-12 23:39:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

刚从vasp的说明书中看到了nband的一些有意义的信息。具体内容如下:
In iterative matrix-diagonalization schemes eigenvectors close to the top of the calculated number of vectors converge much slower than the lowest eigenvectors. This might result in a significant performance loss if not enough empty bands are included in the calculation.
Therefore we recommend to set NBANDS to NELECT/2 + NIONS/2. This setting is safe in most cases. In some cases, it is also possible to decrease the number of additional bands to NIONS/4 for large systems without performance loss, but on the other hand transition metals do require a much larger number of empty bands (up to 2*NIONS).
To check this parameter perform several calculations for a fixed potential with an increasing number of bands (e.g. starting from NELECT/2 + NIONS/2). An accuracy of 10−6 should be obtained in 10-15 iterations. Mind that the RMM-DIIS scheme is more sensible to the number of bands than the default CG algorithm.
由于时间过久,原帖已禁止我改动。希望管理员可以将这些补充加到原帖后面,谢谢:)

[ Last edited by xirainbow on 2010-3-30 at 18:19 ]
17楼2010-03-30 18:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

风轻云淡7

金虫 (小有名气)

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
nbnd一般是要大于电子数的一般的,因为那样是不能输出费米面的……所以楼主还是去大一点好,但是最好不要大于价电子书数……
18楼2010-08-29 21:32:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xirainbow 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 08:24 by byEJ59EAM6hR
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见