We consider the elastic constants of phases with nonspherical nuclei, so-called pasta phases, predicted to occur in the inner crust of a neutron star. First, we treat perfectly ordered phases and give numerical estimates for lasagna and spaghetti when the pasta elements are spatially uniform: the results are in order-of-magnitude agreement with the numerical simulations of Caplan, Schneider, and Horowitz, Phys. Rev. Lett. 121, 132701 (2018). We then turn to pasta phases without long-range order and calculate upper (Voigt) and lower (Reuss) bounds on the effective shear modulus and find that the lower bound is zero, but the upper bound is nonzero. To obtain better estimates, we then apply the self-consistent formalism and find that this predicts that the shear modulus of the phases without long-range order is zero if the pasta elements are spatially uniform. In numerical simulations, the pasta elements are found to be modulated spatially and we show that this modulation is crucial to obtaining a nonzero elastic moduli for pasta phases without long-range order. In the self-consistent formalism we find that, for lasagna, the effective shear modulus is linear in the elastic constants that do not vanish when the pasta elements are spatially uniform while, for spaghetti, it varies as the square root of these elastic constants. We also consider the behavior of the elastic constant associated with a homologous strain (hydrostatic compression) of the structure of the pasta phases without long-range order.

The interaction between nuclei in the inner crust of neutron stars consists of two contributions, the so-called "direct" interaction and an "induced" one due to density changes in the neutron fluid. For large nuclear separations $r$ the contributions from nuclear forces to each of these terms are shown to be nonzero. In the static limit they are equal in magnitude but have opposite signs and they cancel exactly. We analyze earlier results on effective interactions in the light of this finding. We consider the properties of long-wavelength collective modes and, in particular, calculate the degree of mixing between the lattice phonons and the phonons in the neutron superfluid. Using microscopic theory, we calculate the net non-Coulombic contribution to the nucleus--nucleus interaction and show that, for large $r$, the leading term is due to exchange of two phonons and varies as $1/r^7$: it is an analog of the Casimir--Polder interaction between neutral atoms.

Elastic properties of the solid regions of neutron star crusts and white dwarfs play an important role in theories of stellar oscillations. Matter in compact stars is presumably polycrystalline and, since the elastic properties of single crystals of such matter are very anisotropic, it is necessary to relate elastic properties of the polycrystal to those of a single crystal. We calculate the effective shear modulus of polycrystalline matter with randomly oriented crystallites using a self-consistent theory that has been very successful in applications to terrestrial materials and show that previous calculations overestimate the shear modulus by approximately 28%.

In the standard picture of the crust of a neutron star, matter there is simple: a body-centered-cubic (bcc) lattice of nuclei immersed in an essentially uniform electron gas. We show that at densities above that for neutron drip ($\sim4\times10^11$) g cm$^{-3}$ or roughly one thousandth of nuclear matter density, the interstitial neutrons give rise to an attractive interaction between nuclei that renders the lattice unstable. We argue that the likely equilibrium structure is similar to that in displacive ferroelectric materials such as BaTiO$_3$. As a consequence, properties of matter in the inner crust are expected to be much richer than previously appreciated and we mention consequences for observable neutron star properties.