  Using numerical simulations that mimic recent experiments on hexagonal
colloidal ice, we show that colloidal hexagonal artificial spin ice exhibits an
inner phase within its ice state that has not been observed previously. Under
increasing colloid-colloid repulsion, the initially paramagnetic system crosses
into a disordered ice-regime, then forms a topologically charge ordered state
with disordered colloids, and finally reaches a three-fold degenerate, ordered
ferromagnetic state. This is reminiscent of, yet distinct from, the inner
phases of the magnetic kagome spin ice analog. The difference in the inner
phases of the two systems is explained by their difference in energetics and
frustration.
