GGA+U+SOC

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Post Reply
Message
Author
vrinsom
Newbie
Newbie
Posts: 5
Joined: Sat Dec 14, 2019 4:23 am

GGA+U+SOC

#1 Post by vrinsom » Thu Aug 11, 2022 2:28 pm

Hello,
I am observing some strange behavior when I try to run spin orbit coupling calculations. I am running test calculations on Ni-doped CeO2 and Ni-doped CeO2 with an oxygen vacancy at the nearest-neighbor site to the Ni dopant. The POSCARs for both the runs are identical, except that the latter run has an oxygen vacancy next to the Ni.

My settings for both the calculations are exactly the same. I generate the WAVECAR and CHGCAR from a collinear run using 1 k-point at Gamma (I use the vasp_std binary file). I do not relax the structure. I then run the non-collinear calculation reading the WAVECAR and CHGCAR from the collinear run and at 1 k-point at Gamma.

However, for the non-collinear calculation in the Ni-doped case (without oxygen vacancy), I receive the error
ERROR: while reading WAVECAR, plane wave coefficients changed 31318
15659
But for the non-collinear calculation in the Ni-doped case with the oxygen vacancy, VASP is able to read the WAVECAR correctly and successfully completes the calculation.

This is strange, since the only difference in the calculations is the presence of an oxygen vacancy ! I was wondering why this occurs, and how do I get VASP to read the WAVECAR for the runs without an oxygen vacancy?

alexey.tal
Global Moderator
Global Moderator
Posts: 228
Joined: Mon Sep 13, 2021 12:45 pm

Re: GGA+U+SOC

#2 Post by alexey.tal » Fri Aug 12, 2022 8:32 am

Hi,

Could you please provide all relevant input and output files for the steps you have described according to the forum guidelines.

vrinsom
Newbie
Newbie
Posts: 5
Joined: Sat Dec 14, 2019 4:23 am

Re: GGA+U+SOC

#3 Post by vrinsom » Thu Aug 18, 2022 1:23 pm

Sure, here are the required input and output files. Sorry for the delay, thank you for looking into this!
You do not have the required permissions to view the files attached to this post.

alexey.tal
Global Moderator
Global Moderator
Posts: 228
Joined: Mon Sep 13, 2021 12:45 pm

Re: GGA+U+SOC

#4 Post by alexey.tal » Thu Aug 18, 2022 4:38 pm

This problem occurs because of the number of NBANDS. If NBANDS is not set, it will be calculated based on the initial magnetic moments, number of electrons, number of atoms and CPU cores.

In order to be able to read WAVECAR produced in a collinear calculation, the non-collinear calculation should include exactly double the number of bands.

vrinsom
Newbie
Newbie
Posts: 5
Joined: Sat Dec 14, 2019 4:23 am

Re: GGA+U+SOC

#5 Post by vrinsom » Fri Aug 19, 2022 11:00 pm

Got it ! Thank you so much !

Post Reply