计算化学公社

标题: 替同学发帖：如何用VASP算XPS谱? [打印本页]

作者
Author: 幻七熏 时间: 2015-11-19 16:27
标题: 替同学发帖：如何用VASP算XPS谱?
我想用ICORELEVEL=2来算。
第一： CLNT=species，这个值怎么确定的，是POSCAR中元素的排位顺序吗？排在第几，就是几？
第二， CLL值是只能等于或者大于0的吗？算Co 2p的时候，该值等于1还是-1？
第三，算XPS的原子需要relax,其他原子需要固定吗？还是跟正常算吸附时一样，表面几层relax，下边几层Fix？
第四：CLZ=0.5是激发几个电子啊？为什么书册上建议是0.5或者1?
请大神们指点迷津，非常感谢！

作者
Author: 卡开发发 时间: 2015-11-19 18:56
没算过XPS，不过：
1、一般vasp的species应该是按照POSCAR元素顺序；
2、角量子数应该是非负整数；
3、按照“the other core electrons are not allowed to relax, which might cause a slight error in the calculated energies”，这种情况结构弛豫可能没啥意义，我的理解是ibrion=-1就行；
4、“Usually 1 or 0.5 (Slaters transition state) are sensible choices”这个激发能估计是利用Koopmans定理或是Janak定理得到的。

作者
Author: 起就起 时间: 2015-11-24 14:09

卡开发发发表于 2015-11-19 18:56
没算过XPS，不过：
1、一般vasp的species应该是按照POSCAR元素顺序；
2、角量子数应该是非负整数；

谢谢，之前去验证你说的了，所以回复晚了。非常感谢，你说的都是对的，太厉害了。不过我用VASP算了，跟实验都对不上，连趋势都对不上，我也是无语了。还有两个问题，第一，那个IBRION，你通常都用什么值啊，是2，还是1呢？我看1是best choice, 2是reliable backup routines.我一直都用的2；第二，我用IBRION=-1和用IBRION不为零，但是设置NSW=0是一样的吧？都不会动的哈。

作者
Author: 卡开发发 时间: 2015-11-24 14:23

起就起发表于 2015-11-24 14:09
谢谢，之前去验证你说的了，所以回复晚了。非常感谢，你说的都是对的，太厉害了。不过我用VASP算了，跟实 ...

和实验对不上我不太好说，不知道您有没有进行收敛性测试，另外不知道你的体系泛函、方法的选择影响有多大。
1、IBRION我自己一般使用2，通常这种方法还行；
2、我的理解IBRION=-1和IBRION>0，NSW=0没啥区别。

作者
Author: 起就起 时间: 2015-11-24 15:51

卡开发发发表于 2015-11-24 14:23
和实验对不上我不太好说，不知道您有没有进行收敛性测试，另外不知道你的体系泛函、方法的选择影响有多大 ...

恩，谢谢。请问什么叫收敛性测试，是测试参数吗？

作者
Author: 万里云 时间: 2015-11-24 16:01
本帖最后由万里云于 2015-11-24 16:05 编辑

XPS涉及芯电子激发，用全电子程序算得准一些。

推荐exciting/elk，exciting专门有个例子就是算XPS的。
PS：刚才去看了，exciting那个例子是用BSE算得X-ray absorption spectrum，估计计算量惊人。

XPS实际上类似于ARPES，要用GW方法算谱函数。要算到上百eV范围的谱函数，估计计算量更惊人。

作者
Author: 起就起 时间: 2015-11-24 16:21

万里云发表于 2015-11-24 16:01
XPS涉及芯电子激发，用全电子程序算得准一些。

推荐exciting/elk，exciting专门有个例子就是算XPS的。

我要算的不是全电子谱，而是core level shift.所以不用算谱，就算个相对值。这个是可以用VASP算的吧？但是我都对不上实验。

作者
Author: 起就起 时间: 2015-11-24 16:24

起就起发表于 2015-11-24 16:21
我要算的不是全电子谱，而是core level shift.所以不用算谱，就算个相对值。这个是可以用VASP算的吧？但 ...

PS. 这个exciting/elk是什么啊？VASP中的方法吗？我不知道

作者
Author: 万里云 时间: 2015-11-24 16:29

起就起发表于 2015-11-24 16:24
PS. 这个exciting/elk是什么啊？VASP中的方法吗？我不知道

和VASP一样的软件。

作者
Author: 万里云 时间: 2015-11-24 16:29

起就起发表于 2015-11-24 16:21
我要算的不是全电子谱，而是core level shift.所以不用算谱，就算个相对值。这个是可以用VASP算的吧？但 ...

全电子只是让这个值算得准。

作者
Author: 起就起 时间: 2015-11-24 18:15

万里云发表于 2015-11-24 16:29
全电子只是让这个值算得准。

谢谢各位哈

作者
Author: 卡开发发 时间: 2015-11-24 19:17

万里云发表于 2015-11-24 16:29
全电子只是让这个值算得准。

@起就起原则上vasp处理core level可能都得把内层电子引进来，BSE确实计算量惊人啊

作者
Author: 万里云 时间: 2015-11-24 19:39

卡开发发发表于 2015-11-24 19:17
@起就起原则上vasp处理core level可能都得把内层电子引进来，BSE确实计算量惊人啊

VASP的GW和BSE资源占用都十分惊人不说，关于计算参数还有很多古怪的规定。

写程序的人不花时间处理程序逻辑，用程序的人就得花时间了。

作者
Author: 卡开发发 时间: 2015-11-24 19:42

万里云发表于 2015-11-24 19:39
VASP的GW和BSE资源占用都十分惊人不说，关于计算参数还有很多古怪的规定。

写程序的人不花时间处理程 ...

最蛋疼的是手册的说明似乎也不全面，使用者多数得靠猜了。

作者
Author: 万里云 时间: 2015-11-24 19:43
http://exciting-code.org/boron-x-ray-absorption-spectra-using-bse
楼主看看这个吧。

作者
Author: 万里云 时间: 2015-11-24 20:23

卡开发发发表于 2015-11-24 19:42
最蛋疼的是手册的说明似乎也不全面，使用者多数得靠猜了。

VASP的BSE计算还是experimental feature，我当初在论坛考了很久的古才找到的。

作者
Author: 卡开发发 时间: 2015-11-24 20:25

万里云发表于 2015-11-24 20:23
VASP的BSE计算还是experimental feature，我当初在论坛考了很久的古才找到的。

其实下面有个DMFT.f90，感觉里面有秘密，不过我始终找不到相关信息。

作者
Author: 万里云 时间: 2015-11-24 20:38

卡开发发发表于 2015-11-24 20:25
其实下面有个DMFT.f90，感觉里面有秘密，不过我始终找不到相关信息。

没有对应的说明文档，源码就像天书一样。

作者
Author: 快乐小王子 时间: 2024-7-27 16:23

万里云发表于 2015-11-24 19:43
http://exciting-code.org/boron-x-ray-absorption-spectra-using-bse
楼主看看这个吧。

这个链接打不开啊

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