The script appears to work. Are you sure the location of the output NC file is correct?
BTW, 2 comments:
- You should use a pseudopotential+basis set that is suitable for spin-orbit, i.e. SG15SO.
- GeometryOptimization is usually done without spin-orbit.