24小时热门版块排行榜    

查看: 1510  |  回复: 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

木虫 (著名写手)

你是说版本不同造成的?
5楼2010-11-05 20:19:37
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703化学调剂 ,六级已过,有科研经历 +7 曦熙兮 2026-03-15 7/350 2026-03-16 16:34 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
信息提示
请填处理意见