The phase diagram of star polymer solutions in a good solvent is obtained o
ver a wide range of densities and arm numbers by Monte Carlo simulations. T
he effective interaction between the stars is modeled by an ultrasoft pair
potential which is logarithmic in the core-core distance. Among the stable
phases are a fluid as well as body-centered cubic, face-centered cubic, bod
y-centered orthogonal, and diamond crystals. In a limited range of arm numb
ers, reentrant melting and reentrant freezing transitions occur for increas
ing density.