Couplings of torsional and shear oscillations in a neutron star crust

Mature neutron stars are thought to be sufficiently cold that nuclei in the outer layers freeze, solidifying a crust. Crustal elasticity allows the star to support a set of seismic modes, such as torsional oscillations. These axial-parity modes can couple to the polar sector in a number of ways, for example via rotation or a magnetic field. Even in a static, spherically symmetric star, however, these modes can couple at nonlinear order. In this study, such couplings in the crust are examined for the first time: we derive the axisymmetric perturbation equations for second-order axial eigenfunctions, which are sourced by axial-polar couplings at first order, and solve the resulting equations in the time domain. Through our studies, we find that the second-order spectrum contains additional oscillation modes, not predicted by the linear analysis with either axial or polar perturbations, which can be excited to relatively large amplitudes.