Hi,
I am trying to perform berry phase polarization calculation for a 2D material. I used the tag LCALCPOL, but every time i am getting some warning in the output, which i attached in the form of a screenshot here. I tried to use a dense K-grid and also used the other method by using LBERRY tags, but every time I get the same warning. I also use the DIPOL tag with LCALCPOL (like DIPOL = 0.5 0.5 0.144) and tried to vary these values, but I am still getting the same warning. Can you please help me with this problem? The interesting thing in that is it's printing the dipole values in OUTCAR, but that is incorrect. So, I am performing a calculation for a trilayer with centrosymmetric stacking, and the polarization is actually zero, but in OUTCAR I am getting some values, which I think is incorrect. This I am doing to compare with two bilayer polarizations. I hope you understand what i am trying to do.
Also, this is a nonmagnetic material. I do not understand why this is writing "spin channel" in the output. I started with the calculation for the example that is given with the LCALCPOL tag. The calculation was completed, but it did not print any dipole moment in the INCAR file. Can you please correct me if I missed something in this calculation?
Besides that, I have a query, do we need to perform SCF calculation before the Berry phase calculation? I just want to know the proper procedure to get the polarization of a 2D material.
Thank you and kind regards,
Pooja