Can one distinguish a binary black hole undergoing a merger from a binary neutron star if the individual compact companions have masses that fall inside the so-called mass gap of $3-5\ M_\odot$? For neutron stars, achieving such masses typically requires extreme compactness and in this work we present initial data and evolutions of binary neutron stars initially in quasiequilibrium circular orbits having a compactness $C=0.336$. These are the most compact, nonvacuum, quasiequilibrium binary objects that have been constructed and evolved to date, including boson stars. The compactness achieved is only slightly smaller than the maximum possible imposed by causality, $C_{\rm max}=0.355$, which requires the sound speed to be less than the speed of light. By comparing the emitted gravitational waveforms from the late inspiral to merger and postmerger phases between such a binary neutron star vs a binary black hole of the same total mass we identify concrete measurements that serve to distinguish them. With that level of compactness, the binary neutron stars exhibit no tidal disruption up until merger, whereupon a prompt collapse is initiated even before a common core forms. Within the accuracy of our simulations the black hole remnants from both binaries exhibit ringdown radiation that is not distinguishable from a perturbed Kerr spacetime. However, their inspiral leads to phase differences of the order of $\sim 5$ rad over an $\sim 81$ km separation (1.7 orbits) while typical neutron stars exhibit phase differences of $\geq 20$ rad. Although a difference of $\sim 5$ rad can be measured by current gravitational wave laser interferometers (e.g., aLIGO/Virgo), uncertainties in the individual masses and spins will likely prevent distinguishing such compact, massive neutron stars from black holes.

We construct the first dynamically stable ergostars (equilibrium neutron stars that contain an ergoregion) for a compressible, causal equation of state. We demonstrate their stability by evolving both strict and perturbed equilibrium configurations in full general relativity for over a hundred dynamical timescales ($\gtrsim 30$ rotational periods) and observing their stationary behavior. This stability is in contrast to earlier models which prove radially unstable to collapse. Our solutions are highly differentially rotating hypermassive neutron stars with a corresponding spherical compaction of $C=0.3$. Such ergostars can provide new insights into the geometry of spacetimes around highly compact, rotating objects and on the equation of state at supranuclear densities. Ergostars may form as remnants of extreme binary neutron star mergers and possibly provide another mechanism for powering short gamma-ray bursts.

We perform long-term simulations of spinning binary neutron stars, with our highest dimensionless spin being $\chi \sim 0.32$. To assess the importance of spin during the inspiral we vary the spin, and also use two equations of state, one that consists of plain nuclear matter and produces compact stars (SLy), and a hybrid one that contains both nuclear and quark matter and leads to larger stars (ALF2). Using high resolution that has grid spacing $\Delta x\sim 98$ m on the finest refinement level, we find that the effects of spin in the phase evolution of a binary system can be larger than the one that comes from tidal forces. Our calculations demonstrate explicitly that although tidal effects are dominant for small spins ($\lesssim 0.1$), this is no longer true when the spins are larger, but still much smaller than the Keplerian limit.

We present a new initial data formulation to solve the full set of Einstein equations for spacetimes that contain a black hole under general conditions. The method can be used to construct complete initial data for spacetimes (the full metric) that contain a black hole. Contrary to most current studies the formulation requires minimal assumptions. For example, rather than imposing the form of the spatial conformal metric we impose 3 gauge conditions adapted to the coordinates describing the system under consideration. For stationary, axisymmetric spacetimes our method yields Kerr-Schild black holes in vacuum and rotating equilibrium neutron stars. We demonstrate the power of our new method by solving for the first time the whole system of Einstein equations for a nonaxisymmetric, self-gravitating torus in the presence of a black hole. The black hole has dimensionless spin $J_{\rm bh}/M_{\rm bh}^2=0.9918$, a rotation axis tilted at a $30^\circ$ angle with respect to the angular momentum of the disk, and a mass of $\sim 1/5$ of the disk.

Quasi-equilibrium models of uniformly rotating axisymmetric and triaxial quark stars are computed in general relativistic gravity scenario. The Isenberg-Wilson-Mathews (IWM) formulation is employed and the Compact Object CALculator (COCAL) code is extended to treat rotating stars with finite surface density and new equations of state (EOSs). Besides the MIT bag model for quark matter which is composed of de-confined quarks, we examine a new EOS proposed by Lai and Xu that is based on quark clustering and results in a stiff EOS that can support masses up to $3.3M_\odot$ in the case we considered. We perform convergence tests for our new code to evaluate the effect of finite surface density in the accuracy of our solutions and construct sequences of solutions for both small and high compactness. The onset of secular instability due to viscous dissipation is identified and possible implications are discussed. An estimate of the gravitational wave amplitude and luminosity based on quadrupole formulas is presented and comparison with neutron stars is discussed.

Targets for ground-based gravitational wave interferometers include continuous, quasiperiodic sources of gravitational radiation, such as isolated, spinning neutron stars. In this work we perform evolution simulations of uniformly rotating, triaxially deformed stars, the compressible analogues in general relativity of incompressible, Newtonian Jacobi ellipsoids. We investigate their stability and gravitational wave emission. We employ five models, both normal and supramassive, and track their evolution with different grid setups and resolutions, as well as with two different evolution codes. We find that all models are dynamically stable and produce a strain that is approximately one-tenth the average value of a merging binary system. We track their secular evolution and find that all our stars evolve towards axisymmetry, maintaining their uniform rotation, kinetic energy, and angular momentum profiles while losing their triaxiality.

Carter and Lichnerowicz have established that barotropic fluid flows are conformally geodesic and obey Hamilton's principle. This variational approach can accommodate neutral, or charged and poorly conducting, fluids. We show that, unlike what has been previously thought, this approach can also accommodate perfectly conducting magnetofluids, via the Bekenstein-Oron description of ideal magnetohydrodynamics. When Noether symmetries associated with Killing vectors or tensors are present in geodesic flows, they lead to constants of motion polynomial in the momenta. We generalize these concepts to hydrodynamic flows. Moreover, the Hamiltonian descriptions of ideal magnetohydrodynamics allow one to cast the evolution equations into a hyperbolic form useful for evolving rotating or binary compact objects with magnetic fields in numerical general relativity. Conserved circulation laws, such as those of Kelvin, Alfv\'en and Bekenstein-Oron, emerge simply as special cases of the Poincar\'e-Cartan integral invariant of Hamiltonian systems. We use this approach to obtain an extension of Kelvin's theorem to baroclinic (non-isentropic) fluids, based on a temperature-dependent time parameter. We further extend this result to perfectly or poorly conducting baroclinic magnetoflows. Finally, in the barotropic case, such magnetoflows are shown to also be geodesic, albeit in a Finsler (rather than Riemann) space.

We consider the influence of potential quark matter existing at high densities in neutron star interiors on gravitational waves (GW) emitted in a binary neutron star merger event. Two types of equations of state (EoS) at zero temperatures are used, one describing pure nuclear matter, the other nuclear matter with a phase transition to quark matter at very high densities. Binary equilibrium sequences close to the innermost stable circular orbit (ISCO) are calculated to determine the GW frequencies just before merger. It is found that EoS effects begin to play a role for gravitational masses larger than $M_\infty\simeq1.5M_\odot$. The difference in the gravitational wave frequency at the ISCO grows to up to $\simeq 10%$ for the maximal allowed mass given by the EoSs used. Then, we perform 3D hydrodynamic simulations for each EoS varying the initial mass and determine the characteristic GW frequencies of the merger remnants. The implications of quark matter show up mainly in a different collapse behaviour of the merger remnant. If the collapse does not take place immediately after merger, we find a phase difference between two EoS's in the post-merger GW signal. We also compare the GW frequencies emitted by the merger remnant to values from simulations using a polytropic EoS and find an imprint of the non-constant adiabatic index of our EoSs. All calculations are based on the conformally flat (CF) approximation to general relativity and the GW signal from the merger simulation is extracted up to quadrupole order.