| 查看: 6986 | 回复: 17 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
用castep进行结构优化的遇到的一些问题
|
|||
|
各位虫友： 本人最近在做有关层状硅酸盐的第一性原理计算，应用的软件是MS中的castep，由于接触第一性原理时间不长，在结构优化过程中遇到了一些问题，希望得到同行前辈的指点，小弟先谢过了！ 1、根据我看到的相关资料，在利用castep进行计算时，需要对参数（主要是截断能和K点）进行收敛性测试。在论坛中一个帖子上，有同行提到“将体系总能量和晶格常数组为几何优化的目标物理量，但能量截断值大于某一能量值或k点取样达到一定密度时，总能量和晶格常数将趋于一稳定值”，请问各位虫友，收敛性测试具体如何进行？比如，通过什么具体的方法得到晶格常数对K点采样的稳定值。 2、在一篇相关的文献中看到：“The pseudopotentials we used for Si, Al, and O were generated with the core electrons treated using the same Perdew -Wang 91 GGA as for the valence electrons in contrast to the common practice of generating pseudopotentials using the LDA then using the GGA approximation for valence electrons。” 此论文也用的是castep方法，请问各位虫友，在castep中如何选择各元素的赝势，是不是如下图，在potentials选项卡中选择，那么同一元素不同的赝势有什么区别，以及如何选择芯电子和价电子的赝势？ 3、我所要计算的模型叶腊石的分子式为[Al4Si8O20(OH)4]，请问此结构在做结构优化时是否需要考虑自旋极化，如果需要考虑的话，是否需要用LDA+U的方法以及initial spin如何设定？或者能够推荐相关的资料或者书籍，谢谢。 4、结构优化后的结果和实验结果对比的话，根据大家的经验，一般模拟的晶格参数与实验结果相差多少时就认为模拟的结果是合理的？  potentials.jpg |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
面上意见出来了
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
yangyangwlhx
木虫 (正式写手)
少校
- 应助: 4 (幼儿园)
- 金币: 2953.3
- 散金: 458
- 红花: 5
- 帖子: 710
- 在线: 419.3小时
- 虫号: 2001586
- 注册: 2012-09-14
- 专业: 电化学
|
楼主 您好 请问您在计算时 加U是从计算的什么阶段开始加的 是收敛性测试 还是几何结构优化 还是后面计算其他性质时加的 我知道的有3种方法 不知道合适不合适 求指导 一种是 收敛性测试时不加U 然后几何结构优化时 也不加U 在几何结构优化后 ，计算能带和态密度以及其他性质时，加U ；另外一种方法是 收敛性测试时 不加U 几何结构优化时 加U 此时不计算任何性质 只计算单点能， 在几何结构优化结束后，然后计算其他性质时，继续加U，获得计算结果；最后一种方法是 收敛性测试时不加U，结合结构优化时加U并且同时选择计算能带结构和态密度以及其他性质。 我做的体系中有过渡金属，所以向楼主请教 |
15楼2013-11-30 20:31:53
【答案】应助回帖
★
感谢参与，应助指数 +1
WDD880227: 金币+1, 多谢交流 2012-11-07 19:03:47
感谢参与，应助指数 +1
WDD880227: 金币+1, 多谢交流 2012-11-07 19:03:47
|
我也是刚刚接触，说下我的理解 1、收敛测试是根据一系列几何优化，做出能量曲线，固定K点，测试不同截断能下的能量，然后再固定截断能，测试不同K点的能量。这样做出曲线，找能量低的点或相对稳定的点，不一定是能量最低的参数。 2、赝势的选取不是固定的，你可以对同一体系进行几个函数的测试，然后选取与实验值相对接近的，也可根据文献中的做，看看别人都用什么，差别不是很大，主要是因为近似方法不同。另外，同一个体系用的是相同的函数，不是对每个原子进行函数设置。 3、自旋极化是针对有磁性的体系，如果体系没有磁性不必选择。+U是对结果进行修正的，当理论结果与实验结果相差较大时才考虑+U，而且U值也不是随便加的。 4、晶格尽量接进吧，这个没什么标准，相差越小，说服力越强。 只是个人愚见  |
» 本帖已获得的红花（最新10朵）
2楼2012-11-07 13:27:41
wangliangzhe
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2463.6
- 散金: 5
- 帖子: 598
- 在线: 252.9小时
- 虫号: 302656
- 注册: 2006-12-02
- 专业: 区域环境化学
3楼2012-11-07 14:39:28
送鲜花一朵 |
谢谢您的指点，可我还有一些困惑： 1、有些资料上说收敛性测试是改变截断能或K点对系统进行单点能计算，你说是对系统进行结构优化，到底哪种方法是正确的？而且，我做过比如固定截断能，改变K点对同一系统进行结构优化，如果总能趋于稳定，但晶格常数却变化很大，那该如何处理？就是说通过什么途径的到晶格常数的稳定值。 2、上图给出的选项卡中，应该就是对每个元素可以进行赝势的选择，关键是什么时候需要选择，选择的标准是什么？ 3、您说自旋极化是针对有磁性的体系，请问如何判断体系是否有磁性（不好意思，本人原本不是学化学的，可能问的问题比较低级），或者我给出的这个体系[Al4Si8O20(OH)4]是否具有磁性？ 4、晶格参数和实验越接近就越好的话，除了调节截断能、K点以及赝势的种类，还可以通过调节哪些参数提高计算的准确度，以及选择这些参数有什么标准？ 请各位虫友多多指教！！！ |
4楼2012-11-08 09:22:42
投票:
30