| 查看: 2125 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】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计算总能等自洽计算时力算不准呢 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有13人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有26人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助Siesta计算熟手，想学习此软件
已经有8人回复
- siesta 计算电子占据数的问题
已经有3人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- 静态自洽计算遇到的问题
已经有5人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- 【求助】siesta进行ordern计算出错求助
已经有4人回复
- 【求助】我在用ATK做自洽计算的时候，出现如下的错误：
已经有3人回复
- 【求助】VASP中，“结构优化”和“电子自洽计算”的关系？
已经有8人回复
- 【求助】结构优化之后的静态自洽计算
已经有12人回复
- 【讨论】自洽计算与非自洽计算的含义、所做的工作
已经有10人回复
- 【原创】并行自洽计算后得到和节点数一样多的波函数文件wfc
已经有4人回复
- 【求助】静态自洽计算读WAVECAR有用吗？
已经有9人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
- 【求助】态密度计算究竟是自洽计算，还是非自洽计算？
已经有24人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+2):谢谢回复！ 2010-11-08 12:48:01
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+2):谢谢回复！ 2010-11-08 12:48:01
|
即使是同一个版本，在用siesta的时候如果收敛了，再次运算的话，如果用上一个DM，会在1到2个SCF之内自洽。但是如果是做输运用transiesta的话，即使是收敛过的DM，如果再次用，还会发现不收敛，并且经过的SCF比较多。这个确实是一个问题。但是你的问题我在使用siesta的时候没有遇到过。尤其是你的几何结构是按照绝对的收敛条件来的。按照这个结构再去算一次，即使是从头猜测DM，只要收敛了，当然DM的收敛标准一定要和优化的时候用的一样的收敛标准。那么应该收敛后的受力应该不会差别很大，最起码还是满足几何结构优化的时候的那个条件的。你的这个问题可以写信给siesta-l@uam.es，这个是siesta的一个讨论区，但是需要注册才能给这个地址发邮件。有可能你的问题，在这个能找到答案。 |
8楼2010-11-08 10:23:43
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
2楼2010-11-05 17:18:21
|
这就是自洽最后的结果： 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
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
4楼2010-11-05 19:46:34
50