【求助】CASTEP模块几何优化失败，求解决方案。失败原因详见本帖内容 - 第一性原理 - MS

版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

chdx0618

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2362.6
散金: 796
红花: 6
帖子: 229
在线: 82.9小时
虫号: 1196285
注册: 2011-01-26
性别: GG
专业: 无机非金属类光电信息与功

引用回帖:

Originally posted by acridine at 2011-03-24 20:40:07:
楼主莫慌

降低精度之后的这个错误信息出现在什么位置？是一开始计算就出现了，还是在计算了几个能量之后出现的，亦或是在完成了一次scf收敛之后出现的？我猜是第三种

能否把错误出现之前的一 ...

你好！看来就是你说的第三种情分况，后半部分文件如下：


+---------------- MEMORY AND SCRATCH DISK ESTIMATES PER NODE -----------------+
|                                                    Memory       Disk |
| Model and support data                            404.4 MB    409.0 MB |
| Electronic energy minimisation requirements       426.3 MB    185.5 MB |
|                                              ----------------------------- |
| Approx. total storage required per node          830.7 MB    594.5 MB |
|                                                                            |
| Requirements will fluctuate during execution and may exceed these estimates |
+-----------------------------------------------------------------------------+
Calculating finite basis set correction with  3 cut-off energies.
Calculating total energy with cut-off of  340.000eV.
------------------------------------------------------------------------ <-- SCF
SCF loop    Energy          Fermi          Energy gain    Timer <-- SCF
                           energy       per atom       (sec) <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -2.43247486E+004  4.83713644E+001                      108.58  <-- SCF
   1  -3.25071705E+004  3.23760606E+000 1.36373699E+002    481.31  <-- SCF
   2  -3.35041494E+004  3.30066966E+000 1.66163142E+001    837.30  <-- SCF
   3  -3.35353740E+004  3.21807494E+000 5.20411473E-001 1129.96  <-- SCF
   4  -3.35170771E+004  1.94314821E+000  -3.04949337E-001 1448.58  <-- SCF
   5  -3.34999662E+004  2.50333578E+000  -2.85182193E-001 1798.88  <-- SCF
   6  -3.34907770E+004  2.51791447E+000  -1.53152541E-001 2164.25  <-- SCF
   7  -3.35587570E+004  1.56213346E+000 1.13300021E+000 2481.02  <-- SCF
   8  -3.34954701E+004  2.89507610E+000  -1.05478204E+000 2822.10  <-- SCF
   9  -3.34838323E+004  2.74292404E+000  -1.93963688E-001 3175.07  <-- SCF
   10  -3.34811649E+004  2.49483488E+000  -4.44560478E-002 3499.60  <-- SCF
   11  -3.34820324E+004  2.45781979E+000 1.44581381E-002 3815.71  <-- SCF
   12  -3.34824619E+004  2.50274836E+000 7.15895446E-003 4134.60  <-- SCF
   13  -3.34829532E+004  2.47632455E+000 8.18691249E-003 4452.47  <-- SCF
   14  -3.34829828E+004  2.46503714E+000 4.94533643E-004 4771.07  <-- SCF
------------------------------------------------------------------------ <-- SCF

2*Integrated Spin Density =    41.2082
2*Integrated |Spin Density| =    44.8577

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

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

Calculating total energy with cut-off of  345.000eV.
------------------------------------------------------------------------ <-- SCF
SCF loop    Energy          Fermi          Energy gain    Timer <-- SCF
                           energy       per atom       (sec) <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -3.34773315E+004  2.58332620E+000                   4949.79  <-- SCF
   1  -3.35621441E+004  2.57283318E+000 1.41354302E+000 5252.58  <-- SCF
   2  -3.35622668E+004  2.57253163E+000 2.04620281E-003 5609.02  <-- SCF
   3  -3.38328350E+004  1.50218567E+000 4.50947001E+000 5971.61  <-- SCF
   4  -3.38361142E+004  1.48411522E+000 5.46522223E-002 6295.31  <-- SCF
   5  -3.35750493E+004  2.88090445E+000  -4.35108115E+000 6658.47  <-- SCF
   6  -3.35758944E+004  2.87679861E+000 1.40850889E-002 6945.71  <-- SCF
   7  -3.35073800E+004  2.74869630E+000  -1.14190711E+000 7301.81  <-- SCF
   8  -3.34909281E+004  2.83067048E+000  -2.74198151E-001 7683.81  <-- SCF
   9  -3.34832807E+004  2.37640441E+000  -1.27456432E-001 8042.82  <-- SCF
   10  -3.34782473E+004  2.38656033E+000  -8.38894192E-002 8371.45  <-- SCF
   11  -3.34810172E+004  2.72760682E+000 4.61641968E-002 8696.10  <-- SCF
   12  -3.34780946E+004  2.61682826E+000  -4.87104438E-002 9018.32  <-- SCF
   13  -3.34775561E+004  2.67276882E+000  -8.97462792E-003 9356.19  <-- SCF
   14  -3.34776280E+004  2.68039011E+000 1.19930141E-003 9676.66  <-- SCF
------------------------------------------------------------------------ <-- SCF

2*Integrated Spin Density =    22.2209
2*Integrated |Spin Density| =    27.5608

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

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


Error in calculate_finite_basis:
   previous calculation had total energy = -33482.98282286    eV
   and this calculation has total energy = -33477.62803476    eV
   => something has gone seriously wrong with minimiser

Checkpoint file cannot be written.
Error in calculate_finite_basis: total energies are not variational - aborting
Current trace stack:
calculate_finite_basis_corr
check_elec_ground_state
castep
万望不吝赐教！

[ Last edited by chdx0618 on 2011-3-25 at 15:00 ]

赞一下

回复此楼

21楼2011-03-25 14:59:32

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

galiha

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1581.1
帖子: 135
在线: 163.1小时
虫号: 1195405
注册: 2011-01-24
专业: 极端条件下使用的金属材料

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

我也遇到了此类问题~同求解

赞一下(1人)

回复此楼

22楼2011-03-25 15:11:45

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

acridine

木虫 (正式写手)

1ST强帖: 11
应助: 13 (小学生)
金币: 4901.6
红花: 13
帖子: 706
在线: 389.5小时
虫号: 305478
注册: 2006-12-10
性别: GG
专业: 凝聚态物性 II ：电子结构

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gzqdyouxia(金币+2): 谢谢提示 2011-03-25 18:33:18

现在的问题应该是：虽然收敛了，但是由于精度低，在进行finite_basis修正的时候Ecut的改变导致能量变化过大，软件判定出现了问题 (当然有问题了: 精度过低)。不过好消息是收敛了，这也是当初选择降低精度的目的。

你现在可以考虑不管finite_basis修正而仅仅考虑低精度的scf计算，利用低精度得到的电子密度作为初试值再进行高精度的计算（castep应该有这个功能吧，你可以在help查一下restart calculation等关键词），这样可能有助于高精度计算时候的收敛。

赞一下(5人)

回复此楼

23楼2011-03-25 16:45:08

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

chdx0618

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2362.6
散金: 796
红花: 6
帖子: 229
在线: 82.9小时
虫号: 1196285
注册: 2011-01-26
性别: GG
专业: 无机非金属类光电信息与功

引用回帖:

Originally posted by acridine at 2011-03-25 16:45:08:
现在的问题应该是：虽然收敛了，但是由于精度低，在进行finite_basis修正的时候Ecut的改变导致能量变化过大，软件判定出现了问题 (当然有问题了: 精度过低)。不过好消息是收敛了，这也是当初选择降低精度的目的。 ...

专业、高手！我先试试看！

赞一下

回复此楼

24楼2011-03-26 09:55:09

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

dali2006

木虫 (著名写手)

小木虫酱油店店主

应助: 0 (幼儿园)
金币: 3136
散金: 132
帖子: 1165
在线: 151.7小时
虫号: 295050
注册: 2006-11-11
性别: GG
专业: 材料

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1): 谢谢建议 2011-03-26 12:09:55

内存溢出的话可以考虑用job control里面的MEMORY模式

赞一下(2人)

回复此楼

25楼2011-03-26 10:45:57

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

awmc2008

至尊木虫 (文坛精英)

学习者

应助: 62 (初中生)
金币: 15981.6
红花: 10
帖子: 12168
在线: 732.6小时
虫号: 600677
注册: 2008-09-12
性别: GG
专业: 凝聚态物性 II ：电子结构

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢建议 2011-03-26 12:10:17

SCF＝200，cut-off energies设大些，Empty bands=44, Runtime optimization: memory。这样计算速度会慢一下。建议增大内存，大于14G。知道计算的体系，才可提进一步的建议。

赞一下(2人)

回复此楼

xuexizhe

26楼2011-03-26 11:09:48

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

chdx0618

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2362.6
散金: 796
红花: 6
帖子: 229
在线: 82.9小时
虫号: 1196285
注册: 2011-01-26
性别: GG
专业: 无机非金属类光电信息与功

引用回帖:

Originally posted by awmc2008 at 2011-03-26 11:09:48:
SCF＝200，cut-off energies设大些，Empty bands=44, Runtime optimization: memory。这样计算速度会慢一下。建议增大内存，大于14G。知道计算的体系，才可提进一步的建议。

你好！我计算的体系是Ag掺杂LaBaMnO体系，这个体系的掺杂量大约是0.08333（Ag：LBMO），整个体系共60个原子，其中一个La位被Ag替代。
另外电脑的物理内存已不能再增加，两个插槽，每个4G。望不吝赐教！

[ Last edited by chdx0618 on 2011-3-27 at 00:33 ]

赞一下

回复此楼

27楼2011-03-26 23:38:12

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