求助-vasp该怎么计算滑移铁电的二维能量差异图

1070548940

figure1 d图是figure1b图中沿白色实线方向滑移的能量，其中横坐标0到1/3处的情况，是铁电材料中的双势阱。我清楚双势阱这个曲线如何用vasp计算得到，但是figure1b图这种二维的滑移能量总体情况，我该怎么用vasp计算得到。难道是类似算AB实线的这情况，用沿着不同的滑移路径将整个二维平面覆盖，然后算每个滑移路径的能量最后组合成这个二维能量差的情况吗？
文献链接:https://pubs.acs.org/doi/10.1021/acs.nanolett.2c03373?ref=pdf
如有表述不清还请见谅，恳请大家指点迷津！