In this work, we compute rates of merging neutron stars (MNS) in galaxies of different morphological type, as well as the cosmic MNS rate in a unitary volume of the Universe adopting different cosmological scenarios. Our aim is to provide predictions of kilonova rates for future observations both at low and high redshift. In the adopted galaxy models, we take into account the production of r-process elements either by MNS or core-collapse supernovae. In computing the MNS rates we adopt either a constant total time delay for merging (10 Myr) or a distribution function of such delays. Our main conclusions are: i) the observed present time MNS rate in our Galaxy is well reproduced either with a constant time delay or a distribution function $\propto t^{-1}$. The [Eu/Fe] vs. [Fe/H] relation in the Milky Way can be well reproduced with only MNS, if the time delay is short and constant. If the distribution function of delays is adopted, core-collapse supernovae as are also required. ii) The present time cosmic MNS rate can be well reproduced in any cosmological scenario, either pure luminosity evolution or a typical hierarchical one, and spirals are the main contributors to it. iii) The spirals are the major contributors to the cosmic MNS at all redshifts in hierarchical scenarios. In the pure luminosity evolution scenario, the spirals are the major contributors locally, whereas at high redshift ellipticals dominate. iv) The predicted cosmic MNS rate well agrees with the cosmic rate of short Gamma Ray Bursts if the distribution function of delays is adopted, in a cosmological hierarchical scenario observationally derived. v) Future observations of Kilonovae in ellipticals will allow to disentangle among constant or a distribution of time delays as well as among different cosmological scenarios.

Probing the origin of r-process elements in the universe represents a multi-disciplinary challenge. We review the observational evidence that probe the properties of r-process sites, and address them using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contributed to the total mass of r-process elements present in our Galaxy. We found discrepancies with the neutron star (NS-NS) merger scenario. Assuming they are the only site, the decreasing trend of [Eu/Fe] at [Fe/H]\,$>-1$ in the disk of the Milky Way cannot be reproduced while accounting for the delay-time distribution (DTD) of coalescence times ($\propto~t^{-1}$) derived from short gamma-ray bursts and population synthesis models. Steeper DTD functions ($\propto~t^{-1.5}$) or power laws combined with a strong burst of mergers before the onset of Type~Ia supernovae can reproduce the [Eu/Fe] trend, but this scenario is inconsistent with the similar fraction of short gamma-ray bursts and Type~Ia supernovae occurring in early-type galaxies, and reduces the probability of detecting GW170817 in an early-type galaxy. One solution is to assume an extra production site of Eu that would be active in the early universe, but would fade away with increasing metallicity. If this is correct, this extra site could be responsible for roughly 50% of the Eu production in the early universe, before the onset of Type~Ia supernovae. Rare classes of supernovae could be this additional r-process source, but hydrodynamic simulations still need to ensure the conditions for a robust r-process pattern.

The merging of two neutron stars (MNS) is thought to be the source of short gamma-ray bursts (SGRB) and gravitational wave transients, as well as the main production site of r-process elements like Eu. We have derived a new delay time distribution (DTD) for MNS systems from theoretical considerations and we have tested it against (i) the SGRB redshift distribution and (ii) the Galactic evolution of Eu and Fe, in particular the [Eu/Fe] vs [Fe/H] relation. For comparison, we also tested other DTDs, as proposed in the literature. To address the first item, we have convolved the DTD with the cosmic star formation rate, while for the second we have employed a detailed chemical evolution model of the Milky Way. We have also varied the DTD of Type Ia SNe (the main Fe producers), the contribution to Eu production from core-collapse SNe, as well as explored the effect of a dependence on the metallicity of the occurrence probability of MNS. Our main results can be summarized as follows: (i) the SGRB redshift distribution can be fitted using DTDs for MNS that produce average timescales of 300-500 Myr; (ii) if the MNS are the sole producers of the Galactic Eu and the occurrence probability of MNS is constant the Eu production timescale must be on the order of <30 Myr; (iii) allowing for the Eu production in core-collapse SNe, or adopting a metallicity-dependent occurrence probability, allow us to reproduce both observational constraints, but many uncertainties are still present in both assumptions.