24小时热门版块排行榜    

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

木虫 (著名写手)

谢谢,不知道这些差别会不会影响我得总能计算?因为感觉后者的原子受力比较大,好想不是最稳定的结构?
7楼2010-11-06 10:06:43
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 289求调剂 +6 步川酷紫123 2026-03-11 6/300 2026-03-17 10:23 by Sammy2
[考研] 0703化学336分求调剂 +3 zbzihdhd 2026-03-15 4/200 2026-03-16 23:52 by zbzihdhd
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[教师之家] 焦虑 +7 水冰月月野兔 2026-03-13 9/450 2026-03-16 10:00 by Quakerbird
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
请填处理意见