Dear all, I am reading the codes. The way to calcualte XYZMAT implemented in opeigr.F seems not being described in PRB-61-10040. Is there anybody can give some hints to understand how opergr.F works? Thanks, Zhenyu