<div dir="ltr"><div>Dear WIEN2k users</div><div><br></div><div>Hi. I am studying some materials</div><div>which have strong spin-orbit coupling (SOC) effects.</div><div><br></div><div>What I want to do is</div><div>a modeling with simple SOC form, such as</div><div>H_so → ξvec[L]·vec[S].</div><div><br></div><div>I read some lecture notes, and noticed that</div><div>WIEN2k works with Dirac relativistic treatment in LAPWSO.</div><div><br></div><div>Question1)</div><div>Can LAPWSO print out vec[L] (i.e., with three components) for each k-point?</div><div><br></div><div>==================================================</div><div><br></div><div>I also realized that one can get vec[L] for each k-point from LAPWDM (with some modification of the code).</div><div><br></div><div>Question2)</div><div>With 「vec[L] & vec[S] from LAPWDM」 and 「ξ from LAPWSO」,</div><div>is it possible to make a model, H_so → ξvec[L]·vec[S]?</div><div><br></div><div>==================================================</div><div><br></div><div>Any comment will be very helpful for me.</div><div>Thank you for reading this e-mail,</div><div>and have a nice day.</div><div><br></div><div>- Kyohn</div></div>