| 查看: 1673 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
加入“nspin、 starting_magnetization”关键词后scf过程不收敛，求解释 已有1人参与
|
||
|
对同一个模型，不加入“nspin、 starting_magnetization”关键词时候scf过程收敛，但是对所有原子考虑自旋，加入关键词后scf不收敛了。我的输入文件如下： &control calculation='scf' restart_mode='from_scratch', pseudo_dir = './', outdir='./' prefix='NiOH2' tstress = .true. tprnfor = .true. / &system ibrav = 14, celldm(1) =11.829688914, celldm(2) =1.0, celldm(3) =3.194888179, celldm(4) =0.0, celldm(5) =0.0, celldm(6) =-0.5, nat= 20, ntyp= 3, ecutwfc = 25.0, ecutrho = 300.0 occupations='smearing', smearing='gaussian', degauss=0.02 nspin=2 starting_magnetization(1)=0.1 starting_magnetization(2)=0.1 starting_magnetization(3)=0.1 / &electrons diagonalization='david' conv_thr = 1.0e-8 mixing_beta = 0.7 / ATOMIC_SPECIES H 1.0 H.pbe-kjpaw_psl.0.1.UPF O 16.0 O.pbe-n-kjpaw_psl.0.1.UPF Ni 59.0 Ni.pbe-n-kjpaw_psl.0.1.UPF ATOMIC_POSITIONS crystal H 0.1666666666666667 0.3333333333333334 0.4699999988079069 H 0.3333333333333334 0.1666666666666667 -0.4699999988079068 H 0.6666666666666669 0.3333333333333335 0.4699999988079069 H 0.8333333333333334 0.1666666666666669 -0.4699999988079068 H 0.1666666666666667 0.8333333333333334 0.4699999988079069 H 0.3333333333333334 0.6666666666666669 -0.4699999988079068 H 0.6666666666666669 0.8333333333333336 0.4699999988079069 H 0.8333333333333334 0.6666666666666670 -0.4699999988079068 O 0.1666666666666667 0.3333333333333334 0.2399999946355819 O 0.3333333333333334 0.1666666666666667 -0.2399999946355818 O 0.6666666666666669 0.3333333333333335 0.2399999946355819 O 0.8333333333333334 0.1666666666666669 -0.2399999946355818 O 0.1666666666666667 0.8333333333333334 0.2399999946355819 O 0.3333333333333334 0.6666666666666669 -0.2399999946355818 O 0.6666666666666669 0.8333333333333336 0.2399999946355820 O 0.8333333333333334 0.6666666666666670 -0.2399999946355818 Ni 0.0000000000000000 0.0000000000000000 0.0000000000000000 Ni 0.5000000000000001 0.0000000000000001 0.0000000000000000 Ni 0.0000000000000000 0.5000000000000000 0.0000000000000000 Ni 0.5000000000000001 0.5000000000000002 0.0000000000000001 K_POINTS (automatic) 8 8 8 0 0 0 同时终端出现： [customer@node01 ~]$ mpirun -np 16 pw.x<NiOH2.scf.in>NiOH2.scf.out 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 这是什么原因，请各位大神解答一下... |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
3楼2015-07-29 16:04:16
pangrui1985
铜虫 (小有名气)
- 应助: 25 (小学生)
- 金币: 23.2
- 红花: 6
- 帖子: 142
- 在线: 89.2小时
- 虫号: 2808500
- 注册: 2013-11-18
- 专业: 凝聚态物性I:结构、力学和
2楼2015-07-29 14:50:31
20