Universality and intermittency in relativistic turbulent flows of a hot plasma

With the aim of determining the statistical properties of relativistic turbulence and unveiling novel and non-classical features, we resent the results of direct numerical simulations of driven turbulence in an ultrarelativistic hot plasma using high-order numerical schemes. We study the statistical properties of flows with average Mach number ranging from $\sim 0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find that flow quantities, such as the energy density or the local Lorentz factor, show large spatial variance even in the subsonic case as compressibility is enhanced by relativistic effects. The velocity field is highly intermittent, but its power-spectrum is found to be in good agreement with the predictions of the classical theory of Kolmogorov. Overall, our results indicate that relativistic effects are able to significantly enhance the intermittency of the flow and affect the high-order statistics of the velocity field, while leaving unchanged the low-order statistics, which instead appear to be universal and in good agreement with the classical Kolmogorov theory. To the best of our knowledge, these are the most accurate simulations of driven relativistic turbulence to date.Comment: 5 pages, 4 figures. Minor changes to match the version accepted on ApJ

Universality and intermittency in relativistic turbulent flows of a hot gas

With the aim of determining the statistical properties of relativistic turbulence and unveiling novel and non-classical features, we present the results of direct numerical simulations of driven turbulence in an ultrarelativistic hot plasma using high-order numerical schemes. We study the statistical properties of flows with average Mach number ranging from $\sim 0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find that flow quantities, such as the energy density or the local Lorentz factor, show large spatial variance even in the subsonic case as compressibility is enhanced by relativistic effects. The velocity field is highly intermittent, but its power-spectrum is found to be in good agreement with the predictions of the classical theory of Kolmogorov.Comment: Talk given at the ASTRONUM2012 conference on the 25th of June 201

Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes

We have developed the formalism necessary to employ the discontinuous-Galerkin approach in general-relativistic hydrodynamics. The formalism is firstly presented in a general 4-dimensional setting and then specialized to the case of spherical symmetry within a 3+1 splitting of spacetime. As a direct application, we have constructed a one-dimensional code, EDGES, which has been used to asses the viability of these methods via a series of tests involving highly relativistic flows in strong gravity. Our results show that discontinuous Galerkin methods are able not only to handle strong relativistic shock waves but, at the same time, to attain very high orders of accuracy and exponential convergence rates in smooth regions of the flow. Given these promising prospects and their affinity with a pseudospectral solution of the Einstein equations, discontinuous Galerkin methods could represent a new paradigm for the accurate numerical modelling in relativistic astrophysics.Comment: 24 pages, 19 figures. Small changes; matches version to appear in PR

Multimessenger Parameter Estimation of GW170817

We combine gravitational wave (GW) and electromagnetic (EM) data to perform a Bayesian parameter estimation of the binary neutron star (NS) merger GW170817. The EM likelihood is constructed from a fit to a large number of numerical relativity simulations which we combine with a lower bound on the mass of the remnant's accretion disk inferred from the modeling of the EM light curve. In comparison with previous works, our analysis yields a more precise determination of the tidal deformability of the binary, for which the EM data provide a lower bound, and of the mass ratio of the binary, with the EM data favoring a smaller mass asymmetry. The 90\% credible interval for the areal radius of a $1.4\ M_\odot$ NS is found to be $12.2^{+1.0}_{-0.8} \pm 0.2\ {\rm km}$ (statistical and systematic uncertainties).Comment: 7 pages, 3 figures, accepted to the EPJA Topical Issue: The first Neutron Star Merger Observation - Implications for Nuclear Physic

Binary Neutron Star Merger Simulations with a Calibrated Turbulence Model

Magnetohydrodynamic (MHD) turbulence in neutron star (NS) merger remnants can impact their evolution and multimessenger signatures, complicating the interpretation of present and future observations. Due to the high Reynolds numbers and the large computational costs of numerical relativity simulations, resolving all the relevant scales of the turbulence will be impossible for the foreseeable future. Here, we adopt a method to include subgrid-scale turbulence in moderate resolution simulations by extending the large-eddy simulation (LES) method to general relativity (GR). We calibrate our subgrid turbulence model with results from very-high-resolution GRMHD simulations, and we use it to perform NS merger simulations and study the impact of turbulence. We find that turbulence has a quantitative, but not qualitative impact on the evolution of NS merger remnants, on their gravitational wave signatures, and on the outflows generated in binary NS mergers. Our approach provides a viable path to quantify uncertainties due to turbulence in NS mergers.Comment: 20 pages, 6 figures. Submitted to the special issue "Numerical Relativity and Gravitational Wave" of Symmetry. Minor changes, matches accepted versio

High-Order Numerical-Relativity Simulations of Binary Neutron Stars

We report simulations of the inspiral and merger of binary neutron stars performed with \texttt{WhiskyTHC}, the first of a new generation of numerical relativity codes employing higher than second-order methods for both the spacetime and the hydrodynamic evolution. We find that the use of higher-order schemes improves substantially the quality of the gravitational waveforms extracted from the simulations when compared to those computed using traditional second-order schemes. The reduced de-phasing and the faster convergence rate allow us to estimate the phase evolution of the gravitational waves emitted, as well as the magnitude of finite-resolution effects, without the need of phase- or time-alignments or rescalings of the waves, as sometimes done in other works. Furthermore, by using an additional unpublished simulation at very high resolution, we confirm the robustness of our high convergence order of $3.2$.Comment: Submitted for the ASTRONUM-2014 proceedings. Includes a previously unpublished high-resolution simulatio

Critical Phenomena in Neutron Stars II: Head-on Collisions

We consider the head-on collision of equal-mass neutron stars boosted towards each other and we study the behavior of such systems near the threshold of black-hole formation. In particular, we confirm the previous findings by [1] that a type-I critical phenomenon can be observed by fine-tuning the initial mass of the two neutron stars. At the same time, we argue against the interpretation that the critical solution is not a perturbed spherical star and show instead that the metastable star corresponds to a (perturbed) equilibrium solution on the unstable branch of the equilibrium configurations. As a result, the head-on collision of two neutron stars near the critical threshold can be seen as a transition in the space of configurations from an initial stable solution over to a critical metastable one which can either migrate to a stable solution or collapse to a black hole. The critical exponent for this process shows a fine structure which was already observed in the case of the critical collapse of scalar fields but never before for perfect fluids.Comment: 15 pages, 19 figure