24小时热门版块排行榜    

查看: 1575  |  回复: 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的回帖

pangrui1985

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
正常啊,加了自旋自由度更多,更难收敛,此时一般要改小mixing_beta
另外,我敢肯定你的截断能不够,你用的氧赝势的推荐截断能是47李德宝,所以就更难收敛了
2楼2015-07-29 14:50:31
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 赵奇一 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[考博] 申博26年 +4 八6八68 2026-03-19 4/200 2026-03-24 15:49 by 小Ben呵呵
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 327求调剂 +5 prayer13 2026-03-23 5/250 2026-03-23 22:11 by 星空星月
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
信息提示
请填处理意见