【求助】castep做几何优化时其参数如何设置？ - 第一性原理 - MS

11楼2010-11-03 09:27:56

zgt2007

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3224.3
散金: 37
红花: 1
帖子: 549
在线: 306.8小时
虫号: 480229
注册: 2007-12-15
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

Originally posted by 明月李 at 2010-10-28 08:45:17:
k点的选取要进行测试，看能量值的差别小于0.001，选取较小的k点，不可以随便取

在两次几何优化完以后，得到能量差已经在0.001范围内，但是感觉计算的时间太常，这时想让k-point取少些，请问在这种情况下k-point应怎么选取！例如，之前优化选择的是10x10x4，现在想让k-point变少，通常的选取原则是什么？是否三个轴向的k值都要变，还是只变某一轴向或某两个轴向的K值？另外，对actual spacing以多大为宜，请赐教，不胜感激！

自力更生，艰苦奋斗！

12楼2010-11-03 09:51:11

zgt2007

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3224.3
散金: 37
红花: 1
帖子: 549
在线: 306.8小时
虫号: 480229
注册: 2007-12-15
性别: GG
专业: 凝聚态物性I:结构、力学和

★
cenwanglai(金币+1):呵呵，多是这么过来的。每天进步一点，很快就会懂得多了。 2010-11-03 10:23:52

引用回帖:

Originally posted by nodding at 2010-11-03 09:27:56:
“k点的选取要进行测试，看能量值的差别小于0.001，选取较小的k点，不可以随便取”

能量值的差别小于0.001，是指不同k点对应的总能吗？

这个难道不看体系大小么？ 4个原子体系的总能，与 32个原子体系的总 ...

你的问题我回答不了，我也是刚刚接触MS，也是似懂非懂。

赞一下(1人)

自力更生，艰苦奋斗！

13楼2010-11-03 09:53:22

明月李

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 733.5
散金: 9
帖子: 96
在线: 43.9小时
虫号: 850433
注册: 2009-09-17
性别: MM
专业: 理论和计算化学

zgt2007(金币+1): 2010-11-03 11:15:11

我们都只是用VASP上进行自行迭代优化，你是用哪种方式优化参数

14楼2010-11-03 10:12:51

zgt2007

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3224.3
散金: 37
红花: 1
帖子: 549
在线: 306.8小时
虫号: 480229
注册: 2007-12-15
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

Originally posted by 明月李 at 2010-11-03 10:12:51:
我们都只是用VASP上进行自行迭代优化，你是用哪种方式优化参数

我用的是castep！

自力更生，艰苦奋斗！

15楼2010-11-03 10:20:26

snowmami

新虫 (正式写手)

应助: 1 (幼儿园)
金币: 1573.1
散金: 78
红花: 2
帖子: 352
在线: 171.6小时
虫号: 896353
注册: 2009-11-07
专业: 凝聚态物性 II ：电子结构

zgt2007(金币+1): 2010-11-03 16:19:36

引用回帖:

Originally posted by zgt2007 at 2010-10-28 10:23:55:

一般用castep做几何优化时，通过改变k-piont，cutoff energy来看能量变化差别多少，若两次计算优化能量变化很少，正如你所说当两次能量值的差别小于0.001，即认为收敛。但是castep里可调节的选项很多，除了 ...

请问高手：
我现在想优化氧化锌单胞，为什么我算来算去都没有达到两次能量差小于0.001eV啊！k=4*4*2,K=4*4*3,k=4*4*4,k=4*4*5等等我试了好多，都要大于0.001eV,最小的也有0.00152eV,这是为什么呢？
迷惑中，期待解答！

16楼2010-11-03 13:46:11

zgt2007

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3224.3
散金: 37
红花: 1
帖子: 549
在线: 306.8小时
虫号: 480229
注册: 2007-12-15
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

Originally posted by snowmami at 2010-11-03 13:46:11:

请问高手：
我现在想优化氧化锌单胞，为什么我算来算去都没有达到两次能量差小于0.001eV啊！k=4*4*2,K=4*4*3,k=4*4*4,k=4*4*5等等我试了好多，都要大于0.001eV,最小的也有0.00152eV,这是为什么呢？
迷惑中 ...

增加截断能或者K-point试试，一般保持一个不变，改变另一个另一个看看能量差，祝你成功！

自力更生，艰苦奋斗！

17楼2010-11-03 16:22:27

我爱砖头

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 157.5
散金: 3
帖子: 26
在线: 7.7小时
虫号: 3596250
注册: 2014-12-15

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

3楼: Originally posted by 明月李 at 2010-10-28 08:45:17
k点的选取要进行测试，看能量值的差别小于0.001，选取较小的k点，不可以随便取

我想问一下，K点测试是测试什么

赞一下(1人)

18楼2015-04-01 16:59:49

我爱砖头

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 157.5
散金: 3
帖子: 26
在线: 7.7小时
虫号: 3596250
注册: 2014-12-15

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

10楼: Originally posted by zgt2007 at 2010-11-03 09:13:26
每次优化完结构以后就得到一个enthalpy（打开project explorer里面的xxx.castep文件，结果在最下面），两次的差值如果变换很小，就可以做为收敛的标准。还有本人也是刚接触MS，很多也不清楚。
...

我想问一下，如果对之前的模型掺杂，那K-point和Cotoff energy还要做收敛测试吗

赞一下(1人)

19楼2015-04-01 17:07:54

一个坏小孩

禁虫

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

本帖内容被屏蔽

20楼2020-03-18 12:57:25