24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2094  |  回复: 19

sunyang1988

至尊木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1):感谢回复 2010-09-15 15:58:15
引用回帖:
Originally posted by vasp001 at 2010-09-15 06:52:18:


似乎你的说法有问题啊。
应该是能量收敛,然后才会计算力是否收敛吧
根据 赫尔曼-费曼定理,应该是这个样子。

我驰豫结果一般是能量收敛比力要快
心之所向无惧无悔
11楼2010-09-15 14:08:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)

引用回帖:
Originally posted by lgxyz at 2010-09-15 09:02:03:
多谢各位!如果采用力作为收敛标准,一般要多大,有什么依据呢?

这个不知道,我一般取-0.01
心之所向无惧无悔
12楼2010-09-15 14:09:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

也有人说取-0.02, -0.05的也有,真是不知道该怎么取。
13楼2010-09-15 15:22:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-09-16 15:03:56
引用回帖:
Originally posted by lgxyz at 2010-09-15 15:22:00:
也有人说取-0.02, -0.05的也有,真是不知道该怎么取。

自己试一试,够用就行了
计算资源好的话标准可以高一点啊
心之所向无惧无悔
14楼2010-09-15 16:13:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peng.029

捐助贵宾 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
EDIFFG是离子弛豫收敛,如果默认的话是前后两步能量之差达到EDIFF*10就停止了,如果是力的收敛一般取负值表示,一般参考文献取值的,有对比意义。就我的情况,因为做NEB计算,对力的要求很高,而收敛很难,所以取-0.01时在NSW=300的情况下都没有收敛,所以就考虑提高到-0.03先,这个值是没有绝对的取法的。
上善若水。
15楼2010-09-15 17:40:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

多谢!既然NEB对力要求比较高,而且很难收敛,为何不该用能量?
16楼2010-09-15 18:38:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
229747楼: Originally posted by sunyang1988 at 2010-09-13 16:56:34
OUTCAR用vim打开命令模式下输入?drift可以看最后一个离子步的drift,有每一个原子的受力情况和总和total drift

那看力是否收敛,是是说离子要在xyz三个方向上都得收敛呢?

POSITION                                       TOTAL-FORCE (eV/Angst)
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      0.00000      0.00000      2.68333         0.000000      0.000000      0.016551
      1.57643      2.27340      1.35837         0.000000      0.000000     -0.016551
      1.57643      1.36687      3.43680         0.000000      0.021828     -0.006638
      1.57643      3.17994      3.43680         0.000000     -0.021828     -0.006638
      0.00000      3.64028      0.60490         0.000000      0.021828      0.006638
      0.00000      0.90653      0.60490         0.000000     -0.021828      0.006638
-----------------------------------------------------------------------------------
    total drift:                                0.000000      0.000000      0.000000

从我给出的这个数据看出离子收敛的力标准是多少呢?
臾兮,臾兮,何时是归期!
17楼2012-08-06 10:21:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2405415楼: Originally posted by kmw.8668 at 2012-08-06 10:21:20
那看力是否收敛,是是说离子要在xyz三个方向上都得收敛呢?

POSITION                                       TOTAL-FORCE (eV/Angst)
POSITION                                       TOTAL-FORCE (eV/Ang ...

这个是三个方向的分力,收敛条件要总的里收敛
心之所向无惧无悔
18楼2012-08-07 22:21:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2405580楼: Originally posted by sunyang1988 at 2012-08-07 22:21:13
这个是三个方向的分力,收敛条件要总的里收敛...

那你的意思就是说在OUTCAR文件里看到是每个原子在三个方向上受到的分力大小,而每个原子总的受力没有显示出来。那total drift对应的三个数值指的是所有原子分别在xyz方向受到的力之和还是?total drift 是 0 0 0说明了什么?
臾兮,臾兮,何时是归期!
19楼2012-08-08 14:57:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2405630楼: Originally posted by kmw.8668 at 2012-08-08 14:57:32
那你的意思就是说在OUTCAR文件里看到是每个原子在三个方向上受到的分力大小,而每个原子总的受力没有显示出来。那total drift对应的三个数值指的是所有原子分别在xyz方向受到的力之和还是?total drift 是 0 0 0说 ...

total drift是晶胞内所有原子受力之和,也是分三个方向。0 0 0说明整体不受力
心之所向无惧无悔
20楼2012-08-08 19:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgxyz 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 12:03 by onwj4wpxp2
[基金申请] 精华III评审感受-评审感受-评审感受 +14 ferrarichen 2026-05-11 18/900 2026-05-15 11:12 by cmhchen
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +9 dd921ww 2026-05-12 10/500 2026-05-15 10:41 by muyiliuhui
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 09:23 by onwj4wpxp2
[基金申请] 青C资助名额大幅增加! +10 西葫芦炒鸡蛋 2026-05-13 14/700 2026-05-15 09:07 by gy116024
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 cjf4bx70cj 2026-05-14 4/200 2026-05-15 09:03 by gagyerk94e
[文学芳草园] 风把牡丹吹跑了 +4 myrtle 2026-05-12 7/350 2026-05-14 23:58 by myrtle
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 5/250 2026-05-14 20:15 by 一切都是空工
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[考博] 申博自荐 +4 食品的橙子 2026-05-09 6/300 2026-05-14 16:05 by great1919
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[基金申请] 请问大佬b0816评完了吗 +3 市民华南虎 2026-05-12 7/350 2026-05-14 07:41 by 市民华南虎
[论文投稿] 有带发论文的吗 +3 山楂之术 2026-05-09 3/150 2026-05-13 17:56 by Cyhcl2629
[硕博家园] 导师各种操作恶心咋办 +11 苍白的小青天 2026-05-09 13/650 2026-05-13 17:11 by 六两废铜
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
[考博] 西南大学考核制博士 +3 lijunjie84 2026-05-11 6/300 2026-05-12 18:09 by lijunjie84
[文学芳草园] 窗边初夏的小雨 +7 阿美_Lml888 2026-05-09 10/500 2026-05-12 15:27 by 阿美_Lml888
[考博] 现在不知道怎么办,感觉很痛苦 +4 qweww 2026-05-11 5/250 2026-05-11 20:23 by Oversize
信息提示
请填处理意见