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 ]

回复此楼

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

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

非天之翼蓝羽

金虫 (初入文坛)

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

根据大家的意见，更新了输入文件：
natfix 8
  iatfix 1 1 2 3 3 3 3 3
指定了spgroup 55
xred仍然采用了参考文献里的实验值坐标——我所找的晶体结构cif格式可以在ms导出来，用写字板打开就得到晶体的结构数据——但是log里还是出错了：
symatm : WARNING -
  Trouble finding symmetrically equivalent atoms
  Applying inv of symm number  2 to atom number 1  of typat  1
  gives tratom= -7.0850E-01 -1.3000E-01  0.0000E+00.
  This is further away from every atom in crystal than the allowed tolerance.
  The inverse symmetry matrix is -1  0  0
                              0 -1  0
                              0  0  1
  and the nonsymmorphic transl. tnons = 0.0000000 0.0000000 0.0000000
  The nearest coordinate differs by -4.170E-01 -2.600E-01  0.000E+00
  for indsym(nearest atom)= 1

  This indicates that when symatm attempts to find atoms symmetrically
  related to a given atom, the nearest candidate is further away than some
  tolerance.  Should check atomic coordinates and symmetry group input data.

请教为什么用实验值还是会出错呢？原cif文件里，pbam空间群跟xred共存得很好啊。。。