替同学发帖：如何用VASP算XPS谱?

幻七熏

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

卡开发发

没算过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定理得到的。

起就起

卡开发发 发表于 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是一样的吧？都不会动的哈。

卡开发发

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

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

起就起

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

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

万里云

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

推荐exciting/elk，exciting专门有个例子就是算XPS的。
PS：刚才去看了，exciting那个例子是用BSE算得X-ray absorption spectrum，估计计算量惊人。
XPS实际上类似于ARPES，要用GW方法算谱函数。要算到上百eV范围的谱函数，估计计算量更惊人。

起就起

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

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

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

起就起

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

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

万里云

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

和VASP一样的软件。

万里云

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

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

起就起

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

谢谢各位哈

卡开发发

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

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

万里云

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

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

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

卡开发发

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

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

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

万里云

http://exciting-code.org/boron-x-ray-absorption-spectra-using-bse
楼主看看这个吧。