24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者非天之翼蓝羽的 5 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

非天之翼蓝羽

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 794.9
帖子: 44
在线: 1.6小时
虫号: 263537
注册: 2006-07-02
性别: MM
专业: 复合材料

[交流] 【求助】RF计算出错，参数设置问题？

以下为输入文件：
ndtset 2
# Set 1 :
  ionmov1 2
ntime1 5
  tolmxf1 5.0d-5
  natfix1 8
  iatfix1 1 1 2 3 3 3 3 3
# Set 2 :
dilatmx2 1.5
getxred2 -1
  getwfk2 -1
  ionmov2 2
ntime2 10
optcell2 2
  tolmxf2 5.0d-5
strfact2 100
strprecon 0.1
  natfix2 8
  iatfix2 1 1 2 3 3 3 3  3
#Common input data
  acell 11.1275  22.2447  15.5455 ＃Bhor单位
  angdeg  90 90 90
spgroup 55
ntypat 3
  znucl 82 40 8
  natom 8
  typat 1 1 2 3 3 3 3 3
xred 0.7085  0.1300  0.0000
      0.7117  0.1246  0.5000
      0.2423  0.1242  0.2536
      0.2912  0.1009  0.0000
      0.2761  0.1585  0.5000
      0.0284  0.2606  0.2207
      0.5000  0.0000  0.3022
      0.0000  0.0000  0.2822
  nband 25
ecut 8.0
ecutsm 0.5
kptopt 1
  ngkpt 4 4 4
nshiftk 1
shiftk 0.5 0.5 0.5
diemac 3.0
iscf 5
  nstep 60
tolvrs 1.0d-14

这几天进行了粗浅的尝试和反思，但运行还是出现了问题：

首先是软件貌似不识别我输入的参数，acell，angdeg，xred三者应该可以确定空间群了吧？我想应该输出的是pbam，但是output中：
DATASET 1 : space group P1 (#  1); Bravais aP (primitive triclinic)

另外输出文件还有warning：
ntime= 5 was not enough Broyd/MD steps to converge gradients:
  max grad (force/stress) = 2.9312E-02 > tolmxf= 5.0000E-05 ha/bohr (free atoms)
nstep= 60 was not enough SCF cycles to converge;
  potential residual=  4.822E-10 exceeds tolvrs=  1.000E-14
ntime= 10 was not enough Broyd/MD steps to converge gradients:
  max grad (force/stress) = 5.9921E-02 > tolmxf= 5.0000E-05 ha/bohr (free atoms)
不知道在这种情况下是应该增加ntime，nstep还是降低tolvrs跟tolmxf呢

?

求教了。。。

[ Last edited by 非天之翼蓝羽 on 2010-3-29 at 22:19 ]

回复此楼

» 猜你喜欢

1楼 2010-03-22 18:42:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)

应助: 0 (幼儿园)
金币: 6963.1
散金: 100
红花: 2
帖子: 2433
在线: 536.4小时
虫号: 932458
注册: 2009-12-23
专业: 凝聚态物性 II ：电子结构

非天之翼蓝羽(金币+2):感谢讨论。。 2010-03-26 14:59

有一点不明白
natom 5
xred
怎么有8个原子坐标

#Definition of the unit cell
acell 7.3775 7.3775 7.8570 #relax it
rprim
1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.0

#Definition of the atoms
natom 5
typat 1 2 3 3 3
xred
0.0 0.0 0.0
0.5 0.5 0.5
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5

[ Last edited by lono75 on 2010-3-26 at 21:49 ]