Suppression of composition $g$-modes in chemically-equilibrating warm neutron stars

We investigate the impact of chemical equilibration and the resulting bulk viscosity on non-radial oscillation modes of warm neutron stars at temperatures up to $T\approx 5$ MeV, relevant for protoneutron stars and neutron-star post-merger remnants. In this regime, the relaxation rate of weak interactions becomes comparable to the characteristic frequencies of composition $g$-modes in the core, resulting in resonant damping. To capture this effect, we introduce the dynamical sound speed, a complex, frequency-dependent generalization of the adiabatic sound speed that encodes both the restoring force and the dissipative effects of bulk compression. Using realistic weak reaction rates and three representative equations of state, we compute the complex frequencies of composition $g$-modes with finite-temperature profiles. We find that bulk viscous damping becomes increasingly significant with temperature and can completely suppress composition $g$-modes. In contrast, the $f$-mode remains largely unaffected by bulk viscosity due to its nearly divergence-free character. Our results highlight the sensitivity of $g$-mode behavior to thermal structure, weak reaction rates, and the equation of state, and establish the dynamical sound speed as a valuable descriptor characterizing oscillation properties in dissipative neutron star matter.