Seismology of adolescent neutron stars: Accounting for thermal effects and crust elasticity

We study the oscillations of relativistic stars, incorporating key physics associated with internal composition, thermal gradients and crust elasticity. Our aim is to develop a formalism which is able to account for the state-of-the-art understanding of the complex physics associated with these systems. As a first step, we build models using a modern equation of state including composition gradients and density discontinuities associated with internal phase-transitions (like the crust-core transition and the point where muons first appear in the core). In order to understand the nature of the oscillation spectrum, we carry out cooling simulations to provide realistic snapshots of the temperature distribution in the interior as the star evolves through adolescence. The associated thermal pressure is incorporated in the perturbation analysis, and we discuss the presence of $g$-modes arising as a result of thermal effects. We also consider interface modes due to phase-transitions and the gradual formation of the star's crust and the emergence of a set of shear modes.