极性slab模型偶极矩矫正的问题

wernea

请教一下各位大侠，极性slab模型偶极矩矫正的问题，我所了解的大概步骤如下：
①切面
②第一次优化表面
③第一次静态自洽，同时设置LVTOT = .TRUE.，输出静态势，画图，看真空层是否倾斜，平整代表体系无极化，否则有极化，需进行下列矫正偶极矩的步骤（输出静态势是静态自洽这一步吗？据说是用一个程序取出c方向的静态势的，这个程序有哪位亲知道，还请告知。）
④第二次优化表面（cp contcar poscar），这次优化需要偶极矫正，即设LDIPOL= .TRUE. ；IDIPOL = 3 （真空层加在c方向）；DIPOL = center of cell（看manual中是说要同时设置这个DIPOL = center of cell，不知道有没有理解对，有的人说不需设置DIPOL这个参数，vasp会自己找这个DIPOL的位置？）
⑤第二次静态自洽，也需要设置偶极矫正参数。这样到这一步可得到可靠的表面能？
以上是我所了解的，有不正确的地方还请帮我指出，拜谢！