Charge ordering in the three-dimensional Coulomb glass at finite temperatures and low disorders.

In this paper, we have studied the three dimensional Coulomb glass lattice model at half-filling using Monte Carlo Simulations. Annealing of the system shows a second-order transition from paramagnetic to charge-ordered phase for zero as well as small disorders. We have also calculated the critical exponents and transition temperature using a finite sizing scaling approach. The Monte Carlo simulation is done using the Metropolis algorithm, which allowed us to study larger lattice sizes. The transition temperature and the critical exponents at zero disorder matched the previous studies within numerical error. We found that the transition temperature of the system decreased as the disorder is increased. The values of critical exponents $\alpha$ and $\gamma$ were less and value of $\nu$ more than the corresponding zero disorder values. The use of large system sizes led to the correct variation of critical exponents with the disorder.