Elastic properties play an important role in neutron stars and white dwarfs. They are crucial for modeling stellar oscillations and different processes in magnetars and in degenerate stars which enter compact binary systems. Using electrostatic energy of deformed lattices, we calculate elastic moduli of ordered binary body-centered cubic (sc2) and face-centered cubic (fc2) lattices. We use two methods to determine the effective shear modulus $\mu_{\rm eff}$. We show that $\mu_{\rm eff}$ calculated as an average of the shear stiffness over all possible wavevectors (for polarization vectors perpendicular to the wave vector) agrees with the results obtained from the linear mixing rule. For the sc2 lattice, our calculations are also consistent with the results of numerical simulations of disordered binary body-centered cubic lattice.

We have studied phonon and thermodynamic properties of a body-centered cubic (bcc) Coulomb crystal of ions with weakly polarized electron background in a uniform magnetic field ${\bf B}$. At $B=0$ the difference between phonon moments calculated using the Thomas-Fermi (TF) and random phase (RPA) approximations is always less than one percent and for description of phonon properties of a crystal TF formalism was used. This formalism was successfully applied to investigate thermodynamic properties of magnetized Coulomb crystals. It was shown that the influence of the polarization of the electron background is significant only at $\kappa_{\rm TF}a > 0.1$ and $T \ll T_{\rm p}(1+h^2)^{-1/2}$, where $\kappa_{\rm TF}$ is the Thomas-Fermi wavenumber, $a$ is the ion sphere radius, $T_{\rm p}\equiv\hbar \omega_{\rm p}$ is the ion plasma temperature, $h\equiv \omega_B/\omega_{\rm p}$, $\omega_B$ is the ion cyclotron frequency and $\omega_{\rm p}$ is the ion plasma frequency.

We study electrostatic and phonon properties of Yukawa crystals. It is shown that in the harmonic approximation these systems which is in use in the theory dusty plasma can be described analytically by the model from the theory of neutron stars and white dwarfs. Using this approximation we consider properties of body-centred cubic (bcc), face-centred cubic (fcc), hexagonal close-packed (hcp), and MgB$_2$ lattices. MgB$_2$ and hcp lattices have never been studied earlier in the context of Yukawa systems. It is shown that they never possess the smallest potential energy and the phase diagram of stable Yukawa crystals contains bcc and fcc lattices only. However, corrections to the charge density $\propto (\kappa a)^4$ can noticeably change the structural diagram of Yukawa systems. The developed analytical model also allows to describe low-temperature effects where numerical simulations are difficult.

Although crystallized neutron star crust is responsible for many fascinating observational phenomena, its actual microscopic structure in tremendous gravitational and magnetic fields is not understood. Here we show that in a non-uniform magnetic field, three-dimensional ionic Coulomb crystals comprising the crust may stretch or shrink while their electrostatic pressure becomes anisotropic. The pressure depends non-linearly on the magnitude of the stretch, so that a continuous magnetic field evolution may result in an abrupt crystal elongation or contraction. This may provide a trigger for magnetar activity. A phonon mode instability is revealed, which sets the limits of magnetic field variation beyond which the crystal is destroyed. These limits sometimes correspond to surprisingly large deformations. It is not known what happens to crust matter subject to a pressure anisotropy exceeding these limits. We hypothesize that the ion system then possesses a long-range order only in one or two dimensions, that is becomes a liquid crystal.