We perform general-relativistic simulations of charged black holes targeting GW150914. We show that the early inspiral is most efficient for detecting black hole charge through gravitational waves and that GW150914 places only weak constraints on the charge-to-mass ratio, ${Q}/{M}< 0.4$. Our work applies to electric and magnetic charge, and to theories with black holes endowed with U(1) (hidden or dark) charges. With our results, we place an upper bound on the deviation of the so-called theory of MOdified Gravity (MOG) from general relativity in the dynamical, strong-field regime.

A general relativistic, stationary and axisymmetric black hole in a four-dimensional asymptotically-flat spacetime is fully determined by its mass, angular momentum and electric charge. The expectation that astrophysically relevant black holes do not posses charge has resulted in a limited number of investigations of moving and charged black holes in the dynamical, strong-field gravitational (and electromagnetic) regime, where numerical studies are necessary. Apart from having a theoretical interest, the advent of multimessenger astronomy with gravitational waves offers new ways to think about charged black holes. In this work, we initiate an exploration of charged binary black holes by generating valid initial data for general relativistic simulations of black hole systems that have generic electric charge, linear and angular momenta. We develop our initial data formalism within the framework of the conformal transverse-traceless (Bowen-York) technique using the puncture approach, and apply the theory of isolated horizons to attribute physical parameters (mass, charge and angular momentum) to each hole. We implemented our formalism in the case of a binary system by modifying the publicly available TwoPunctures and QuasiLocalMeasures codes. We demonstrate that our code can recover existing solutions and that it has excellent self-convergence properties for a generic configuration of two black holes.

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.

A binary neutron star merger produces a rapidly and differentially rotating compact remnant whose lifespan heavily affects the electromagnetic and gravitational emissions. Its stability depends on both the equation of state (EOS) and the rotation law and it is usually investigated through numerical simulations. Nevertheless, by means of a sufficient criterion for secular instability, equilibrium sequences can be used as a computationally inexpensive way to estimate the onset of dynamical instability, which, in general, is close to the secular one. This method works well for uniform rotation and relies on the location of turning points: stellar models that are stationary points in a sequence of equilibrium solutions with constant rest mass or angular momentum. Here, we investigate differentially rotating models (using a large number of equations of state and different rotation laws) and find that several universal relations between properly scaled gravitational mass, rest mass and angular momentum of the turning-point models that are valid for uniform rotation, are insensitive to the degree of differential rotation, to high accuracy.