Charge density functional plus U calculation of lacunar spinel GaM4Se8(M = Nb, Mo, Ta, and W). (25th March 2019)