2D polarization calculation

[pooja_rani]

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

[zahedzx]

Hi Pooja,

The warning

Code: Select all

PEAD_POLARIZATION_CALC: WARNING: contribution from string ... not well clustered

indicates that the Berry phases from different k-point strings are not sufficiently consistent. This is usually related to k-point sampling, insufficient wavefunction convergence, or a very small band gap (near-degenerate occupied states), rather than the DIPOL tag itself.

The "spin channel 1" message is normal even for a non-spin-polarized calculation (ISPIN=1) and is not an indication that the calculation is magnetic.

Regarding the non-zero polarization: in the modern Berry-phase formalism, the absolute polarization is only defined modulo a polarization quantum. Therefore, a centrosymmetric structure does not necessarily give a printed value of exactly zero. The physically meaningful quantity is usually the change in polarization between two structures.

For polarization calculations, it is recommended to first perform a fully converged SCF calculation and then run a static calculation with LCALCPOL = .TRUE. using the converged charge density/wavefunctions.

Best,
Zahed