24小时热门版块排行榜    

查看: 2032  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

奋力拼搏

至尊木虫 (文坛精英)

[求助] 如何查看计算精度?

用CASTEP计算单点能时,如何查看能量的计算精度?十分感谢,请赐教!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

天道酬勤
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2013-07-09 11:51:29
奋力拼搏: 金币+10, 有帮助, 计算完成后的能量精度应该在哪里查看呢?十分感谢!! 2013-07-10 17:06:56
个人认为是 calculation-electronic-more-SCF-SCF tolerance中标示的数值。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2013-07-08 17:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖


csgt0: 金币+1, 谢谢 2013-07-11 10:59:43
引用回帖:
2楼: Originally posted by 月只蓝 at 2013-07-08 17:21:42
个人认为是 calculation-electronic-more-SCF-SCF tolerance中标示的数值。

calculation-electronic-more-SCF-SCF tolerance设置的数值,是计算的收敛指标,达到这个数值,计算会停止。
而具体的残差数值 在status.txt中会实时显示
在结果文档.castep中,也能找到。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2013-07-10 17:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
fegg7502: 金币+5, 应助指数+1, 3ks 2013-07-17 08:23:13
引用回帖:
5楼: Originally posted by 奋力拼搏 at 2013-07-11 11:43:02
status.txt文件里面显示的计算精度并没有达到calculation-electronic-more-SCF-SCF tolerance设置的计算精度,但是却显示succesfully-completed,这是为什么呢?如何解决呢?非常渴望您的赐教!!...

其一,显示同步问题
status显示的情况会有延迟的,要看.castep中的记录;做计算steps很长的计算的时候,常常发现,比如说,提示succesfully-completed的时候,status中显示的是做了115步了,残差还没到设定值,但是你打开.castep就会发现,实际已经做到118步,四个残差标准全部达到了,而且甚至会发现,status中的显示残差,不是对应于117步的,而是第97步的,见附图1.

其二,关于计算精度
.casrep文档有最详细的信息,这些信息和你在CASTEP Calculation中设置的是完全一致的。如附图2。达到该精度,计算会收敛,停止。
当然,这个时候只能说,精度已经达到设置数值,具体残差,要看到最后一步的详细信息,以上述计算118步的.castep文档为例,如附图3。
如何查看计算精度?
附图1.jpg


如何查看计算精度?-1
附图2.jpg


如何查看计算精度?-2
附图3.jpg

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 附图1.jpg
  • 2013-07-16 15:29:29, 131.21 K
  • 附件 2 : 附图2.jpg
  • 2013-07-16 15:29:30, 168.94 K
  • 附件 3 : 附图3.jpg
  • 2013-07-16 15:29:31, 61.6 K
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
9楼2013-07-16 15:29:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖


fegg7502: 金币+1, 3ks 2013-07-17 08:23:26
引用回帖:
7楼: Originally posted by 奋力拼搏 at 2013-07-16 14:38:25
在status.txt文件和结果文档.castep文件中具体什么地方?请举例说明,十分感谢!!!...

9楼中描述“status中的显示残差,不是对应于117步的,而是第97步的,见附图1.” 应当改为:
“status中的显示残差,不是对应于115步的或者118步,而是第117步的,见附图1。”
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
10楼2013-07-16 15:33:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖


fegg7502: 金币+1, 鼓励交流 2013-07-17 08:24:04
引用回帖:
4楼: Originally posted by 奋力拼搏 at 2013-07-11 11:39:16
在status.txt和.castep里面的什么地方可以查到计算后的收敛精度?为什么计算单点能的输出文件中没有status.txt文件呢?请您多多指教,十分感谢!!...

我认为是,对于那些计算很快的过程,如你说单点能计算,status文件可能会不显示,等你计算做得更深入了,你还会发现,对于迭代次数很少的计算,可能会显示status,但是convergence.xcd文件可能不显示。
这些文件(status和.xcd文件),都是用来描述计算过程的,我们通常关心计算结果,所以对于这些问题,不必过于纠结。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
11楼2013-07-16 15:43:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
奋力拼搏: 金币+45, ★★★很有帮助, 十分感谢!! 2013-07-16 20:40:22
fegg7502: 金币+1, 3ks 2013-07-17 08:24:21
引用回帖:
7楼: Originally posted by 奋力拼搏 at 2013-07-16 14:38:25
在status.txt文件和结果文档.castep文件中具体什么地方?请举例说明,十分感谢!!!...

从图中可以看到,计算完成,用了13步,status中文件显示为11步,能量值和残差对于的都是第11步的。对于该单点能计算,体系的能量值和残差应当按第13步所示来取。
如何查看计算精度?-3
附图4.jpg

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 附图4.jpg
  • 2013-07-16 15:52:11, 135.13 K
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
12楼2013-07-16 15:52:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
奋力拼搏: 金币+10, ★★★很有帮助, 十分感谢 2013-07-17 17:13:10
csgt0: 金币+1, 谢谢 2013-07-17 17:59:37
引用回帖:
14楼: Originally posted by 奋力拼搏 at 2013-07-17 17:04:30
您好!请问附图4中Energy gain per atom的值为负值是怎么回事?求赐教!!...

本次计算值与上一次计算值的差。本次计算值若比上一次的大,差值就为正;反之,则为负。判断计算是否收敛,是看此差值的绝对值。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
15楼2013-07-17 17:09:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
16楼: Originally posted by 奋力拼搏 at 2013-07-17 17:13:32
请问您做过磁性合金吗?...

我没有做过磁性材料的计算,建议你去第一性原理版块看看
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
17楼2013-07-17 17:29:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 奋力拼搏 的主题更新
信息提示
请填处理意见