计算化学公社

标题: CO态密度很怪？ [打印本页]

作者
Author: panger 时间: 2016-8-27 00:49
标题: CO态密度很怪？
用CASTEP做点东西，突然发现CO态密度和文献中相差太多了？有经验的同行能帮忙看看这是怎么回事吗？谢谢大家。

作者
Author: panger 时间: 2016-8-27 00:57
自己算的那张图费米能级以上什么都没有？是哪里设的不对？

作者
Author: 卡开发发 时间: 2016-8-27 03:32

panger 发表于 2016-8-27 00:57
自己算的那张图费米能级以上什么都没有？是哪里设的不对？

性质计算中设置的空带数目太少了，导致HOMO以上的能级都出不来，还有CASTEP会把HOMO平移到0eV的位置。-20eV那个深能级要么是文献没放出来，要么是赝势的差异导致的。

作者
Author: panger 时间: 2016-8-27 09:32
谢谢您。那把空带设到100%可以吗？另外HOMO的位置能调整吗？

作者
Author: panger 时间: 2016-8-27 10:08
本帖最后由 panger 于 2016-8-27 10:10 编辑

卡开发发发表于 2016-8-27 03:32
性质计算中设置的空带数目太少了，导致HOMO以上的能级都出不来，还有CASTEP会把HOMO平移到0eV的位置。-20 ...

非常感谢您。还得再麻烦您看看我上层楼的问题，能麻烦指导一下吗？谢谢。

作者
Author: 卡开发发 时间: 2016-8-27 10:47

panger 发表于 2016-8-27 10:08
非常感谢您。还得再麻烦您看看我上层楼的问题，能麻烦指导一下吗？谢谢。

你之前设置多大？这个值不是按照百分数的，而是一个整数，你应该会错意了，是性质计算DOS里面指定的空带数，你可以先设置24（已测试）。

HOMO的位置可以自己读取Fermi能级然后吧图形拷贝到Origin里面去平移。

作者
Author: panger 时间: 2016-8-27 15:13

卡开发发发表于 2016-8-27 10:47
你之前设置多大？这个值不是按照百分数的，而是一个整数，你应该会错意了，是性质计算DOS里面指定的空带 ...

谢谢您。关于HOMO平移的问题。我看输出文件里面费米能级计算的都是负值。那么平移的时候是不是应该把算出来的HOMO[castep默认的位置]向负值方向移动？

作者
Author: 卡开发发 时间: 2016-8-27 16:22

panger 发表于 2016-8-27 15:13
谢谢您。关于HOMO平移的问题。我看输出文件里面费米能级计算的都是负值。那么平移的时候是不是应该把算出 ...

对。
Fermi能级移动到0，意味着原有参考真空能级为0得到的能级E_n-E_Fermi=E'_n。
现在CASTEP给出的是平移后的E'_n，反过来求E_n=E'_n+E_E_Fermi。

作者
Author: panger 时间: 2016-8-27 21:54

卡开发发发表于 2016-8-27 16:22
对。
Fermi能级移动到0，意味着原有参考真空能级为0得到的能级E_n-E_Fermi=E'_n。
现在CASTEP给出的是 ...

真心感谢，谢谢。

作者
Author: panger 时间: 2016-9-1 12:10

卡开发发发表于 2016-8-27 16:22
对。
Fermi能级移动到0，意味着原有参考真空能级为0得到的能级E_n-E_Fermi=E'_n。
现在CASTEP给出的是 ...

还有个问题得麻烦您一下。计算d带中心的时候是应该用平移以后的数据还是平移以前的数据？发现很多结果平移以后基本都在费米能级以下了。那积分到费米能级也就没意义了吧？另外，平移以后计算的值好像偏大很多？希望您给点拨一下。谢谢您。

作者
Author: 卡开发发 时间: 2016-9-1 13:08

panger 发表于 2016-9-1 12:10
还有个问题得麻烦您一下。计算d带中心的时候是应该用平移以后的数据还是平移以前的数据？发现很多结果平 ...

d带中心理论其实我并不是很熟悉，因为我不是做催化的，不过J. Chem. Phys,120(21)里面：

The d-band center was calculated as the first moment of the projected d-band density of states on the surface atoms referenced to the Fermi level.

这个意味着需要进行平移。这个平移是整个能带都平移，不是单独只移动Fermi能级，所以不会改变各能级与Fermi能级的相对位置，千万不要理解错误，因此先积分再平移或者先平移再积分原则上都可以，说白了就是积分时做了个变量代换后边界也动了一下而已。

平移后肯定会变大一些啊，如前面所说

我看输出文件里面费米能级计算的都是负值

Fermi能级移动到0，意味着原有参考真空能级为0得到的能级E_n-E_Fermi=E'_n。

这样移动后E'_n肯定变大。

PS：之前的问题其实建立在某种不严格的讨论，之前说bulk系统的真空能级一般是按照0进行参考的，但是3D周期性的slab模型并非如此，因为z方向还是有周期性的，而z实际计算中并不是无穷大的，所以导致真空能级不是0，而是需要通过求真空静电势来得到。

还有，直接对PDOS积分的话可能有三个因素会影响准确性，1、构成PDOS曲线的点的个数，2、展宽值，3、分割原子球的WS半径，应该尽可能避免。

作者
Author: panger 时间: 2016-9-2 10:30

卡开发发发表于 2016-9-1 13:08
d带中心理论其实我并不是很熟悉，因为我不是做催化的，不过J. Chem. Phys,120(21)里面：
这个意味着需要 ...

谢谢您，那您说的这些影响因素在MS中能够考虑吗？尽量使结果更可信一些。

作者
Author: 卡开发发 时间: 2016-9-2 12:10
本帖最后由卡开发发于 2016-9-2 12:16 编辑

panger 发表于 2016-9-2 10:30
谢谢您，那您说的这些影响因素在MS中能够考虑吗？尽量使结果更可信一些。

哦，这样，我还以为用的是vasp。
1、PDOS绘制能够调整number of per eV，默认是100，你不妨提高看看对结果影响如何；
2、展宽的话，SCF用多大绘制的时候就用多大，这是绘制PDOS时候能选择的最小展宽了吧，SCF处理的时候看看展宽对Entropy的影响，如果Entropy比较小的话，一般对结果影响不是很大；
3、无论是castep还是dmol3处理PDOS的时候分割都是用的Mulliken的方式（这个问题我严格考证过，MS的help文档里面有提到

PDOS calculations are based on Mulliken population analysis, which allows the contribution from each energy band to a given atomic orbital to be calculated. Summation of these contributions over all bands produces a weighted DOS. DMol3 and CASTEP allow you to select the type of weighting required. It is possible, for example, to generate LDOS by adding together all the contributions due to orbitals on a given atom.

），就没有太多选择了吧。

作者
Author: panger 时间: 2016-9-2 12:26

卡开发发发表于 2016-9-2 12:10
哦，这样，我还以为用的是vasp。
1、PDOS绘制能够调整number of per eV，默认是100，你不妨提高看看对结 ...

谢谢您。目前组里目前主要用CASTTEP比较多。不过我总想尝试开发一下CP2K。不知道这款软件做这方面的计算是否合适？您接触过吗？

作者
Author: 卡开发发 时间: 2016-9-2 15:30

panger 发表于 2016-9-2 12:26
谢谢您。目前组里目前主要用CASTTEP比较多。不过我总想尝试开发一下CP2K。不知道这款软件做这方面的计算 ...

cp2k能算，但是我没用过，现在都是再用QE。

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