| 查看: 3546 | 回复: 10 | |||
[求助]
收敛测试和结构优化
|
|||
|
使用CAStep的时候，建完模型之后需要进行结构优化，而选择合适的参数又需要收敛测试 是不是收敛测试就是正交不同的参数计算单点能，画趋势图呢？ 收敛测试后确定好参数再进行结构优化，那么结构优化好之后，测试好的参数在计算性能的过程中会不会又不收敛？ 收敛测试和结构优化那个先进行呢？ |
回复此楼» 猜你喜欢
0703化学321分求调剂
已经有10人回复
-
0703化学求调剂
已经有7人回复
-
一志愿厦门大学材料工程专硕354找调剂！！！
已经有5人回复
-
266分，求材料相关专业调剂
已经有5人回复
-
292求调剂
已经有13人回复
-
材料化工340求调剂
已经有3人回复
-
295材料工程专硕求调剂
已经有9人回复
-
085600 材料与化工 329分求调剂
已经有18人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有23人回复
-
08工科，295，接受跨专业调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- siesta结构优化收敛标准
已经有5人回复
- 用VASP优化结构收敛不了，求指教
已经有10人回复
- 优化结构，如何加强收敛
已经有13人回复
- 如何做vasp掺杂
已经有15人回复
- （文献学习）一篇计算HfO2表面 结构与原子性质的文章
已经有26人回复
- 用Dmol优化收敛的结构，用vasp不收敛？
已经有9人回复
- abinit计算晶体生成焓的问题讨论
已经有7人回复
- 求助：vasp结构优化不能收敛
已经有18人回复
- 分子表面吸附 空间结构如何确定
已经有4人回复
- 【求助】K点取得越大越不收敛，什么原因！
已经有20人回复
- 【求助】结构优化不收敛
已经有9人回复
- 【求助】HF方法优化配合物结构无法收敛
已经有3人回复
- 【原创】关于结构弛豫设置及收敛判据选择的一点体会
已经有140人回复
- 【求助】还是收敛测试和几何优化的问题，又来麻烦大家了...
已经有3人回复
- 【求助】收敛测试怎么做？？？
已经有17人回复
- 【求助】收敛性测试的问题
已经有40人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
zzjjungle
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1215.6
- 红花: 1
- 帖子: 78
- 在线: 225.2小时
- 虫号: 1254708
- 注册: 2011-04-04
- 性别: GG
- 专业: 金属功能材料
★ ★ ★ ★
fzx2008: 金币+4, 辛苦了 :) 2013-03-13 14:44:40
fzx2008: 金币+4, 辛苦了 :) 2013-03-13 14:44:40
|
以前看的帖子，不知道对楼主是否有帮助.... 收敛测试当然是必要的，不然你怎么能确定你的参数是合理的呢。 之所以你会觉得“不同的截断能算出来的能量也差不多” 是因为castep自己会根据你所选择的赝势、晶胞尺寸，来提供预设参数。选择fine精度时，这些参数一般是可以满足能量收敛性的。所以在这些参数附近测试，能量值自然变化不大，因为已经收敛了。 但有两个缺点 1.默认参数可能精度不够；尤其是涉及应力的计算，比如力学性质或晶格优化，其需要对应力进行收敛测试，通常其对参数的要求要高于能量收敛测试。 2.默认参数可能精度过高；对于小体系并没有什么影响，但大体系则很明显。因为截断能和k点数目的不同是可以成倍增加或减少计算量的，在大体系中可能会造成几天的计算差别，所以花个几小时找到适合的参数其实是很划算的。 我说的收敛就是“收敛性测试”里的收敛，并非某一次计算的时候scf的收敛；这是两个概念，“收敛性测试”和几何优化是否成功是没有关系的。 举个例子，比如说你用截断能100eV几何优化某个晶胞，得到晶格参数为3.0A，使用200eV优化为3.1A，使用300eV为3.13A，使用400eV为3.14A...... 随着截断能不断增加，得到的晶格参数逐收敛于某一数值，这个过程就是收敛测试。你要根据你的需要，比如你想让精度达到0.01A，那就要使用400eV以上的截断能，如果是精度0.1A就够，用200eV的截断能就可以了。 但上述这些优化过程都是成功的，但并不代表其参数就是收敛的 顺便帖一个自己测试模守恒赝势的结果图～拟合EOS可以使用ELK自带的eos工具，所以整个过程只需要写个脚本后等上一段时间来收数据和画图就足够了（如果使用QE自带的工具，由于ev.x属于人机交互型的程序，需要人工干预，比较缺乏效率）～ 比较好的测试收敛的方法不是直接在一个结构上调整cutoff计算，还是在取比较高的kpt之后，每个cutoff做一个EOS拟合，然后对比EOS拟合后的总能和晶格常数。这是因为对于不同的cutoff，对应的最低能量点也有很小量的偏移~ 首先取一个比较大的k点密度，这样保证k点是足够的。然后保持这个k点密度不变，对cutoff值做测试。具体的测试办法是，对应一个cutoff值，不断的改变晶格常数的值，找到总能量和晶格常数的变化关系，然后利用下面的公式： 拟合出公式里面的V0和E0,这样就得到了这个cutoff值对应的V0和E0了。反复这样测试，最后得到cutoff变化时V0和E0的变化了。从而找到一个收敛的cutoff值。是这样的吧？ 公式中的B0和B0导数可以从程序输出中得到是吗？ 我想做这个分子吸附在表面上，所以我是不是该做表面的收敛测试，而不是分子的收敛测试。做好表面的收敛测试直接把分子加上就可以了吧 截断能我一般都是按照参考文献取的，你找一下类似的体系，看一下别人的设置（赝势的类型也要和你一样），这个取大一些肯定没问题，关键是k点要测试到收敛（以我的经验，用参考文献的值能够确保收敛性，比如文献提到的截断能450，k点5*5*5。我会测试一下截断能400,450,500，此时k点取文献值。确定完截断能后，再试k点，比如4*4*4，5*5*5，6*6*6。你会发现确实是收敛的）。对于你算吸附，先做好表面的截断能、k点测试，把表面的结构弛豫好，再加上要吸附的分子，再弛豫。因为弛豫分子主要是为了得到分子能量，计算吸附能。所以我觉得你没必要用分子去测试截断能。在同一个计算中截断能要保持一致，一般氧的截断能比较高，当你的体系有氧的话，截断能一般是450或500。还有收敛性是看两个计算的能量差值，不是看绝对大小的，楼主没必要为此纠结，很容易解决的。 |
6楼2013-03-12 16:05:44
gy850501
木虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 1802.6
- 散金: 20
- 红花: 4
- 帖子: 268
- 在线: 71.5小时
- 虫号: 1036520
- 注册: 2010-06-06
- 性别: GG
- 专业: 金属材料的微观结构
2楼2013-03-12 11:13:24
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
英文电子书: 金币+50, ★★★很有帮助 2013-03-13 10:13:05
fzx2008: 金币+3, 谢谢指导 2013-03-13 14:44:02
感谢参与，应助指数 +1
英文电子书: 金币+50, ★★★很有帮助 2013-03-13 10:13:05
fzx2008: 金币+3, 谢谢指导 2013-03-13 14:44:02
|
1.一般来说，像平面波截断能动、k点网格等影响计算精度的关键参数，用单点能做测试即可。即，增加Ecut或者k点，体系的总能变化很小。 2.确定后的参数再进行结构优化 3.MS比较简单，一般采取Fine的精度，用程序自己确定的参数计算即可。除一些有关磁性的计算外，一般都能很好的跟实验比较。 |
3楼2013-03-12 12:53:44
飞行鸟
至尊木虫 (文坛精英)
- 应助: 2135 (讲师)
- 金币: 22804.4
- 散金: 5054
- 红花: 377
- 沙发: 6
- 帖子: 16032
- 在线: 1486.8小时
- 虫号: 638457
- 注册: 2008-10-27
- 专业: 宇宙学
4楼2013-03-12 14:09:35
5楼2013-03-12 15:00:56
zhoupeng87
版主 (文学泰斗)
- 应助: 6148 (副教授)
- 贵宾: 8.516
- 金币: 47410.3
- 散金: 58001
- 红花: 566
- 沙发: 2391
- 帖子: 64964
- 在线: 9636.5小时
- 虫号: 962700
- 注册: 2010-03-06
- 性别: GG
- 专业: 无机材料化学
- 管辖: 精细化工
7楼2013-03-12 16:35:52
8楼2013-03-13 10:12:22
9楼2013-03-13 11:53:34
10楼2013-03-13 12:34:14
