Inner Phases of Colloidal Hexagonal Spin Ice

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.