Weldingspock 发表于 2024-4-26 16:45
老哥请问一下，怎么用berry phase计算极化强度呀。我有两个体系，一个是slab一个是bulk结构。真心求教

你好，请问您解决了吗，我也真心请教，如果您有时间能否帮我？15293064229

老哥请问一下，怎么用berry phase计算极化强度呀。我有两个体系，一个是slab一个是bulk结构。真心求教

NUAAcxk 发表于 2023-9-11 15:38
请教楼主后来知道跳变的原因了吗

程序bug 换vasp 6.2之后的版本就行

请教楼主后来知道跳变的原因了吗

exatude 发表于 2023-7-4 13:47
echo -e "400 /n 1 1 0 /n -1 1 0 /n 0 0 1 /n" | vaspkit
如果装了vaspkit的话，用这行命令可以把三方 ...

好谢谢。 算出来确实是根号三关系

不过这样变成正交晶胞和之前好像没什么区别。。 还是要x轴结果除以cos30°得到原来对角线方向极化强度。如果我把b轴变到y轴上，这样ab轴是在xy轴上，但这样再移动原子会比之前麻烦不少。。 用vtst的nebmake.pl给的过渡态不对

helpcal 发表于 2023-7-4 09:03
我的a轴就是沿着x轴的，这样应该可以？对角线方向是用x方向的结果除以cos30°。

我之前尝试过重新弄 ...

echo -e "400 /n 1 1 0 /n -1 1 0 /n 0 0 1 /n" | vaspkit
如果装了vaspkit的话，用这行命令可以把三方晶胞重新设置为正交的。
你应该也算出来x方向和y方向极化强度绝对值不一样，而且比值大概是根号三附近吧，我说的a,b方向和x,y方向定义不一样就是这个意思。

Deepast 发表于 2023-7-4 09:38
DIPOL 有什么讲究吗，我试过自己设会不收敛

不清楚。手册解释的没有特别清楚。不过统一就行。不收敛是位置选择问题，测试几个选个收敛的就行。

DIPOL 有什么讲究吗，我试过自己设会不收敛

exatude 发表于 2023-7-3 14:23
最近也在算这种东西。如果你用的是图中这种三方晶胞来算的话，VASP给出的极化强度里的x，y方向和晶胞的a，b ...

我的a轴就是沿着x轴的，这样应该可以？对角线方向是用x方向的结果除以cos30°。

我之前尝试过重新弄成a b轴正交的晶胞但没找到好方法弄起来很麻烦就后来没弄了。步骤大概就是基矢变换一下然后手动添加原子。。想问下有没有更有效的方法。

最近也在算这种东西。如果你用的是图中这种三方晶胞来算的话，VASP给出的极化强度里的x，y方向和晶胞的a，b方向并不对应，建议先redefine lattice成正交的晶胞再算。