We investigate the dynamical stability of relativistic, differentially rotating, quasi-toroidal models of neutron stars through hydrodynamical simulations in full general relativity. We find that all quasi-toroidal configurations studied in this work are dynamically unstable against the growth of non-axisymmetric modes. Both one-arm and bar mode instabilities grow during their evolution. We find that very high rest mass configurations collapse to form black holes. Our calculations suggest that configurations whose rest mass is less than the binary neutron star threshold mass for prompt collapse to black hole transition dynamically to spheroidal, differentially rotating stars that are dynamically stable, but secularly unstable. Our study shows that the existence of extreme quasi-toroidal neutron star equilibrium solutions does not imply that long-lived binary neutron star merger remnants can be much more massive than previously found. Finally, we find models that are initially supra-Kerr ($J/M^2>1$) and undergo catastrophic collapse on a dynamical timescale, in contrast to what was found in earlier works. However, cosmic censorship is respected in all of our cases. Our work explicitly demonstrates that exceeding the Kerr bound in rotating neutron star models does not imply dynamical stability.

We construct equilibrium models of uniformly and differentially rotating hybrid hadron-quark stars using equations of state (EOSs) with a first-order phase transition that gives rise to a third family of compact objects. We find that the ratio of the maximum possible mass of uniformly rotating configurations - the supramassive limit - to the Tolman-Oppenheimer-Volkoff (TOV) limit mass is not EOS-independent, and is between 1.15 and 1.31,in contrast with the value of 1.20 previously found for hadronic EOSs. Therefore, some of the constraints placed on the EOS from the observation of the gravitational wave event GW170817 do not apply to hadron-quark EOSs. However, the supramassive limit mass for the family of EOSs we treat is consistent with limits set by GW170817, strengthening the possibility of interpreting GW170817 with a hybrid hadron-quark EOSs. We also find that along constant angular momentum sequences of uniformly rotating stars, the third family maximum and minimum mass models satisfy approximate EOS-independent relations, and the supramassive limit of the third family is approximately 16.5 % larger than the third family TOV limit. For differentially rotating spheroidal stars, we find that a lower-limit on the maximum supportable rest mass is 123 % more than the TOV limit rest mass. Finally, we verify that the recently discovered universal relations relating angular momentum, rest mass and gravitational mass for turning-point models hold for hybrid hadron-quark EOSs when uniform rotation is considered, but have a clear dependence on the degree of differential rotation.

We study the solution space of general relativistic, axisymmetric, equilibria of differentially rotating neutron stars with realistic, nuclear equations of state. We find that different types of stars, which were identified by earlier works for polytropic equations of state, arise for realistic equations of state, too. Scanning the solution space for the sample of realistic equations of state we treat, we find lower limits on the maximum rest masses supported by cold, differentially rotating stars for each type of stars. We often discover equilibrium configurations that can support more than 2 times the mass of a static star. We call these equilibria "\"ubermassive", and in our survey we find \"ubermassive stars that can support up to 2.5 times the maximum rest mass that can be supported by a cold, non-rotating star with the same equation of state. This is nearly two times larger than what previous studies employing realistic equations of state had found, and which did not uncover \"ubermassive neutron stars. Moreover, we find that the increase in the maximum rest mass with respect to the non-spinning stellar counterpart is larger for moderately stiff equations of state. These results may have implications for the lifetime and the gravitational wave and electromagnetic counterparts of hypermassive neutron stars formed following binary neutron star mergers.