24小时热门版块排行榜    

查看: 1502  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
请填处理意见