求助：vasp berry phase算slab面内极化强度突变

helpcal

请问有没有人遇到过这种极化强度？这是slab的面内极化强度(x方向）。中间态很平稳，但是到2个参考态铁电和顺电相极化强度突变（斜率突增）。按中间态的变化趋势粗略估算得到的面内极化强度和文献差别大概4.7%，但直接铁电减顺电结果都2倍了。。

附上结构图以及从铁电到顺电原子移动方向和berry phase算的极化强度结果。polarization quantum用的是e*x方向晶格常数。（这里统一用dipole moment 单位）

补充一下新细算的结果，在99.2-99.3%处发生突变，变化值接近polarization quantum的一半（和其一半还是有约0.08e*A的不同，约为其4%的值）不清楚这是不是只是巧合还是有什么其他原因。

相关参数：
# Dipole correction
LDIPOL = .TRUE.
IDIPOL = 3
DIPOL = 0.5 0.5 0.6

LCALCPOL = .TRUE.


（之前搜索过DIPOL，网上对其没有特别统一的方案，而且也有说真正的dipol位置是在自己设置的值+0.5z，但vaspwiki上我没看到这点。就随便选了个0.6。但不管怎么说这里用berry phase算极化强度应该需要设置统一的DIPOL）

zhenjun

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

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

Weldingspock

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

helpcal

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

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

NUAAcxk

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

helpcal

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给的过渡态不对

exatude

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方向定义不一样就是这个意思。

helpcal

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

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

Deepast

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

helpcal

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

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

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

exatude

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