In xc-pot/exchpbe.f90?ref_type=heads, I found that it is not designed for a spin polarized system. However, when I choose to do a noncollinear calculation, how could this work correctly? I mean that I only find this code about pbe exchange energy, but it seems to be incorrect. I didn't add libxc command in my input document, so I assume that the functional of my calculation uses functionals from xcpot.
Best wishes,
Yule

Dear Gregor,
Thank you for your answer, and problems are solved.
Sincerely,
Yule