Structure of multicomponent Coulomb crystals

Coulomb plasmas crystallize in a number of physical systems, such as dusty plasmas, neutron star crusts, and white dwarf cores. The crystal structure of the one-component and binary plasma has received significant attention in the literature, though the less studied multicomponent plasma may be most relevant for many physical systems which contain a large range of particle charges. We report on molecular dynamics simulations of multicomponent plasmas near the melting temperature with mixtures taken to be realistic x-ray burst ash compositions. We quantify the structure of the crystal with the bond order parameters and radial distribution function. Consistent with past work, low charge particles form interstitial defects and we argue that they are in a quasiliquid state within the lattice. The lattice shows screening effects which preserves long-range order despite the large variance in particle charges, which may impact transport properties relevant to astrophysics.