24小时热门版块排行榜    

查看: 2902  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见