Probing the origin of r-process elements in the universe represents a multidisciplinary challenge. We review the observational evidence that probes 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. When we assume that 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 (∝t ⁻¹) derived from short gamma-ray bursts (GRBs) and population synthesis models. Steeper DTD functions (∝t ^−1.5) or power laws combined with a strong burst of mergers before the onset of supernovae (SNe) Ia can reproduce the [Eu/Fe] trend, but this scenario is inconsistent with the similar fraction of short GRBs and SNe Ia occurring in early-type galaxies, and it reduces the probability of detecting GW170817 in an early-type galaxy. One solution is to assume an additional production site of Eu that would be active in the early universe, but would fade away with increasing metallicity. If this is correct, this additional site could be responsible for roughly 50% of the Eu production in the early universe before the onset of SNe Ia. 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.

We study the evolution of europium (Eu) and barium (Ba) abundances in Local Group dwarf spheroidal and ultrafaint dwarf galaxies by means of detailed chemical evolution models and compare our results with new sets of homogeneous abundances. The adopted models include gas infall and outflow and have been previously tested. We investigate several production scenarios for r-process elements: merging neutron stars and magnetorotational-driven supernovae. Production of Ba through the main s-process acting in low- and intermediate-mass stars is considered as well. We also test different sets of nucleosynthesis yields. For merging neutron stars we adopt either a constant and short delay time for merging or a delay time distribution function. Our simulations show that (i) if r-process elements are produced only by a quick source, it is possible to reproduce the [Eu/Fe] versus [Fe/H], but those models fail in reproducing the [Ba/Fe] versus [Fe/H]. (ii) If r-process elements are produced only with longer delays the opposite happens. (iii) If both a quick source and a delayed one are adopted, such as magnetorotational-driven supernovae and merging neutron stars with a delay time distribution, the [Eu/Fe] abundance pattern is successfully reproduced, but models still fail in reproducing the [Ba/Fe]. (iv) On the other hand, the characteristic abundances of Reticulum II can be reproduced only if both the Eu and the r-process fraction of Ba are produced on short and constant time delays during a single merging event. We discuss also other possible interpretations, including an inhomogeneous mixing of gas that might characterize this galaxy.

We compute rates of merging neutron stars (MNS) in different galaxies, as well as the cosmic MNS rate in 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, the production of r-process elements either by MNS or core-collapse supernovae is taken into account. To compute the MNS rates, we adopt either a constant total time delay for merging (10 Myr) or a distribution function of such delays. We conclude (i) the observed present time MNS rate in our Galaxy is well reproduced either with a constant time delay or a distribution function ∝t−1. (ii) The [Eu/Fe] versus [Fe/H] relation can be well reproduced with only MNS, if the time delay is short and constant. If a distribution function of delays is adopted, core-collapse supernovae are also required. (iii) The present time cosmic MNS rate can be well reproduced in several cosmological scenarios. (iv) Spiral galaxies 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. (v) 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. (vi) Future observations of kilonovae in ellipticals will allow us to disentangle among constant or a distribution of time delays and among different cosmological scenarios.

Rocky planets with temperate conditions provide the best chance for discovering habitable worlds and life outside the Solar system. In the last decades, new instrumental facilities and large observational campaigns have been driven by the quest for habitable worlds. Climate models aimed at studying the habitability of rocky planets are essential tools to pay off these technological and observational endeavours. In this context, we present eos-estm, a fast and flexible model aimed at exploring the impact on habitability of multiple climate factors, including those unconstrained by observations. eos-estm is built on estm, a seasonal-latitudinal energy balance model featuring an advanced treatment of the meridional and vertical transport. The novel features of eos-estm include: (1) parametrizations for simulating the climate impact of oceans, land, ice, and clouds as a function of temperature and stellar zenith distance and (2) a procedure (eos) for calculating the radiative transfer in atmospheres with terrestrial and non-terrestrial compositions illuminated by solar- and non-solar-type stars. By feeding eos-estm with Earth’s stellar, orbital, and planetary parameters, we derive a reference model that satisfies a large number of observational constraints of the Earth’s climate system. Validation tests of non-terrestrial conditions yield predictions that are in line with comparable results obtained with a hierarchy of climate models. The application of eos-estm to planetary atmospheres in maximum greenhouse conditions demonstrates the possibility of tracking the snowball transition at the outer edge of the HZ for a variety of planetary parameters, paving the road for multiparametric studies of the HZ.

The merging rate of double neutron stars (DNS) has a great impact on many astrophysical issues, including the interpretation of gravitational waves signals, of the short gamma-ray bursts (GRBs), and of the chemical properties of stars in galaxies. Such rate depends on the distribution of the delay times (DDT) of the merging events. In this paper, we derive a theoretical DDT of merging DNS following from the characteristics of the clock controlling their evolution. We show that the shape of the DDT is governed by a few key parameters, primarily the lower limit and the slope of the distribution of the separation of the DNS systems at birth. With a parametric approach, we investigate on the observational constraints on the DDT from the cosmic rate of short GRBs and the europium-to-iron ratio in Milky Way stars, taken as tracer of the products of the explosion. We find that the local rate of DNS merging requires that $\sim \! 1 {{\ \rm per\ cent}}$ of neutron stars progenitors live in binary systems which end their evolution as merging DNS within a Hubble time. The redshift distribution of short GRBs does not yet provide a strong constraint on the shape of the DDT, although the best-fitting models have a shallow DDT. The chemical pattern in Milky Way stars requires an additional source of europium besides the products from merging DNS, which weakens the related requirement on the DDT. At present both constraints can be matched with the same DDT for merging DNS.