Background: The inner crust of a neutron star is believed to consist of Coulomb-frustrated complex structures known as "nuclear pasta" that display interesting and unique low-energy dynamics. Purpose: To elucidate the structure and composition of the neutron-star crust as a function of temperature, density, and proton fraction. Methods: A new lattice-gas model, the "Charged-Ising Model" (CIM), is introduced to simulate the behavior of neutron-star matter. Preliminary Monte Carlo simulations on 30^3 lattices are performed for a variety of temperatures, densities, and proton fractions. Results: Results are obtained for the heat capacity, pair-correlation function, and static structure factor for a variety of conditions appropriate to the inner stellar crust. Conclusions: Although relatively simple, the CIM captures the essence of Coulomb frustration that is required to simulate the subtle dynamics of the inner stellar crust. Moreover, the computationally demanding long-range Coulomb interactions have been pre-computed at the appropriate lattice sites prior to the start of the simulation resulting in enormous computational gains. This work demonstrates the feasibility of future CIM simulations involving a large number of particles as a function of density, temperature, and proton fraction.

In this letter we show that the presence of the long-range Coulomb force in dense stellar matter implies that the total charge cannot be associated with a chemical potential, even if it is a conserved quantity. As a further consequence, the analytical properties of the partition sum are modified, changing the order of the phase transitions and affecting the possible occurrence of critical behaviours. The peculiar thermodynamic properties of the system can be understood introducing a model hamiltonian in which each charge is independently neutralized by a uniform background of opposite charge.