Using the recently developed effective interaction potentials between polyelectrolyte stars, we examine the structure and phase behavior of solutions of the same. The effective interaction is ultrasoft and density dependent, owing to the integration of the counterionic degrees of freedom. The latter contribute extensive volume terms that must be taken into account in drawing the phase diagram of the system. The structural behavior of the uniform fluid is characterized by anomalous structure factors, akin to those found previously for solutions of uncharged star polymers. The phase diagram of the system is very rich, featuring a fluid phase at low arm numbers of the stars, two reentrant melting regions, as well as a variety of crystal structures with unusual symmetry. The physical origin of these features can be traced back to the ultrasoft nature of the effective interaction potential. (C) 2004 American Institute of Physics.