24小时热门版块排行榜    

查看: 1514  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

njuxyh

木虫 (著名写手)

[交流] 【求助】SIESTA计算总能等自洽计算时力准不准阿 已有3人参与

siesta: Atomic forces (eV/Ang):
     1    0.002280    0.001885    0.000011
     2   -0.000823   -0.002001    0.000006
     3   -0.002785    0.000413   -0.000004
     4   -0.003841   -0.003846    0.000010
     5   -0.003414    0.002833   -0.000021
     6    0.002270    0.000616    0.000007
     7   -0.000137    0.002821    0.000011
     8    0.000639    0.000494   -0.000006
     9   -0.004456    0.004070    0.000019
上面是优化后的受力,已经达收敛标准了
我用这个结构,再做scf计算
发现siesta: Atomic forces (eV/Ang):
   1    0.002875   -0.084867    0.035258
   2   -0.000382    0.039908    0.042227
   3   -0.003277    0.002771    0.036465
   4   -0.001384   -0.004595    0.035799
   5   -0.001858    0.005840    0.035947
   6    0.006556   -0.000592    0.036087
   7    0.001250   -0.038206    0.042518
   8    0.000911    0.087377    0.035318
   9   -0.003271   -0.001110    0.034860
  10    0.002253   -0.013611    0.034873
  11   -0.001831    0.004953    0.035918
  12    0.001754   -0.005937    0.035983
可以看出,自洽计算时的输出的力和优化后输出的力不一样阿
现在的问题:是不是 SIESTA计算总能等自洽计算时力算不准呢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

njuxyh

木虫 (著名写手)

这就是自洽最后的结果:
1    0.002875   -0.084867    0.035258
   2   -0.000382    0.039908    0.042227
   3   -0.003277    0.002771    0.036465
   4   -0.001384   -0.004595    0.035799
   5   -0.001858    0.005840    0.035947
   6    0.006556   -0.000592    0.036087
   7    0.001250   -0.038206    0.042518
   8    0.000911    0.087377    0.035318
   9   -0.003271   -0.001110    0.034860
  10    0.002253   -0.013611    0.034873
  11   -0.001831    0.004953    0.035918
  12    0.001754   -0.005937    0.035983

本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2565072&page=1#pid1801284

[ Last edited by njuxyh on 2010-11-5 at 17:56 ]
3楼2010-11-05 17:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

zhangguangping

木虫 (著名写手)


zxzj05(金币+1):感谢交流! 2010-11-05 21:04:27
引用回帖:
Originally posted by njuxyh at 2010-11-05 08:13:37:
siesta: Atomic forces (eV/Ang):
     1    0.002280    0.001885    0.000011
     2   -0.000823   -0.002001    0.000006
     3   -0.002785    0.000413   -0.000004
     4   -0.003841   -0.003846   ...

如果你新的计算SCF自洽结束的话,应该差别不是很大。
弘德明志博学笃行
2楼2010-11-05 17:18:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)


zxzj05(金币+1):感谢交流! 2010-11-05 21:04:44
njuxyh(金币+5): 2010-11-08 10:12:43
引用回帖:
Originally posted by njuxyh at 2010-11-05 10:55:28:
这就是自洽最后的结果:
1    0.002875   -0.084867    0.035258
   2   -0.000382    0.039908    0.042227
   3   -0.003277    0.002771    0.036465
   4   -0.001384   -0.004595    0.035799
   5   - ...

好像这个就是传说中的不同的初始DM猜测得到不同的结果吧。我在siesta-3.0-b和siesta-3.0-rc2之间遇到过这个问题。
弘德明志博学笃行
4楼2010-11-05 19:46:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

njuxyh

木虫 (著名写手)

你是说版本不同造成的?
5楼2010-11-05 20:19:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 7/350 2026-03-17 19:09 by 捣蛋猪猪
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +3 小李0854 2026-03-16 3/150 2026-03-17 13:40 by 热情沙漠
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
请填处理意见