It depends on what kind of system you are dealing with and what you want to calculate, but for example the DFT codes VASP (plane wave basis) and the GREEN-SIESTA implementation of the SIESTA code (atomic orbital basis) can deal quite successfully with SOC. The KKR code by the Budapest-Vienna group, for example, is fully relativistic and solves the Dirac equation directly, so the SOC is properly addressed there. Many other modern DFT codes also include SOC implementations (WIEN2k, Quantum Espresso, ABINIT, etc.), but I'm not so familiar with them.

https://www.vasp.at/

http://www.icmm.csic.es/jcerda/EHT_TB/Green/Guide/green-guide.html#green-siesta

http://journals.aps.org/prb/abstract/10.1103/PhysRevB.51.9552

Hi Maryam,

Spin-orbit coupling can be applied fully self-consistently (two-component ZORA Hamiltonian) in both the molecular DFT code ADF and the periodic DFT code BAND.

Furthermore, SOC in ADF can be used for many spectroscopic properties (NMR, EPR, UV/ VIS, XAS) and one can calculate spin-orbit coupling matrix elements between excited states.

You can try out ADF & BAND (the whole ADF modeling suite) for 30 days from the SCM website through the link below.

Hope this helps, let me know if you have further questions

Fedor

https://www.scm.com/trial

ORCA can handle it both on DFT and ab initio levels, if I'm not mistaking

Thanks every body

I will try the links you provided me with

@T. P M Goumans, How to find the result of SOC from ADF out file? Would you like to help me please? Regards

There are many tutorials and examples online. You can also contact support if you can't figure it out.

Does NWChem support TDDFT with ZORA correction?

GAMESS-US is the best choice!

Have a try with Molcas