Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large
variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes.
With the recent advent of multi-messenger gravitational-wave astronomy and in anticipation of more sensitive, next-generation gravitational-wave detectors, we investigate the dynamics, gravitational-wave emission, and nucleosynthetic yields of numero
us eccentric binary neutron-star mergers having different equations of state. For each equation of state we vary the orbital properties around the threshold of immediate merger, as well as the binary mass ratio. In addition to a study of the gravitational-wave emission including $f$-mode oscillations before and after merger, we couple the dynamical ejecta output from the simulations to the nuclear-reaction network code texttt{SkyNet} to compute nucleosynthetic yields and compare to the corresponding results in the case of a quasi-circular merger. We find that the amount and velocity of dynamically ejected material is always much larger than in the quasi-circular case, reaching maximal values of $M_{rm ej, max} sim 0.1 , M_{odot}$ and $v_{rm max}/c sim 0.75$. At the same time, the properties of this material are rather insensitive to the details of the orbit, such as pericenter distance or post-encounter apoastron distance. Furthermore, while the composition of the ejected matter depends on the orbital parameters and on the equation of state, the relative nucleosynthetic yields do not, thus indicating that kilonova signatures could provide information on the orbital properties of dynamically captured neutron-star binaries.
We explore in a parameterized manner a very large range of physically plausible equations of state (EOSs) for compact stars for matter that is either purely hadronic or that exhibits a phase transition. In particular, we produce two classes of EOSs w
ith and without phase transitions, each containing one million EOSs. We then impose constraints on the maximum mass, ($M < 2.16 M_{odot}$), and on the dimensionless tidal deformability ($tilde{Lambda} <800$) deduced from GW170817, together with recent suggestions of lower limits on $tilde{Lambda}$. Exploiting more than $10^9$ equilibrium models for each class of EOSs, we produce distribution functions of all the stellar properties and determine, among other quantities, the radius that is statistically most probable for any value of the stellar mass. In this way, we deduce that the radius of a purely hadronic neutron star with a representative mass of $1.4,M_{odot}$ is constrained to be $12.00!<!R_{1.4}/{rm km}!<!13.45$ at a $2$-$sigma$ confidence level, with a most likely value of $bar{R}_{1.4}=12.39,{rm km}$; similarly, the smallest dimensionless tidal deformability is $tilde{Lambda}_{1.4}!>!375$, again at a $2$-$sigma$ level. On the other hand, because EOSs with a phase transition allow for very compact stars on the so-called `twin-star branch, small radii are possible with such EOSs although not probable, i.e. $8.53!<!R_{1.4}/{rm km}!<!13.74$ and $bar{R}_{1.4}=13.06,{rm km}$ at a $2$-$sigma$ level, with $tilde{Lambda}_{1.4}!>!35.5$ at a $3$-$sigma$ level. Finally, since these EOSs exhibit upper limits on $tilde{Lambda}$, the detection of a binary with total mass of $3.4,M_{odot}$ and $tilde{Lambda}_{1.7}!>!461$ can rule out twin-star solutions.
We investigate the effect of large magnetic fields on the $2+1$ dimensional reduced-magnetohydrodynamical expansion of hot and dense nuclear matter produced in $sqrt{s_{rm NN}}$ = 200 GeV Au+Au collisions. For the sake of simplicity, we consider the
case where the magnetic field points in the direction perpendicular to the reaction plane. We also consider this field to be external, with energy density parametrized as a two-dimensional Gaussian. The width of the Gaussian along the directions orthogonal to the beam axis varies with the centrality of the collision. The dependence of the magnetic field on proper time ($tau$) for the case of zero electrical conductivity of the QGP is parametrized following [Deng 2012], and for finite electrical conductivity following [Tuchin 2013]. We solve the equations of motion of ideal hydrodynamics for such an external magnetic field. For collisions with non-zero impact parameter we observe considerable changes in the evolution of the momentum eccentricities of the fireball when comparing the case when the magnetic field decays in a conducting QGP medium and when no magnetic field is present. The elliptic-flow coefficient $v_2$ of $pi^{-}$ is shown to increase in the presence of an external magnetic field and the increment in $v_2$ is found to depend on the evolution and the initial magnitude of the magnetic field.
We show how gravitational-wave observations with advanced detectors of tens to several tens of neutron-star binaries can measure the neutron-star radius with an accuracy of several to a few percent, for mass and spatial distributions that are realist
ic, and with none of the sources located within 100 Mpc. We achieve such an accuracy by combining measurements of the total mass from the inspiral phase with those of the compactness from the postmerger oscillation frequencies. For estimating the measurement errors of these frequencies we utilize analytical fits to postmerger numerical-relativity waveforms in the time domain, obtained here for the first time, for four nuclear-physics equations of state and a couple of values for the mass. We further exploit quasi-universal relations to derive errors in compactness from those frequencies. Measuring the average radius to well within 10% is possible for a sample of 100 binaries distributed uniformly in volume between 100 and 300 Mpc, so long as the equation of state is not too soft or the binaries are not too heavy.
A numerical-relativity calculation yields in general a solution of the Einstein equations including also a radiative part, which is in practice computed in a region of finite extent. Since gravitational radiation is properly defined only at null infi
nity and in an appropriate coordinate system, the accurate estimation of the emitted gravitational waves represents an old and non-trivial problem in numerical relativity. A number of methods have been developed over the years to extract the radiative part of the solution from a numerical simulation and these include: quadrupole formulas, gauge-invariant metric perturbations, Weyl scalars, and characteristic extraction. We review and discuss each method, in terms of both its theoretical background as well as its implementation. Finally, we provide a brief comparison of the various methods in terms of their inherent advantages and disadvantages.
A number of works have shown that important information on the equation of state of matter at nuclear density can be extracted from the gravitational waves emitted by merging neutron-star binaries. We present a comprehensive analysis of the gravitati
onal-wave signal emitted during the inspiral, merger and post-merger of 56 neutron-star binaries. This sample of binaries, arguably the largest studied to date with realistic equations of state, spans across six different nuclear-physics equations of state and ten masses, allowing us to sharpen a number of results recently obtained on the spectral properties of the gravitational-wave signal. Overall we find that: (i) for binaries with masses differing no more than $20%$, the frequency at gravitational-wave amplitudes maximum is related quasi-universally with the tidal deformability of the two stars; (ii) the spectral properties vary during the post-merger phase, with a transient phase lasting a few millisecond after the merger and followed by a quasi-stationary phase; (iii) when distinguishing the spectral peaks between these two phases, a number of ambiguities in the identification of the peaks disappear, leaving a simple and robust picture; (iv) using properly identified frequencies, quasi-universal relations are found between the spectral features and the properties of the neutron stars; (v) for the most salient peaks analytic fitting functions can be obtained in terms of the stellar tidal deformability or compactness. Altogether, these results support the idea that the equation of state of nuclear matter can be constrained tightly when a signal in gravitational waves from binary neutron stars is detected.
In the initial stage of relativistic heavy-ion collisions, strong magnetic fields appear due to the large velocity of the colliding charges. The evolution of these fields appears as a novel and intriguing feature in the fluid-dynamical description of
heavy-ion collisions. In this work, we study analytically the one-dimensional, longitudinally boost-invariant motion of an ideal fluid in the presence of a transverse magnetic field. Interestingly, we find that, in the limit of ideal magnetohydrodynamics, i.e., for infinite conductivity, and irrespective of the strength of the initial magnetization, the decay of the fluid energy density $e$ with proper time $tau$ is the same as for the time-honored Bjorken flow without magnetic field. Furthermore, when the magnetic field is assumed to decay $sim tau^{-a}$, where $a$ is an arbitrary number, two classes of analytic solutions can be found depending on whether $a$ is larger or smaller than one. In summary, the analytic solutions presented here highlight that the Bjorken flow is far more general than formerly thought. These solutions can serve both to gain insight on the dynamics of heavy-ion collisions in the presence of strong magnetic fields and as testbeds for numerical codes.
Extending previous work by a number of authors, we have recently presented a new approach in which the detection of gravitational waves from merging neutron star binaries can be used to determine the equation of state of matter at nuclear density and
hence the structure of neutron stars. In particular, after performing a large number of numerical-relativity simulations of binaries with nuclear equations of state, we have found that the post-merger emission is characterized by two distinct and robust spectral features. While the high-frequency peak was already shown to be associated with the oscillations of the hypermassive neutron star produced by the merger and to depend on the equation of state, we have highlighted that the low-frequency peak is related to the merger process and to the total compactness of the stars in the binary. This relation is essentially universal and provides a powerful tool to set tight constraints on the equation of state. We here provide additional information on the extensive analysis performed, illustrating the methods used, the tests considered, as well as the robustness of the results. We also discuss additional relations that can be deduced when exploring the data and how these correlate with various properties of the binary. Finally, we present a simple mechanical toy model that explains the main spectral features of the post-merger signal and can even reproduce analytically the complex waveforms emitted right after the merger.
Determining the equation of state of matter at nuclear density and hence the structure of neutron stars has been a riddle for decades. We show how the imminent detection of gravitational waves from merging neutron star binaries can be used to solve t
his riddle. Using a large number of accurate numerical-relativity simulations of binaries with nuclear equations of state, we find that the postmerger emission is characterized by two distinct and robust spectral features. While the high-frequency peak has already been associated with the oscillations of the hypermassive neutron star produced by the merger and depends on the equation of state, a new correlation emerges between the low-frequency peak, related to the merger process, and the total compactness of the stars in the binary. More importantly, such a correlation is essentially universal, thus providing a powerful tool to set tight constraints on the equation of state. If the mass of the binary is known from the inspiral signal, the combined use of the two frequency peaks sets four simultaneous constraints to be satisfied. Ideally, even a single detection would be sufficient to select one equation of state over the others. We test our approach with simulated data and verify it works well for all the equations of state considered.
