24小时热门版块排行榜    

查看: 1692  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 赵奇一 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 11:32 by l0VvVHGBGRLv
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:30 by l0VvVHGBGRLv
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:18 by l0VvVHGBGRLv
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:10 by l0VvVHGBGRLv
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见