应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » 【求助】RF计算出错，参数设置问题？
111/2返回列表12下一页
查看: 1272  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 非天之翼蓝羽 的 5 个金币

非天之翼蓝羽

金虫 (初入文坛)

[交流] 【求助】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的回帖

非天之翼蓝羽

金虫 (初入文坛)
顶一下，继续求助中，渴盼关注。。。

[ Last edited by 非天之翼蓝羽 on 2010-3-24 at 16:16 ]
赞一下 回复此楼
2楼2010-03-23 10:10:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

非天之翼蓝羽

金虫 (初入文坛)
是问题太简单了么？
赞一下 回复此楼
3楼2010-03-24 16:16:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

非天之翼蓝羽

金虫 (初入文坛)
恳请专家关注……
赞一下 回复此楼
4楼2010-03-26 10:08:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)
非天之翼蓝羽(金币+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 ]
赞一下 回复此楼
由天命,成人事！
5楼2010-03-26 13:28:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

非天之翼蓝羽

金虫 (初入文坛)
引用回帖:
Originally posted by lono75 at 2010-03-26 13:28:06:
有一点不明白
natom 5
xred
怎么有8个原子坐标

我觉得前面出的问题可能来自坐标问题，锆酸铅晶体钙钛矿体心立方结构，晶胞内是有5个原子，但坐标是直接参考文献上的fractional coordinates。沃克夫位置：
Pb1 (4g) 0.6991 0.1228 0.0000
Pb2 (4h)  0.7056 0.1294 0.5000
Zr (8i) 0.2414 0.1248 0.2486
O1 (4g) 0.2756 0.1560 0.0000
O18 (4g) 0.3011 0.0956 0.5000
O2 (8i) 0.0317 0.2622 0.2798
O3 (4f) 0.0000 0.5000 0.2026
O4 (4e) 0.0000 0.0000 0.2293
我也很是疑惑。。唉。固体物理知识太匮乏。。。
赞一下 回复此楼
6楼2010-03-26 14:58:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
非天之翼蓝羽(金币+5):非常感谢~ 2010-03-27 11:13
那么也应该是
natom 8
typat 1 1  2 3 3 3 3 3


增加ntime，nstep还是降低tolvrs跟tolmxf可以同时进行啊
赞一下 回复此楼
储氢家族欢迎储氢研究者！
7楼2010-03-27 01:46:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

非天之翼蓝羽

金虫 (初入文坛)
引用回帖:
Originally posted by zxzj05 at 2010-03-27 01:46:23:
那么也应该是
natom 8
typat 1 1  2 3 3 3 3 3


增加ntime，nstep还是降低tolvrs跟tolmxf可以同时进行啊

就是说文献里的fractional coordinates可以直接用不需要换算哈？非常感谢，我试试~
赞一下 回复此楼
8楼2010-03-27 11:12:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenyun0876

木虫 (正式写手)

老木虫
★
非天之翼蓝羽(金币+8):非常感谢，那么我该如何指定空间群呢？rprim我用angdeg代替了，应该可行吧？谢谢！ 2010-03-27 16:04
qli0218(金币+1):谢谢 2010-04-01 18:54
不指定空间群的话，对acell、xred、rprim的要求很严，误差大一点系统就会误判或无法识别而认为是1.

[ Last edited by wenyun0876 on 2010-3-27 at 14:31 ]
赞一下(1人) 回复此楼
胸无块垒心常泰，腹有诗书气自华！
9楼2010-03-27 14:27:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

非天之翼蓝羽

金虫 (初入文坛)
根据大家的意见，更新了输入文件：
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共存得很好啊。。。
赞一下 回复此楼
10楼2010-03-29 22:18:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 非天之翼蓝羽 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-16 4/200 2026-03-16 20:13 by Wangjingyue
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 304求调剂 +6 小熊joy 2026-03-14 6/300 2026-03-16 12:59 by Iveryant
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
信息提示
关闭
请填处理意见
关闭