Does anybody here know how to do spin orbital coupling through VASP? Please can anybody send me some important file as INCAR KPOINTS and other files that are needed.
For SOC calculations LSORBIT tag works fine, but this tag will not work with the normal vasp setting in the makefile. Here you need to compile the vasp for full fft mashes not for reduced mashes in some direction or Gamma-point.
So compile your vasp without -DNGXhalf and -DNGZhalf flags and then run your calculations with LSORBIT tag.
How do I compile vasp without _DNGXhalf and -DNGZhalf flags before? While running it says ''non collinear calculations require that VASP is compiled without the flag -DNGXhalf and -DNGZhalf .
For older version of vasp (5.3.5 or less) open your Makefile and look for -DNGXhalf, -DNGZhalf and remove them, only these particular tags you should remove and compile it again. Once compiled you can use it.
For newer version (which you can download from VASP, licensed one) vasp 5.4.1 or 5.4.4 you just compile as "make all" or "make ncl" (non-collinear) you will get an executable in "bin" folder (created inside vasp folder) named as "vasp_ncl" that would be enough for doing SOC calculations.
The last point to say "make all" will give you there executables, vasp_std for normal vasp calculations, vasp_ncl for soc calculations and vasp_gam for gamma point calculations. It's better to do make all, it may take some time to compile all
I have tried to perform the spin orbital coupling calculation in cobalt hcp( 2atoms/cell), using the file INCAR2, but nothing happens, I don't have any datas in the OUTCAR and OSZICAR, what's the problem please?
The code in these INCARs take you to the mag. ground state (ISPIN=1). what if we skip going on ground state and avoid taking extra relaxation steps(NSW=24)?
second, instead of taking LREAL=F, if we consider LREAL=Auto. then ?
What is the right way to do spin orbit coupling (SOC) calculation ? We should relax system with SOC or we should relax without SOC and then single point calculation (nscf) for spin orbit coupling effect.
For doing spin orbit coupling (SOC) calculation we should relax system with spin polarization (ISPIN=2) and LORBIT=11 or without spin we should relax the system.
What about LORBMOM tag for orbital and spin moments. we should give this tag.
Yes, you have to relax with ISPIN = 2 for magnetic systems. Relaxation does not depends on LORBIT. LORBIT tag is to get density of states. If you want to plot DOS, then yes you need LORBIT = 10 or 11. For details, see VASP manual.
LORBMOM = .TRUE. is needed to get orbital moments when SOC is added.
Should the symmetry operations get changed on the inclusion of SOC?
When, i am performing band structure calculations with spin-orbit coupling (SOC) then, its symmetry operations are getting reduced than the one without SOC.
Suppose in calculations without SOC, we have 12 symmetry operations and with SOC, it gets reduced to 2.
So, can somebody please put some light whether symmetry operations should change with the inclusion of SOC or not?
If you want to confirm your system that it is ferro or anti ferro magnetic you need to check the long range and short range magnetic ordering using, AFM and FM calculations. You need to run the calculations with AFM and FM (MAGMOM) tags and compare the stable states from energy comparison.
For the confirmation of Ferro and Anti-ferro magnetic system.
Step 1. You will add LORBIT = 11 in the INCAR-file and relax the system.
Step 2. After the completion of step1 .You will check the OUTCAR-file and find the magnetization of each orbital (see the example below). If tot is positive (+) and have some value also than, it will represent Ferro magnetic or If tot is negative (-) and have some value also than, it will represent Anti-ferro magnetic.
Payal Wadhwa Its normal to see reduced symmetry operations and safe side it to switch off symmetry since number of k-points (which depend on sym) do change for different magnetic orientation which might verse final outcome.
Hello! I tried SOC (scf calculations) calculations after relaxing the structure, I'm only able to get DOS for up spin (even though ISPIN = 2) but if I do the same thing removing SOC, I get DOS for both the spins. Can anyone help me?
Sromona, when you include spin orbit coupling effects, spin is no longer a single quantum number, but is instead coupled with the angular momentum (due to relativistic effects) and so there is no spin up or spin down, just total spin. This might help: http://xbeams.chem.yale.edu/~batista/vvv/node26.html. Note that SOC effects are more pronounced for structures/molecules containing heavy atoms (~Z>50).
How can I get DOS for both the spins for a SOC system? I have seen people doing soc calculation for a spin polarized system. What should I do to do that?