24小时热门版块排行榜    

查看: 1691  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 找导师 +5 yuanjiabao 2026-08-29 6/300 2026-08-30 11:52 by zhouyanli11
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 11:32 by l0VvVHGBGRLv
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:07 by l0VvVHGBGRLv
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 06:11 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 05:40 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +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
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见