| 查看: 1587 | 回复: 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 这是什么原因，请各位大神解答一下... |
回复此楼
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
