I have tried to calculate for a system with spin-orbit coupling (SOC) using SIESTA and Wannier90, but an error
non-collinear spin not implemented yet
appears in a SIESTA output file. I would be very grateful if you have some time or knowledge to solve this error or suggest any clue for that.
Here is what I have tried so far.
- Siesta: 5.0.0
- Wannier90: 3.1.0
- Input files: graphene.fdf, graphene.win from this tutorial
- Pseudopotential file: C.psml (ONCVPSP, PBE, fully relativistic) from Pseudo-Dojo
Graphene is only used as an example to check the spin-orbit calculation. Fisrt, I calculated as the input files are written, then the calculation completed without any errors.
Next, I added Spin spin-orbit
to the fdf file to include SOC. The calculation couldn't be completed and an error
getFileNameRoot: non-collinear spin not implemented yet
was displayed. Comparing the output file without SOC, the error is in the line where SIESTA reads nnkp file for Wannier90. It should have been
siesta2wannier90: Reading the graphene.nnkp file,
instead of the above "non-collinear" error which only occurs when SOC is enabled.
If I calculate with SOC without Siesta2Wannier90
lines (see the last part of the fdf file), there is no error. So my conclusion so far is that the cause is in the Wannier90 or the combination of SIESTA and Wannier90, especially in the nnkp file.
Adding spinors = T
to the win file and recreating nnkp file before starting the SIESTA calculation with SOC did not make any difference.
Is there anyone who can suggest where to correct or which command is needed to complete this calculation? I appreciate your help in advance.