| 查看: 1689 | 回复: 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 这是什么原因，请各位大神解答一下... |
回复此楼» 猜你喜欢
找导师
已经有5人回复
-
国社科又开始会评了，不知道这次命运如何
已经有3人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
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
30