Dear CP2k community,<br><br>  I am wondering how the paramagnetic spin-orbit (PSO) operator is exactly defined in CP2k codes. In text books, the PSO is given by <br><br>                                            (Summation_i)(Summation_j) (alpha*alpha)*(Zi/(Rij^3))*( Rj X Pj )<br>
 <br>  where alpha is the fine structure constant, Zi is the nuclear charge , Rj is the electron position and Pj is the electron momentum.<br><br>  If cp2k uses the same formula, the spin-orbit coupling energy should be something like Tr(Density_matrix * PSO_matrix).<br>
<br>  Is "CP2k_INPUT/FORCE_EVAL/DFT/PRINT/AO_MATRICES/PSO" supposed to print out the PSO_matrix?<br><br> Thank you very much.<br><br><br>Hanning Chen<br><br>Northwestern University<br>Evanston, IL 60208 <br> <big class="uctt"><a href="http://cp2k.berlios.de/manual/CP2K_INPUT.html"></a></big>