We consider the local dynamics of a realistic neutron star core, including
composition gradients, superfluidity and thermal effects. The main focus is on
the gravity g-modes, which are supported by composition stratification and
thermal gradients. We derive the equations that govern this problem in full
detail, paying particular attention to the input that needs to be provided
through the equation of state and distinguishing between normal and superfluid
regions. The analysis highlights a number of key issues that should be kept in
mind whenever equation of state data is compiled from nuclear physics for use
in neutron star calculations. We provide explicit results for a particular
stellar model and a specific nucleonic equation of state, making use of cooling
simulations to show how the local wave spectrum evolves as the star ages. Our
results show that the composition gradient is effectively dominated by the
muons whenever they are present. When the star cools below the superfluid
transition, the support for g-modes at lower densities (where there are no
muons) is entirely thermal. We confirm the recent suggestion that the g-modes
in this region may be unstable, but our results indicate that this instability
will be weak and would only be present for a brief period of the star's life.
Our analysis accounts for the presence of thermal excitations encoded in
entrainment between the entropy and the superfluid component. Finally, we
discuss the complete spectrum, including the normal sound waves and, in
superfluid regions, the second sound.Comment: 29 pages, 9 figures, submitted to MNRA
