24小时热门版块排行榜    

查看: 1534  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
请填处理意见