We continue our study of the binary neutron star parameter space by investigating the effect of the spin orientation on the dynamics, gravitational wave emission, and mass ejection during the binary neutron star coalescence. We simulate seven different configurations using multiple resolutions to allow a reasonable error assessment. Due to the particular choice of the setups, five configurations show precession effects, from which two show a precession (`wobbling') of the orbital plane, while three show a `bobbing' motion, i.e., the orbital angular momentum does not precess, while the orbital plane moves along the orbital angular momentum axis. Considering the ejection of mass, we find that precessing systems can have an anisotropic mass ejection, which could lead to a final remnant kick of about $\sim 40 \rm km/s$ for the studied systems. Furthermore, for the chosen configurations, anti-aligned spins lead to larger mass ejecta than aligned spins, so that brighter electromagnetic counterparts could be expected for these configurations. Finally, we compare our simulations with the precessing, tidal waveform approximant IMRPhenomPv2_NRTidalv2 and find good agreement between the approximant and our numerical relativity waveforms with phase differences below 1.2 rad accumulated over the last $\sim$ 16 gravitational wave cycles.

The gravitational wave and electromagnetic signatures connected to the merger of two neutron stars allow us to test the nature of matter at supranuclear densities. Since the Equation of State governing the interior of neutron stars is only loosely constrained, there is even the possibility that strange quark matter exists inside the core of neutron stars. We investigate how strange quark matter cores affect the binary neutron star coalescence by performing numerical relativity simulations. Interestingly, the strong phase transition can cause a reduction of the convergence order of the numerical schemes to first order if the numerical resolution is not high enough. Therefore, an additional challenge is added in producing high-quality gravitational wave templates for Equation of States with a strong phase transition. Focusing on one particular configuration of an equal mass configuration consistent with GW170817, we compute and discuss the associated gravitational wave signal and some of the electromagnetic counterparts connected to the merger of the two stars. We find that existing waveform approximants employed for the analysis of GW170817 allow describing this kind of systems within the numerical uncertainties, which, however, are several times larger than for pure hadronic Equation of States, which means that even higher resolutions have been employed for an accurate gravitational wave model comparison. We also show that for the chosen Equation of State, quasi-universal relations describing the gravitational wave emission after the moment of merger seem to hold and that the electromagnetic signatures connected to our chosen setup would not be bright enough to explain the kilonova associated to GW170817.

As current gravitational wave (GW) detectors increase in sensitivity, and particularly as new instruments are being planned, there is the possibility that ground-based GW detectors will observe GWs from highly eccentric neutron star binaries. We present the first detailed study of highly eccentric BNS systems with full (3+1)D numerical relativity simulations using consistent initial conditions, i.e., setups which are in agreement with the Einstein equations and with the equations of general relativistic hydrodynamics in equilibrium. Overall, our simulations cover two different equations of state (EOSs), two different spin configurations, and three to four different initial eccentricities for each pairing of EOS and spin. We extract from the simulated waveforms the frequency of the f-mode oscillations induced during close encounters before the merger of the two stars. The extracted frequency is in good agreement with f-mode oscillations of individual stars for the irrotational cases, which allows an independent measure of the supranuclear equation of state not accessible for binaries on quasi-circular orbits. The energy stored in these f-mode oscillations can be as large as $10^{-3}M_\odot \sim 10^{51}$ erg, even with a soft EOS. In order to estimate the stored energy, we also examine the effects of mode mixing due to the stars' offset from the origin on the f-mode contribution to the GW signal. While in general (eccentric) neutron star mergers produce bright electromagnetic counterparts, we find that the luminosity decreases when the eccentricity becomes too large, due to a decrease of the ejecta mass. Finally, the use of consistent initial configurations also allows us to produce high-quality waveforms for different eccentricities which can be used as a testbed for waveform model development of highly eccentric binary neutron star systems.

We present the Computational Relativity CoRe collaboration's public database of gravitational waveforms from binary neutron star mergers. The database currently contains 367 waveforms from numerical simulations that are consistent with general relativity and that employ constraint satisfying initial data in hydrodynamical equilibrium. It spans 164 physically distinct configuration with different binary parameters (total binary mass, mass-ratio, initial separation, eccentricity, and stars' spins) and simulated physics. Waveforms computed at multiple grid resolutions and extraction radii are provided for controlling numerical uncertainties. We also release an exemplary set of 18 hybrid waveforms constructed with a state-of-art effective-one-body model spanning the frequency band of advanced gravitational-wave detectors. We outline present and future applications of the database to gravitational-wave astronomy.

We present the first set of numerical relativity simulations of binary neutron mergers that include spin precession effects and are evolved with multiple resolutions. Our simulations employ consistent initial data in general relativity with different spin configurations and dimensionless spin magnitudes $\sim 0.1$. They start at a gravitational-wave frequency of $\sim392$~Hz and cover more than $1$ precession period and about 15 orbits up to merger. We discuss the spin precession dynamics by analyzing coordinate trajectories, quasi-local spin measurements, and energetics, by comparing spin aligned, antialigned, and irrotational configurations. Gravitational waveforms from different spin configuration are compared by calculating the mismatch between pairs of waveforms in the late inspiral. We find that precession effects are not distinguishable from nonprecessing configurations with aligned spins for approximately face-on binaries, while the latter are distinguishable from a nonspinning configurations. Spin precession effects are instead clearly visible for approximately edge-on binaries. For the parameters considered here, precession does not significantly affect the characteristic postmerger gravitational-wave frequencies nor the mass ejection. Our results pave the way for the modeling of spin precession effects in the gravitational waveform from binary neutron star events.

With the increasing sensitivity of advanced gravitational wave detectors, the first joint detection of an electromagnetic and gravitational wave signal from a compact binary merger will hopefully happen within this decade. However, current gravitational-wave likelihood sky areas span $\sim 100-1000\,\textrm{deg}^2$, and thus it is a challenging task to identify which, if any, transient corresponds to the gravitational-wave event. In this study, we make a comparison between recent kilonovae/macronovae lightcurve models for the purpose of assessing potential lightcurve templates for counterpart identification. We show that recent analytical and parametrized models for these counterparts result in qualitative agreement with more complicated radiative transfer simulations. Our analysis suggests that with improved lightcurve models with smaller uncertainties, it will become possible to extract information about ejecta properties and binary parameters directly from the lightcurve measurement. Even tighter constraints are obtained in cases for which gravitational-wave and kilonovae parameter estimation results are combined. However, to be prepared for upcoming detections, more realistic kilonovae models are needed. These will require numerical relativity with more detailed microphysics, better radiative transfer simulations, and a better understanding of the underlying nuclear physics.

In addition to the emission of gravitational waves (GWs) the coalescence and merger of two neutron stars will produce a variety of electromagnetic (EM) signals. In this work we combine a large set of numerical relativity simulations performed by different groups and we present fits for the mass, kinetic energy, and the velocities of the dynamical ejected material. Additionally, we comment on the geometry and composition of the ejecta and discuss the influence of the stars' individual rotation. The derived fits can be used to approximate the luminosity and lightcurve of the kilonovae (macronovae) and to estimate the main properties of the radio flares. This correlation between the binary parameters and the EM signals allows in case of a GW detection to approximate possible EM counterparts when first estimates of the masses are available. After a possible kilonovae observation our results could also be used to restrict the region of the parameter space which has to be covered by numerical relativity simulations.

We present new (3+1) dimensional numerical relativity simulations of the binary neutron star (BNS) mergers that take into account the NS spins. We consider different spin configurations, aligned or antialigned to the orbital angular momentum, for equal and unequal mass BNS and for two equations of state. All the simulations employ quasiequilibrium circular initial data in the constant rotational velocity approach, i.e. they are consistent with Einstein equations and in hydrodynamical equilibrium. We study the NS rotation effect on the energetics, the gravitational waves (GWs) and on the possible electromagnetic (EM) emission associated to dynamical mass ejecta. For dimensionless spin magnitudes of $\chi\sim0.1$ we find that spin-orbit interactions and also spin-induced-quadrupole deformations affect the late-inspiral-merger dynamics. The latter is, however, dominated by finite-size effects. Spin (tidal) effects contribute to GW phase differences up to 5 (20) radians accumulated during the last eight orbits to merger. Similarly, after merger the collapse time of the remnant and the GW spectrogram are affected by the NSs rotation. Spin effects in dynamical ejecta are clearly observed in unequal mass systems in which mass ejection originates from the tidal tail of the companion. Consequently kilonovae and other EM counterparts are affected by spins. We find that spin aligned to the orbital angular momentum leads to brighter EM counterparts than antialigned spin with luminosities up to a factor of two higher.

We present new (3+1)D numerical relativity simulations of the binary neutron star (BNS) merger and postmerger phase. We focus on a previously inaccessible region of the binary parameter space spanning the binary's mass-ratio $q\sim1.00-1.75$ for different total masses and equations of state, and up to $q\sim2$ for a stiff BNS system. We study the mass-ratio effect on the gravitational waves (GWs) and on the possible electromagnetic emission associated to dynamical mass ejecta. We compute waveforms, spectra, and spectrograms of the GW strain including all the multipoles up to $l=4$. The mass-ratio has a specific imprint on the GW multipoles in the late-inspiral-merger signal, and it affects qualitatively the spectra of the merger remnant. The multipole effect is also studied by considering the dependency of the GW spectrograms on the source's sky location. Unequal mass BNSs produce more ejecta than equal mass systems with ejecta masses and kinetic energies depending almost linearly on $q$. We estimate luminosity peaks and light curves of macronovae events associated to the mergers using a simple approach. For $q\sim2$ the luminosity peak is delayed for several days and can be up to four times larger than for the $q=1$ cases. The macronova emission associated with the $q\sim2$ BNS is more persistent in time and could be observed for weeks instead of few days ($q=1$) in the near infrared. Finally, we estimate the flux of possible radio flares produced by the interaction of relativistic outflows with the surrounding medium. Also in this case a large $q$ can significantly enhance the emission and delay the peak luminosity. Overall, our results indicate that BNS merger with large mass ratio have EM signatures distinct from the equal mass case and more similar to black hole - neutron star binaries.

We constrain scattering length parameters in a Bose-Einstein condensate dark matter model by using galaxy clusters radii, with the implementation of a method previously applied to galaxies. At the present work, we use a sample of 114 clusters radii in order to obtain the scattering lengths associated with a dark matter particle mass in the range $10^{-6}-10^{-4}\, {\rm eV}$. We obtain scattering lengths that are five orders of magnitude larger than the ones found in the galactic case, even when taking into account the cosmological expansion in the cluster scale by means of the introduction of a small cosmological constant. We also construct and compare curves for the orbital velocity of a test particle in the vicinity of a dark matter cluster in both the expanding and the non-expanding cases.

We study equal and unequal-mass neutron star mergers by means of new numerical relativity simulations in which the general relativistic hydrodynamics solver employs an algorithm that guarantees mass conservation across the refinement levels of the computational mesh. We consider eight binary configurations with total mass $M=2.7\,M_\odot$, mass-ratios $q=1$ and $q=1.16$, and four different equation of states (EOSs), and one configuration with a stiff EOS, $M=2.5M_\odot$ and $q=1.5$. We focus on the post-merger dynamics and study the merger remnant, dynamical ejecta and the postmerger gravitational wave spectrum. Although most of the merger remnants form a hypermassive neutron star collapsing to a black hole+disk system on dynamical timescales, stiff EOSs can eventually produce a stable massive neutron star. Ejecta are mostly emitted around the orbital plane; favored by large mass ratios and softer EOS. The postmerger wave spectrum is mainly characterized by non-axisymmetric oscillations of the remnant. The stiff EOS configuration consisting of a $1.5M_\odot$ and a $1.0M_\odot$ neutron star shows a rather peculiar dynamics. During merger the companion star is very deformed; about~$\sim0.03M_\odot$ of rest-mass becomes unbound from the tidal tail due torque; and the merger remnant forms stable neutron star surrounded by a massive accretion disk $\sim0.3M_\odot$. Similar configurations might be particularly interesting for electromagnetic counterparts. Comparing results obtained with and without the conservative mesh refinement algorithm, we find that post-merger simulations can be affected by systematic errors if mass conservation is not enforced in the mesh refinement strategy. However, mass conservation also depends on grid details and on the artificial atmosphere setup. [abridged]

Binary neutron star mergers are studied using nonlinear 3+1 numerical relativity simulations and the analytical effective-one-body (EOB) model. The EOB model predicts quasiuniversal relations between the mass-rescaled gravitational wave frequency and the binding energy at the moment of merger, and certain dimensionless binary tidal coupling constants depending on the stars Love numbers, compactnesses and the binary mass ratio. These relations are quasiuniversal in the sense that, for a given value of the tidal coupling constant, they depend significantly neither on the equation of state nor on the mass ratio, though they do depend on stars spins. The spin dependence is approximately linear for small spins aligned with the orbital angular momentum. The quasiuniversality is a property of the conservative dynamics; nontrivial relations emerge as the binary interaction becomes tidally dominated. This analytical prediction is qualitatively consistent with new, multi-orbit numerical relativity results for the relevant case of equal-mass irrotational binaries. Universal relations are thus expected to characterize neutron star mergers dynamics. In the context of gravitational wave astronomy, these universal relations may be used to constrain the neutron star equation of state using waveforms that model the merger accurately.

An analytical solution of Einstein-Maxwell equations with a static fluid as a source is presented. The spacetime is represented by the axially symmetric Weyl metric and the energy-momentum tensor describes a coupling of a fluid with an electromagnetic field. When appropriate limits are performed we recover the well-known solutions of Gutsunaev-Manko and Schwarzschild. Also, using Eckart's thermodynamics, we calculated the temperature, the mechanical pressure, the charge density and the energy density of the system. The analysis of thermodynamic quantities suggests that the solution can be used to represent a magnetized compact stellar object surrounded by a charged fluid.

One approach in modern cosmology consists in supposing that dark matter and dark energy are different manifestations of a single `quartessential' fluid. Following such idea, this work presents a study of the evolution of perturbations of density in a flat cosmological model with a modified Chaplygin gas acting as a single component. Our goal is to obtain properties of the model which can be used to distinguish it from another cosmological models which have the same solutions for the general evolution of the scale factor of the universe, without the construction of the power spectrum. Our analytical results, which alone can be used to uniquely characterize the specific model studied in our work, show that the evolution of the density contrast can be seen, at least in one particular case, as composed by a spheroidal wave function. We also present a numerical analysis which clearly indicates as one interesting feature of the model the appearence of peaks in the evolution of the density constrast.

The stability of a recently proposed general relativistic model of galaxies is studied in some detail. This model is a general relativistic version of the well known Miyamoto-Nagai model that represents well a thick galactic disk. The stability of the disk is investigated under a general first order perturbation keeping the spacetime metric frozen (no gravitational radiation is taken into account). We find that the stability is associated with the thickness of the disk. We have that flat galaxies have more not-stable modes than the thick ones i.e., flat galaxies have a tendency to form more complex structures like rings, bars and spiral arms.

We study the stability of general relativistic static thick disks, as an application we consider the thick disk generated by applying the ``displace, cut, fill and reflect" method, usually known as the image method, to the Schwarzschild metric in isotropic coordinates. The isotropic Schwarzschild thick disk obtained from this method is the simplest model to describe, in the context of General Relativity, real thick galaxies. The stability under a general first order perturbation of the disk is investigated. The first order perturbation, when applying to the conservation equations, leads to a set of differential equations that are, in general, not self-consistent. This opens the possibility of performing various kinds of perturbations to transform the resulting system of equations into a self-consistent system. We perform a complete classification of the perturbations as well as the stability analysis for all the relevant physical perturbations. We found that, in general, the isotropic Schwarzschild thick disk is stable under these kinds of perturbations.

We describe an implicit procedure for solving linear equation systems resulting from the discretization of the three dimensional (seven variables) linear Fokker-Planck equation. The discretization of the Fokker-Planck equation is performed using a twenty-five point molecule that leads to a coefficient matrix with equal number of diagonals. The method is an extension of Stone's implicit procedure, includes a vast class of collision terms and can be applied to stationary or non stationary problems with different discretizations in time. Test calculations and comparisons with other methods are presented in two stationary examples, including an astrophysical application for the Miyamoto-Nagai disk potential for a typical galaxy.

We find a numerical self-consistent stellar model by finding the distribution function of a thin disk that satisfies simultaneously the Fokker-Planck and Poisson equations. The solution of the Fokker-Planck equation is found by a direct numerical solver using finite differences and a variation of Stone's method. The collision term in the Fokker-Planck equation is found using the local approximation and the Rosenbluth potentials. The resulting diffusion coefficients are explicitly evaluated using a Maxwellian distribution for the field stars. As a paradigmatic example, we apply the numerical formalism to find the distribution function of a Kuzmin-Toomre thin disk. This example is studied in some detail showing that the method applies to a large family of actual galaxies.

The stability of general relativistic thin disks is investigated under a general first order perturbation of the energy momentum tensor. In particular, we consider temporal, radial and azimuthal "test matter" perturbations of the quantities involved on the plane $z=0$. We study the thin disks generated by applying the "displace, cut and reflect" method, usually known as the image method, to the Schwarzschild metric in isotropic coordinates and to the Chazy-Curzon metric and the Zipoy-Voorhees metric ($\gamma$-metric) in Weyl coordinates. In the case of the isotropic Schwarzschild thin disk, where a radial pressure is present to support the gravitational attraction, the disk is stable and the perturbation favors the formation of rings. Also, we found the expected result that the thin disk models generated by the Chazy-Curzon and Zipoy-Voorhees metric with only azimuthal pressure are not stable under a general first order perturbation

A numerical solution of Einstein field equations for a spherical symmetric and stationary system of identical and auto-gravitating particles in phase transition is presented. The fluid possess a perfect fluid energy momentum tensor, and the internal interactions of the system are represented by a van der Walls like equation of state able to describe a first order phase transition of the type gas-liquid. We find that the space-time curvature, the radial component of the metric, and the pressure and density show discontinuities in their radial derivatives in the phase coexistence region. This region is found to be a spherical surface concentric with the star and the system can be thought as a foliation of acronal, concentric and isobaric surfaces in which the coexistence of phases occurs in only one of these surfaces. This kind of system can be used to represent a star with a high energy density core and low energy density mantle in hydrodynamic equilibrium.