<div>Dear CP2K community,</div><div>I have been trying to benchmark/test the spin orbit coupling capabilities in CP2K, by measuring the spin orbit splitting of the valence bands of Silicon at the gamma point. <br />These calculations were performed on CP2K v2025.1 using the GTH_SOC_POTENTIALS and TZV2P MOLOPT basis.</div><div><br />The values I get are much smaller than the reported theoretical literature and experimental values. From this run, CP2K gives 22 meV of splitting. I have attached the input file (si.inp) and the output bandstructure file (bandstructure_si). </div><div>The same structure ran through QE with fully relativistic pseudopotentials gives 48 meV of splitting; which is closer to the experimental value of ~43 meV. <br /><br />It would really be helpful if someone can clear it up if I have made a mistake in the input/calculation or this is a known limitation of the pseudopotential/method used in CP2K?</div><div><br /></div><div>Best regards,</div><div>Nisarg.</div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/45561b15-4253-4176-b982-1f8d105a3e6en%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/45561b15-4253-4176-b982-1f8d105a3e6en%40googlegroups.com</a>.<br />