24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3934  |  回复: 20

carlaty

金虫 (正式写手)

[交流] 【求助】castep收敛性测试问题 已有8人参与

查看ms帮助文档的时候,对收敛性测试其有如下陈述:

The only parameter that must be known in order for this correction term to be evaluated is dEtot/d lnEcut, where Etot is the total energy of the system and Ecut is the cutoff energy. CASTEP can calculate this term automatically or you can input this parameter manually.

最后一句话提到castep能自动计算该项,请问这一项计算结果在哪查?

先谢过啦
回复此楼

» 收录本帖的淘帖专辑推荐

第一原理资料汇编 晶体模型计算分析 模拟知识汇总 第一性原理
Materials-studio

» 猜你喜欢

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

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

gongjijun

木虫 (正式写手)


lvjian8596(金币+1):感谢交流! 2010-07-01 13:05:22
carlaty(金币+1):谢谢 2010-07-02 09:55:15
dEtot/d lnEcut  这个在它的输出文件.castep中有的,好像一般测试三个,你可以看看输出文件。
2楼2010-07-01 12:55:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carlaty

金虫 (正式写手)

引用回帖:
Originally posted by gongjijun at 2010-07-01 12:55:17:
dEtot/d lnEcut  这个在它的输出文件.castep中有的,好像一般测试三个,你可以看看输出文件。

测试三个?能仔细说说测试哪三个么?
3楼2010-07-02 09:55:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carlaty

金虫 (正式写手)

引用回帖:
Originally posted by gongjijun at 2010-07-01 12:55:17:
dEtot/d lnEcut  这个在它的输出文件.castep中有的,好像一般测试三个,你可以看看输出文件。

我今天又看了一下,貌似dEtot/d lnEcut 会在xcd文件中以曲线的形式给出来。
不过怎么我计算的时候,有时有有时没有xcd文件输出,很诡异。这是怎么回事啊
各位大侠请不吝赐教啊,谢过了

[ Last edited by carlaty on 2010-7-2 at 10:02 ]
4楼2010-07-02 09:57:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanny0007

金虫 (小有名气)

搭车同问,怎么自动测试
Dowhatyousay,saywhatyoudo.
5楼2010-07-02 09:58:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongjijun

木虫 (正式写手)


carlaty(金币+1): 2010-07-02 10:39:44
gavinliu7390(金币+1):谢谢交流! 2010-07-02 11:25:53
比方说,你选fine精度,它会首先在370,375,380三个点计算单点能,然后给出dEtot/d lnEcut 的关系,在。castep文件的前面有,可以看看。
6楼2010-07-02 10:31:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carlaty

金虫 (正式写手)

引用回帖:
Originally posted by gongjijun at 2010-07-02 10:31:51:
比方说,你选fine精度,它会首先在370,375,380三个点计算单点能,然后给出dEtot/d lnEcut 的关系,在。castep文件的前面有,可以看看。

不好意思,能不能给出您的结果参考一下,特别是其在哪个部分(也就是用    ************ #### **************      分割的    #####    的名字),以及该变量的名字,我查了,没找到dEtot/d lnEcut。
十分感谢。
7楼2010-07-02 10:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongjijun

木虫 (正式写手)

★ ★
carlaty(金币+20):十分感谢,追加20个金币以表心意 2010-07-02 11:24:59
gavinliu7390(金币+2):谢谢交流! 2010-07-02 11:26:16
franch: 回帖置顶 2012-03-26 17:16:11
引用回帖:
Originally posted by carlaty at 2010-07-02 10:50:18:

不好意思,能不能给出您的结果参考一下,特别是其在哪个部分(也就是用    ************ #### **************      分割的    #####    的名字),以及该变量的名字,我查了,没找到dEtot/d lnEcut。
十分感谢。

看下面,我把测试改成了九步走,比三步走更好看点。中间有步骤省略了。



                                            0.00000
                                              0.00000
  
+---------------- MEMORY AND SCRATCH DISK ESTIMATES PER NODE -----------------+
|                                                     Memory          Disk    |
| Model and support data                               15.4 MB         1.1 MB |
| Electronic energy minimisation requirements           0.5 MB         0.3 MB |
|                                               ----------------------------- |
| Approx. total storage required per node              16.0 MB         1.4 MB |









|                                                                             |
| Requirements will fluctuate during execution and may exceed these estimates |
+-----------------------------------------------------------------------------+
Calculating finite basis set correction with  9 cut-off energies.
Calculating total energy with cut-off of  200.000eV.
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -3.07170523E+002  3.25956370E+001                         2.20  <-- SCF
      1  -8.49527856E+002 -7.43184326E+000   5.42357333E+002       2.49  <-- SCF
      2  -8.57940741E+002 -8.29247555E+000   8.41288542E+000       2.74  <-- SCF
      3  -8.58117902E+002 -8.31780439E+000   1.77160512E-001       2.95  <-- SCF
      4  -8.57631481E+002 -1.07793770E+001  -4.86420644E-001       3.26  <-- SCF
      5  -8.57469860E+002 -9.83463613E+000  -1.61620409E-001       3.57  <-- SCF
      6  -8.57467551E+002 -9.90059122E+000  -2.30898949E-003       3.87  <-- SCF
      7  -8.57464933E+002 -9.78596737E+000  -2.61839184E-003       4.17  <-- SCF
      8  -8.57464773E+002 -9.81369134E+000  -1.60263953E-004       4.47  <-- SCF
      9  -8.57464768E+002 -9.81706461E+000  -4.90571523E-006       4.67  <-- SCF
     10  -8.57464764E+002 -9.82083636E+000  -3.54776991E-006       4.89  <-- SCF
     11  -8.57464765E+002 -9.82015533E+000   1.27376553E-007       5.11  <-- SCF
------------------------------------------------------------------------ <-- SCF

Final energy, E             =  -857.4433224441     eV
Final free energy (E-TS)    =  -857.4647645076     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -857.4540434759     eV

Calculating total energy with cut-off of  205.000eV.
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -8.57464432E+002 -9.82075585E+000                         5.33  <-- SCF
      1  -8.58942171E+002 -9.99653292E+000   1.47773912E+000       5.53  <-- SCF
      2  -8.58954120E+002 -9.99890596E+000   1.19485237E-002       5.80  <-- SCF
      3  -8.58953500E+002 -9.51439735E+000  -6.19523684E-004       6.10  <-- SCF
      4  -8.58935522E+002 -9.85707485E+000  -1.79782142E-002       6.40  <-- SCF
      5  -8.58933712E+002 -9.80681747E+000  -1.81052872E-003       6.70  <-- SCF
      6  -8.58933786E+002 -9.80301217E+000   7.40725410E-005       6.91  <-- SCF
      7  -8.58933783E+002 -9.80686032E+000  -2.37219392E-006       7.13  <-- SCF
      8  -8.58933783E+002 -9.80610478E+000  -7.16792642E-007       7.35  <-- SCF
------------------------------------------------------------------------ <-- SCF

Final energy, E             =  -858.9244304167     eV
Final free energy (E-TS)    =  -858.9337826323     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -858.9291065245     eV

Calculating total energy with cut-off of  210.000eV.
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -8.58934529E+002 -9.80758472E+000                         7.59  <-- SCF
      1  -8.60008811E+002 -1.00423903E+001   1.07428227E+000       7.78  <-- SCF
      2  -8.60016310E+002 -1.00428506E+001   7.49936771E-003       8.06  <-- SCF
      3  -8.60010134E+002 -9.62353478E+000  -6.17609717E-003       8.36  <-- SCF
      4  -8.60004345E+002 -9.91967268E+000  -5.78893350E-003       8.66  <-- SCF
      5  -8.60001749E+002 -9.82792323E+000  -2.59652186E-003       8.96  <-- SCF
      6  -8.60001787E+002 -9.83460072E+000   3.80898267E-005       9.24  <-- SCF
      7  -8.60001784E+002 -9.83116936E+000  -3.28152156E-006       9.42  <-- SCF
      8  -8.60001783E+002 -9.83133816E+000  -4.52992179E-007       9.63  <-- SCF
------------------------------------------------------------------------ <-- SCF

Final energy, E             =  -859.9936543751     eV
Final free energy (E-TS)    =  -860.0017830920     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -859.9977187335     eV

Calculating total energy with cut-off of  215.000eV.
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -8.60003597E+002 -9.83119667E+000                        10.14  <-- SCF
      1  -8.60749334E+002 -9.95938744E+000   7.45737327E-001      10.63  <-- SCF
      2  -8.60754333E+002 -9.95942120E+000   4.99880438E-003      11.47  <-- SCF
      3  -8.60753090E+002 -9.73747586E+000  -1.24314424E-003      12.39  <-- SCF
      4  -8.60757014E+002 -1.00048195E+001   3.92417161E-003      13.27  <-- SCF
      5  -8.60749959E+002 -9.85398066E+000  -7.05533483E-003      14.09  <-- SCF
      6  -8.60750155E+002 -9.86768033E+000   1.95926753E-004      14.86  <-- SCF
      7  -8.60750147E+002 -9.86199704E+000  -7.82980885E-006      15.46  <-- SCF
      8  -8.60750147E+002 -9.86221004E+000   3.45047829E-008      16.06  <-- SCF
------------------------------------------------------------------------ <-- SCF

Final energy, E             =  -860.7392441661     eV
Final free energy (E-TS)    =  -860.7501470593     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -860.7446956127     eV

Calculating total energy with cut-off of  220.000eV.
------------------------------------------------------------------------ <-- SCF
.
.
.
.


Final energy, E             =  -863.0492643017     eV
Final free energy (E-TS)    =  -863.0635159963     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -863.0563901490     eV

Calculating total energy with cut-off of  240.000eV.
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -8.63067779E+002 -9.79659805E+000                        30.71  <-- SCF
      1  -8.63338558E+002 -9.84424156E+000   2.70779591E-001      31.03  <-- SCF
      2  -8.63339808E+002 -9.84444594E+000   1.25031551E-003      31.48  <-- SCF
      3  -8.63339614E+002 -9.72849119E+000  -1.94767268E-004      32.00  <-- SCF
      4  -8.63343264E+002 -9.91332530E+000   3.64982933E-003      32.47  <-- SCF
      5  -8.63338761E+002 -9.78452358E+000  -4.50274257E-003      32.95  <-- SCF
      6  -8.63338815E+002 -9.79725940E+000   5.46300398E-005      33.40  <-- SCF
      7  -8.63338809E+002 -9.79225127E+000  -6.61629972E-006      33.81  <-- SCF
      8  -8.63338805E+002 -9.79134107E+000  -3.38042967E-006      34.19  <-- SCF
------------------------------------------------------------------------ <-- SCF

Final energy, E             =  -863.3235868454     eV
Final free energy (E-TS)    =  -863.3388053914     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -863.3311961184     eV

For future reference: finite basis dEtot/dlog(Ecut) =     -31.408713eV
Total energy corrected for finite basis set =    -863.391946 eV

+---------------- MEMORY AND SCRATCH DISK ESTIMATES PER NODE -----------------+
|                                                     Memory          Disk    |
| Model and support data                               15.4 MB         1.1 MB |
| Geometry minimisation requirements                    0.9 MB         0.3 MB |
|                                               ----------------------------- |
| Approx. total storage required per node              16.4 MB         1.4 MB |
|                                                                             |
| Requirements will fluctuate during execution and may exceed these estimates |




图是我自己做的哦。




+-----------------------------------------------------------------------------+

[ Last edited by gongjijun on 2010-7-2 at 11:21 ]
8楼2010-07-02 11:02:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carlaty

金虫 (正式写手)

引用回帖:
Originally posted by gongjijun at 2010-07-02 11:02:04:


看下面,我把测试改成了九步走,比三步走更好看点。中间有步骤省略了。



                                            0.00000
                                              0.00000
  
+------ ...

感谢您的耐心解答,特别是这个图很明显直观,多谢了
9楼2010-07-02 11:26:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongjijun

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by carlaty at 2010-07-02 11:26:05:

感谢您的耐心解答,特别是这个图很明显直观,多谢了

您太客气了,共同进步吧。
正好我今天也在研究这个问题。
10楼2010-07-02 11:29:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 carlaty 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +19 MAX怅惘 2026-04-09 21/1050 2026-04-14 11:11 by 不我拉绿卡
[考研] 085801电气专硕272求调剂 +13 电气李 2026-04-13 14/700 2026-04-14 08:50 by 木木mumu~
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[考研] 材料相关专业344求调剂双非工科学校或课题组 +19 hualkop 2026-04-12 20/1000 2026-04-14 07:02 by laoshidan
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 调剂 +27 不逢春 2026-04-07 28/1400 2026-04-13 20:19 by xxy88888
[考研] 273求调剂 +3 白居不易. 2026-04-09 5/250 2026-04-13 19:27 by or7ro
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 53/2650 2026-04-13 14:20 by 张zhihao
[考研] 化工学硕 285求调剂 +29 Wisjxn 2026-04-07 29/1450 2026-04-13 13:47 by floriea
[考研] 电子信息270求调剂 +18 terminal469 2026-04-07 18/900 2026-04-12 16:23 by ajpv风雷
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 275求调剂 +9 1624447980 2026-04-08 10/500 2026-04-11 10:20 by Delta2012
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
信息提示
请填处理意见