First order transition in two dimensional coulomb glass

We have studied the ground states of two dimensional lattice model of Coulomb Glass via Monte Carlo annealing. Our results show a possibility of existence of a critical disorder (W_c) below which the system is in the charge ordered phase and above it the system is in the disordered phase. We have used finite size scaling to calculate W_c = 0.2413, the critical exponent of magnetization β = 0 indicating discontinuity in magnetization and the critical exponent of correlation length ν = 1.0. The distribution of staggered magnetization for different disorder strengths shows a three peak structure. We thus predict that two dimensional Coulomb Glass shows a first order transition at T=0.