A stable 4f half-metallic ferromagnetism Eu 4 O 3 N predicted from GGA+U calculations