计算化学公社

标题: MS CASTEP valence electrons 设置怎么一回事？ [打印本页]

作者
Author: jlx_ 时间: 2015-9-3 21:56
标题: MS CASTEP valence electrons 设置怎么一回事？

http://dx.doi.org/10.1016/j.ssc.2012.09.003 这篇文章里提到这么一句话：

The orbitals of Al3s22p1, Cu3d104s1 and Mg2p63s2 are treated as valence electrons

不是很明白，这个怎么设置的呢？

作者
Author: 卡开发发 时间: 2015-9-4 00:41
没法详细设置具体价电子组态，在electronic-more-pseudopotential下面能够看到以及选择一些其他的赝势。

有些赝势文件赝化的电子比较少，计算的时候价电子采用数目比较多包含了次亚层甚至跟多，截断能也比较高，比较精确；有些赝化电子很多，则包含很少的几个价电子电子，截断能也比较低，不精确但计算速度快，但并不总是能够给出合理的结果。

作者
Author: jlx_ 时间: 2015-9-4 08:48

卡开发发发表于 2015-9-4 00:41
没法详细设置具体价电子组态，在electronic-more-pseudopotential下面能够看到以及选择一些其他的赝势。

...

谢谢

作者
Author: jiewei 时间: 2015-9-4 09:53
这个手册中有明确说明是可以调的，你看看手册（自带的help文件）

作者
Author: zhongyuabc 时间: 2016-6-7 09:19
请问楼主您问的这个问题知道是怎么回事了吗？

作者
Author: jlx_ 时间: 2016-6-17 22:15

zhongyuabc 发表于 2016-6-7 09:19
请问楼主您问的这个问题知道是怎么回事了吗？

云里雾里，其实@卡卡开发说的就很有道理，正在考虑不用MS，换别的软件

作者
Author: 卡开发发 时间: 2016-6-17 23:31

jlx_ 发表于 2016-6-17 22:15
云里雾里，其实@卡卡开发说的就很有道理，正在考虑不用MS，换别的软件

这个问题我可能没说全，直白点说，一个赝势多少电子赝化多少电子维持active，这些在生成赝势文件时是可调的。一旦赝势文件生成后，这个参数不可再调。也就是说，如果你希望与程序自带的赝势的active的电子数不相同，只能通过生成赝势的程序重新生成新的赝势。

这点，无论是vasp、quantum espresso还是castep等赝势方法均为如此。生成一个好的赝势，需要好的基础功底和足够的经验，需要一定程度降低计算量并且能够给出与全电子计算接近的性质，成熟的商业程序卖点其中之一（如vasp的paw和castep的otfg）就是如此。

Castep其实还算不错，贵了点，而且因为程序的方法比较保守导致计算速度偏慢些。

作者
Author: zhongyuabc 时间: 2016-6-19 09:14
哦，明白了。那我暂时还是用程序默认的吧，因为目前还没有能力去辨别该用什么样的赝势

作者
Author: 卡开发发 时间: 2016-6-19 10:55

zhongyuabc 发表于 2016-6-19 09:14
哦，明白了。那我暂时还是用程序默认的吧，因为目前还没有能力去辨别该用什么样的赝势

你可以在potential选项卡view这些赝势的内容，一般在开头几行会给出一些赝势产生方式相关的信息，比如O有soft和hard的版本，soft里面就写了尽可能“不要用于分子体系计算”这样的提示。

作者
Author: zhongyuabc 时间: 2016-6-19 11:20
嗯，我去看了一下，看到啦！但是发现里面有好多种赝势，不知道不同的赝势是否是因为交换相关泛函不一样？

作者
Author: 卡开发发 时间: 2016-6-19 14:08

zhongyuabc 发表于 2016-6-19 11:20
嗯，我去看了一下，看到啦！但是发现里面有好多种赝势，不知道不同的赝势是否是因为交换相关泛函不一样？

赝势本身分成模守恒和超软两种；按照泛函分为LDA和GGA的，GGA的PBE单独有相应的赝势；按照赝化电子的数目多少分为大核赝势或是小核的（半芯态赝势），当然（平滑程度的）软、硬也不相同。极化比较严重的体系不适合使用太软的赝势。

PS：上次的计算收敛了没？

作者
Author: zhongyuabc 时间: 2016-6-19 15:38
1.谢谢前辈的解答！
2.还没有计算完，还在计算中。
3.但是我个人觉得不会收敛。

作者
Author: 卡开发发 时间: 2016-6-19 15:48

zhongyuabc 发表于 2016-6-19 15:38
1.谢谢前辈的解答！
2.还没有计算完，还在计算中。
3.但是我个人觉得不会收敛。

看样子估计够呛，实在不行换程序好了，几何优化光是BFGS方法也没有可调整的余地了。

作者
Author: zhongyuabc 时间: 2016-6-19 18:18
好的，谢谢前辈

欢迎光临计算化学公社 (http://bbs.keinsci.com/)