| 查看: 1570 | 回复: 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 这是什么原因,请各位大神解答一下... |
» 猜你喜欢
为什么中国大学工科教授们水了那么多所谓的顶会顶刊,但还是做不出宇树机器人?
已经有9人回复
版面费该交吗
已经有9人回复
体制内长辈说体制内绝大部分一辈子在底层,如同你们一样大部分普通教师忙且收入低
已经有13人回复
面上可以超过30页吧?
已经有4人回复
“人文社科而论,许多学术研究还没有达到民国时期的水平”
已经有5人回复
什么是人一生最重要的?
已经有4人回复
今年春晚有几个节目很不错,点赞!
已经有12人回复

pangrui1985
铜虫 (小有名气)
- 应助: 25 (小学生)
- 金币: 23.2
- 红花: 6
- 帖子: 142
- 在线: 89.2小时
- 虫号: 2808500
- 注册: 2013-11-18
- 专业: 凝聚态物性I:结构、力学和
2楼2015-07-29 14:50:31

3楼2015-07-29 16:04:16













回复此楼