24小时热门版块排行榜    

查看: 1576  |  回复: 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
这是什么原因,请各位大神解答一下...
回复此楼
天使可以飞起来是因为他们把自己看得很轻
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵奇一

新虫 (小有名气)

引用回帖:
2楼: Originally posted by pangrui1985 at 2015-07-29 14:50:31
正常啊,加了自旋自由度更多,更难收敛,此时一般要改小mixing_beta
另外,我敢肯定你的截断能不够,你用的氧赝势的推荐截断能是47李德宝,所以就更难收敛了

谢谢你的回复,那请问你有啥好的建议没?我新手,什么都在慢慢琢磨
天使可以飞起来是因为他们把自己看得很轻
3楼2015-07-29 16:04:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

pangrui1985

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
正常啊,加了自旋自由度更多,更难收敛,此时一般要改小mixing_beta
另外,我敢肯定你的截断能不够,你用的氧赝势的推荐截断能是47李德宝,所以就更难收敛了
2楼2015-07-29 14:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] B区考研调剂 +4 yqdszhdap- 2026-03-22 5/250 2026-03-25 08:51 by baoball
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 一志愿哈工大,085400,320,求调剂 +3 gdlf9999 2026-03-24 3/150 2026-03-24 10:08 by 搏击518
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +5 葵梓卫队 2026-03-18 7/350 2026-03-23 16:26 by lingjue
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
信息提示
请填处理意见