Page 1 of 1

VASP + WANNIER90 + spin orbit coupling

PostPosted: Tue Feb 05, 2019 9:25 am
by ew
Hi all,

I have a question concerning the combination of VASP and WANNIER90.

My aim is to calculate Hamiltonian Matrix elements in the Wannier function basis for different materials with spin orbit coupling.

As an example let's take mercury telluride (HgTe).

The vasp calculations work fine, but when it comes to the Wannier90 calculation (LWANNIER90_RUN = .TRUE. in the INCAR file), things start getting tricky.

In comparison to non-SOC calculations I just double the number of wannier functions (nuw_wann). But then my calculation does not stop properly and I do not get any output file (despite setting write_hr = T).

In my file I choose num_wann = 16 and spinors = .true. with the following projections on orbital basis functions:

Begin Projections
Hg : sp3
Te : sp3
End Projections.

But then in my log file I get a mysterious message:

Projection 26 done.
Projection 27 done.
Projection 28 done.
Projection 29 done.
Projection 30 done.
Projection 31 done.
Projection 32 done.
MLWF_WANNIER90: ERROR: number of projections not equal to num_wann 32

I do not understand why the projections go higher than 16 or the meaning of the error message.

Does anyone of you have experience in including spin orbit coupling in a VASP+WANNIER90 calculation and can explain me what I am doing wrong?

Thank you very much in advance!

Best regards

Re: VASP + WANNIER90 + spin orbit coupling

PostPosted: Mon Feb 25, 2019 2:59 pm
by d.pshenay
I had the same problem (I used vasp.5.4.4 with gnu compillers and wannier 1.2).
Try to switch off parallelization parameters in INCAR.
In my case it helped me (Experimentally, I found that I can only left k-point parallelization).

Re: VASP + WANNIER90 + spin orbit coupling

PostPosted: Fri Mar 01, 2019 4:15 pm
by ew
Hey, thank you very much for your reply.

Recently I also did some tests with wannier 1.2 and achieved some promising results,
I'll keep testing with your parallelization hint.

However, are you sure that your calculations are correct when you turn on spin orbit coupling? There are some posts in the forum that
suggest that VASP and WANNIER90 cannot work together at all (e.g. viewtopic.php?f=3&t=16899).
If that's the case I may have to find an alternative way, maybe with the LOCPROJ-parameter.

Best regards

Re: VASP + WANNIER90 + spin orbit coupling

PostPosted: Mon Mar 04, 2019 5:04 pm
by d.pshenay
Yes, I saw that. But, how can I check if it works OK? - I used it with SOC for band interpolation and compared obtained band structure with DFT results. They fit each other pretty close (to withing meV inside frozen window). My similar calculations using QuantumEspresso+Wannier90 interpolation with SOC also gave the same band structure.