| 查看: 1264 | 回复: 2 | |||
[交流]
【求助】截断能测试不收敛是怎么回事 已有1人参与
|
|
如题，我用pwscf测试了BaTiO3的截断能，结果发现截断能很难收敛，基本要到450Ry才能达到1.0e-4Ry精度，我同样测试了其它参数的影响，发现不论怎么换其他参数都是这个结果。 测试脚本： #!/bin/sh for i in 60 70 80 90 100 110 120 130 140 ...... 570 do cat > scf.in < &control calculation= 'scf' restart_mode='from_scratch' pseudo_dir = './' outdir ='./' prefix = 'BaTiO3' tprnfor = .true. / &system ibrav = 1, celldm(1)=7.57, nat=5, ntyp = 3, ecutwfc = $i, ecutrho = 600,occupations = 'smearing', smearing = 'gaussian', degauss = 0.01 / &electrons diagonalization = 'david' conv_thr = 1.0e-8 mixing_beta = 0.3 / ATOMIC_SPECIES Ba 137.327 Ba.pbe-nsp-van.UPF Ti 47.867 Ti.pbe-sp-van_ak.UPF O 15.9994 O.pbe-van_ak.UPF ATOMIC_POSITIONS Ba 0 0 0 Ti 0.5 0.5 0.5 O 0.5 0.5 0 O 0.5 0 0.5 O 0 0.5 0.5 K_POINTS AUTOMATIC 7 7 7 0 0 0 ! pw.x -inp scf.in E=`grep 'Harris-Foulkes estimate ' output.* |tail -1 |awk '{printf "%7.6f \n", $4}'` echo $i $E >>Evk.dat done 截断能再大，到580已算不下去了，请问这个是怎么回事？为什么很难收敛？ 难道是赝势里Ba为非线性nsp赝势原因还是什么？按理说USPP赝势截断能应该很小啊 [ Last edited by zhang668 on 2010-10-27 at 23:10 ] |
回复此楼» 猜你喜欢
能否退出参与的面上项目解除限项
已经有24人回复
-
号外号外
已经有8人回复
-
filecode，4个jtjc了
已经有17人回复
-
人气不行了
已经有9人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有5人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于收敛性测试问题
已经有6人回复
- 收敛测试，这种情况我该继续测试么,简直就是跳跃式的嘛
已经有24人回复
- 关于dmol3的收敛性测试
已经有12人回复
- 难道每次计算前都咬进行收敛性测试吗？
已经有14人回复
- 【求助】K点取得越大越不收敛，什么原因！
已经有20人回复
- 【求助成功】castep 收敛性测试
已经有22人回复
- 【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
已经有14人回复
- 【其他】关于收敛测试的问题
已经有22人回复
- 【求助】收敛测试怎么做？？？
已经有17人回复
- 【求助】请问收敛测试的标准
已经有15人回复
- 【求助】收敛性测试的问题
已经有40人回复
- 【求助】castep收敛性测试问题
已经有20人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
- 【求助】是不是什么计算都得进行收敛测试？
已经有8人回复
2楼2010-10-27 12:02:14
ylli
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1558.8
- 红花: 3
- 帖子: 267
- 在线: 41.6小时
- 虫号: 277021
- 注册: 2006-09-03
- 专业: 凝聚态物性 II ：电子结构
3楼2010-10-27 12:33:36
20