| 查看: 11407 | 回复: 79 | ||||||||||||||||||||
| 【奖励】 本帖被评价55次，作者fo3mt2增加金币 46.1 个 | ||||||||||||||||||||
[资源]
【讨论】关于针对CASTEP收敛问题的参数调试总结。
|
||||||||||||||||||||
|
大家好： 我根据自己在计算过程中遇到的收敛性的问题，做了一点关于CASTEP计算参数调试的总结。当然，这些只是我的个人心得体会，大家如果觉得有什么问题请提出指正，或是如果有自己独到的见解，请跟帖加上来，供大家一起学习提高。 如果计算在Max. cycle全部用完之后仍然不能收敛，要先查看一下 但多数情况下，我们会发现，能量并没有比刚开始计算时的能量有明显下降，而是围绕着某个平均值在反复地震荡，而且震动幅度也没有减小。那说明是参数设定的问题，此时就算再怎么增加Max. cycle的数值，计算最终还是不会收敛的。 这个时候，最好是增加energy cutoff的数值，降低K point set的数值，增加smearing width数值（当然，如果通过调试前两个数值就可以使计算收敛的话，就最好不要动这个数值）。 再有，如果计算是使用的density mixing的话，那要将mix_charge_amp，mix_spin_amp的数值降低，将mix_energy_cutoff的数值增加到energy cutoff数值的3到4倍，增加mix_history_length的数值。 当然，如果调试上面的参数仍然发现收敛很辛苦的话，就将mixing scheme转成All Band/EDFT,这会解决收敛问题，但是计算时间将会是使用density mixing计算所用时间的3倍以上。 |
回复此楼» 收录本帖的淘贴专辑推荐
 第一性原理计算资源总汇 | 材料计算模拟实用技巧 | MaterialStudio学习使用 | 第一性原理 |
| 第一原理资料汇编 | Materials-studio | 科研软件资料 | 计算MS等 |
| 模拟知识汇总 | first principle | MS资料 | 第一性原理计算 |
| ms and vasp学习 | 科研相关资料 | 第一性原理 | castep问题集合 |
| Ms的经验贴 | MS第一性原理 |
» 本帖已获得的红花（最新10朵）
麓山北客» 猜你喜欢
选择
已经有9人回复
-
三个评委收到的基金包里面的项目是相同的吗？
已经有3人回复
-
讨论贴：好本子一定会被打招呼的淘汰吗
已经有50人回复
-
入职高校3年发表10+SCI，尽人事听天命
已经有42人回复
-
导师不让硕转博，让我去国外读博，能理解吗？
已经有20人回复
-
化学B02口青基 代表作都是什么水平的？向大佬求助
已经有12人回复
-
统计一下：硕士毕业答辩后的谢师宴是学生出钱，还是老师出钱？
已经有9人回复
-
求25博导，金属增材制造方向
已经有3人回复
-
2024杰青和万人领军什么时候会评
已经有3人回复
-
讨厌酒桌文化
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- CASTEP金属优化
已经有7人回复
- 关于周期信号的收敛性的问题
已经有4人回复
- 计算收敛性问题
已经有13人回复
- castep scf不收敛怎么办，怎么办！！ 急急急！！
已经有3人回复
- castep收敛问题
已经有4人回复
- Matalab 程序调试
已经有6人回复
- 封闭腔 连续性收敛问题
已经有5人回复
- 关于收敛性测试问题
已经有6人回复
- Castep计算无法收敛，请高人指教
已经有17人回复
- CASTEP计算问题集
已经有11人回复
- CASTEP计算掺杂不收敛
已经有5人回复
- 第一性原理计算入门操作
已经有5人回复
- GDB单步调试问题
已经有8人回复
- CASTEP关于设置温度参数的计算问题
已经有3人回复
- 个人总结有关MS的各种问题（第一季）
已经有705人回复
- 有关CASTEP Band Structure问题
已经有8人回复
- 关于CASTEP的几何优化的两个疑惑。希望大家一起讨论讨论
已经有6人回复
- CASTEP 能量曲线上升 但是收敛了
已经有13人回复
- CASTEP进行布居分析时遇到的问题
已经有3人回复
- 【求助成功】castep 收敛性测试
已经有22人回复
- 【求助】关于参数估计的问题
已经有4人回复
- 【求助】收敛性测试的问题
已经有40人回复
- 【求助】castep收敛性测试问题
已经有20人回复
- 【求助】初学者咨询安装CASTEP的电脑问题
已经有16人回复
- 【求助】MS 这些收敛参数怎么翻译啊
已经有6人回复
19楼2009-11-22 12:41:11
★★★★★ 五星级,优秀推荐
★ ★ ★
wuli8(金币+3,VIP+0):3ks 11-20 13:00
wuli8(金币+3,VIP+0):3ks 11-20 13:00
|
总结的挺好的，但是指出几点不足，即： 1. 参数调试过于经验化，如果楼主能结合对程序原理的认识， 或许能够把每个参数调节的意义给阐述的更为清晰些； 从优化算法， 体系特性， 物理过程等方面着手来进行分析， 或许可以得到更加全面的认识（如果楼主能够结合上述内容对CASTEP收敛进行详细的分析和总结的话，我会考虑接受楼主的文章，并且发表在小木虫计算与模拟研究RCS期刊上）。 2. 建议楼主在进行平面波DFT优化之前，通过其他的办法， 比如经典力学优化以及DMol中局域波函数进行优化， 这样的话，可以帮助CASTEP在最短的时间内找到合适的最小值， 也就是通常意义上说的加快收敛。 |
4楼2009-11-20 11:13:07
7楼2009-11-20 15:52:44
21楼2009-11-22 22:07:33
| The default density mixing settings use Pulay mixing and conjugate-gradient minimization of each electronic state. You should only attempt to change these parameters if SCF convergence is very poor. Sometimes it helps to reduce the length of the DIIS history from the default value of 20 to a smaller value (5-7). It might also be helpful to decrease the mixing amplitude from the default value of 0.5 to 0.1-0.2. |
50楼2011-04-26 20:01:27
★★★★★ 五星级,优秀推荐
★ ★ ★
wuli8(金币+3,VIP+0):3ks 11-20 13:02
wuli8(金币+3,VIP+0):3ks 11-20 13:02
|
great! 其实在做收敛速度无非是两个目的：１，计算精度要高；２计算速度要快． 一般平面波程序，在计算精度上，主要需要测试ecut(倒空间的最大平面波截断，理论上要求无穷大，因为那样基矢才是完备的．但实际上不需要，也做不到)．还有就是ｋ点取样．理论上要计算无穷个ｋ点，但实际也不可能．我觉得应该先测试ecut,因为你首先要保证每个ｋ点能量是正确的．然后再测试ｋ点．其实到了，最后即使ｋ点取的很大的时候，还是有很大的误差的，这主要来自于系统误差．所以有时候，会造成震荡的情况，这主要还是由于ｅｃｕｔ太小，或者是有的ｋ点根本就难以算准．这个时候，如无特殊需要，基本就认为它收敛了．除非你想要在０．１ｍｅＶ的误差． 但是根据我的想法，我就不认为它能将计算精度控制在０．１ｍｅＶ．即使收敛测试达到了．（当然,这是个人看法！） 第二，就是计算速度．因为目前ＤＦＴ主要还是用自洽的方法，在理论上应该是用新的电荷密度等于旧的电荷密度，但实际上如果这样做，一是收敛很慢，二是容易形成震荡情形。一般的做法是rho_in=mix*tho_old+(1-mix)*rho_old。（具体到每个程序，肯定是做法不同了。）一般认为mix=0.3比较合理。这样收敛的非常快。 所以很多程序在测试各种mix的时候，就是在测试这个。 最后，选择合适的参数。在你所需要的精度范围内，最快的计算。 我们的目的是要解决问题，有时候不需要那么精确。 |
5楼2009-11-20 12:41:51
52楼2011-04-26 20:22:35
3楼2009-11-20 11:09:44
6楼2009-11-20 15:49:31
8楼2009-11-20 15:57:43
9楼2009-11-20 16:15:39
12楼2009-11-21 00:13:38
13楼2009-11-21 09:13:29
14楼2009-11-21 09:26:54
15楼2009-11-21 10:47:32
17楼2009-11-22 10:16:44
18楼2009-11-22 11:16:18
20楼2009-11-22 13:01:03
22楼2009-11-23 01:57:00
23楼2009-11-24 13:08:45
24楼2009-11-26 17:47:45
25楼2009-12-10 10:04:06
26楼2009-12-10 14:50:26
27楼2010-04-22 09:21:05
30楼2010-05-04 22:37:50
31楼2010-05-04 23:56:10
32楼2010-05-05 00:07:01
|
本帖内容被屏蔽 |
33楼2010-08-10 10:52:06
34楼2010-08-10 11:08:35
35楼2010-08-10 18:52:54
36楼2010-08-10 21:08:57
37楼2010-08-11 12:14:31
38楼2010-08-11 21:02:23
39楼2010-12-09 15:04:56
40楼2010-12-09 16:10:23
41楼2010-12-09 16:29:30
42楼2010-12-21 19:52:19
43楼2010-12-21 23:39:25
45楼2011-04-24 15:57:38
46楼2011-04-24 20:08:57
47楼2011-04-25 09:06:02
48楼2011-04-26 10:04:56
49楼2011-04-26 14:33:57
简单回复
wuli82楼
2009-11-20 10:51
回复
zdhlover10楼
2009-11-20 20:21
回复
awmc200811楼
2009-11-20 20:48
回复
linldq16楼
2009-11-21 16:06
回复
顺水28楼
2010-04-22 09:27
回复
56293123529楼
2010-05-04 22:00
回复
hh1985040244楼
2010-12-23 08:54
回复
10