24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

赵奇一

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 1311
帖子: 172
在线: 187.9小时
虫号: 2798451
注册: 2013-11-13
性别: GG
专业: 粒子物理学和场论

[求助] 加入“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
这是什么原因，请各位大神解答一下...

回复此楼