Theoretical study of deuteron polarization in γ+d → π0+d