A systematic search for multicomponent crystal structures is carried out for five different ternary systems of nuclei in a polarizable background of electrons, representative of accreted neutron star crusts and some white dwarfs. Candidate structures are "bred" by a genetic algorithm, and optimized at constant pressure under the assumption of linear response (Thomas-Fermi) charge screening. Subsequent phase equilibria calculations reveal eight distinct crystal structures in the $T=0$ bulk phase diagrams, five of which are complicated multinary structures not before predicted in the context of compact object astrophysics. Frequent instances of geometrically similar but compositionally distinct phases give insight into structural preferences of systems with pairwise Yukawa interactions, including and extending to the regime of low density colloidal suspensions made in a laboratory. As an application of these main results, we self-consistently couple the phase stability problem to the equations for a self-gravitating, hydrostatically stable white dwarf, with fixed overall composition. To our knowledge, this is the first attempt to incorporate complex multinary phases into the equilibrium phase layering diagram and mass-radius-composition dependence, both of which are reported for He-C-O and C-O-Ne white dwarfs. Finite thickness interfacial phases ("interphases") show up at the boundaries between single-component bcc crystalline regions, some of which have lower lattice symmetry than cubic. A second application -- quasi-static settling of heavy nuclei in white dwarfs -- builds on our equilibrium phase layering method. Tests of this nonequilibrium method reveal extra phases which play the role of transient host phases for the settling species.

Static and dynamic properties of low density outer envelopes of neutron stars are calculated within the nonlinear magnetic Thomas-Fermi model, assuming degenerate electrons. A novel domain decomposition enables proper description of lattice symmetry and may be seen as a prototype for the general class of problems involving nonlinear charge screening of periodic, quasi-low-dimensionality structures, e.g. liquid crystals. We describe a scalable implementation of the method using Hypre. Phase velocity of long wavelength transverse phonons is found to be a factor of 5-7 larger than in the corresponding Coulomb crystal model, which could have implications for low temperature phonon-mediated thermal conductivity. Other findings include $c'<0$ elastic instabilities for both bcc and fcc lattices, reminiscent of the situation in some light actinides, and suggestive of a symmetry-lowering transition to a tetragonal or orthorhombic lattice.