We present a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive elements synthesised in the dynamical ejecta of neutron star-neutron star (NS-NS) mergers. We have studied the material ejected from four NS-NS merger systems based on hydrodynamical simulations which handle neutrino effects in an elaborate way by including neutrino equilibration with matter in optically thick regions and re-absorption in optically thin regions. We find that the neutron richness of the dynamical ejecta is significantly affected by the neutrinos emitted by the post-merger remnant, in particular when compared to a case neglecting all neutrino interactions. Our nucleosynthesis results show that a solar-like distribution of r-process elements with mass numbers $A \gtrsim 90$ is produced, including a significant enrichment in Sr and a reduced production of actinides compared to simulations without inclusion of the nucleonic weak processes. The composition of the dynamically ejected matter as well as the corresponding rate of radioactive decay heating are found to be rather independent of the system mass asymmetry and the adopted equation of state. This approximate degeneracy in abundance pattern and heating rates can be favourable for extracting the ejecta properties from kilonova observations, at least if the dynamical component dominates the overall ejecta. Part II of this work will study the light curve produced by the dynamical ejecta of our four NS merger models.

The majority of existing results for the kilonova (or macronova) emission from material ejected during a neutron-star (NS) merger is based on (quasi-) one-zone models or manually constructed toy-model ejecta configurations. In this study, we present a kilonova analysis of the material ejected during the first $\sim 10\,$ ms of a NS merger, called dynamical ejecta, using directly the outflow trajectories from general relativistic smoothed-particle hydrodynamics simulations, including a sophisticated neutrino treatment and the corresponding nucleosynthesis results, which have been presented in Part I of this study. We employ a multidimensional two-moment radiation transport scheme with approximate M1 closure to evolve the photon field and use a heuristic prescription for the opacities found by calibration with atomic-physics-based reference results. We find that the photosphere is generically ellipsoidal but augmented with small-scale structure and produces emission that is about 1.5–3 times stronger towards the pole than the equator. The kilonova typically peaks after $0.7\!-\!1.5\,$ d in the near-infrared frequency regime with luminosities between $3\!-\!7\times 10^{40}\,$ erg s−1 and at photospheric temperatures of $2.2\!-\!2.8\times 10^3\,$ K. A softer equation of state or higher binary-mass asymmetry leads to a longer and brighter signal. Significant variations of the light curve are also obtained for models with artificially modified electron fractions, emphasizing the importance of a reliable neutrino-transport modelling. None of the models investigated here, which only consider dynamical ejecta, produces a transient as bright as AT2017gfo. The near-infrared peak of our models is incompatible with the early blue component of AT2017gfo.

Theoretically predicted yields of elements created by the rapid neutron capture (r-) process carry potentially large uncertainties associated with incomplete knowledge of nuclear properties and approximative hydrodynamical modelling of the matter ejection processes. We present an in-depth study of the nuclear uncertainties by varying theoretical nuclear input models that describe the experimentally unknown neutron-rich nuclei. This includes two frameworks for calculating the radiative neutron capture rates and 14 different models for nuclear masses, β-decay rates and fission properties. Our r-process nuclear network calculations are based on detailed hydrodynamical simulations of dynamically ejected material from NS-NS or NS-BH binary mergers plus the secular ejecta from BH-torus systems. The impact of nuclear uncertainties on the r-process abundance distribution and the early radioactive heating rate is found to be modest (within a factor of ∼20 for individual A > 90 abundances and a factor of 2 for the heating rate). However, the impact on the late-time heating rate is more significant and depends strongly on the contribution from fission. We witness significantly higher sensitivity to the nuclear physics input if only a single trajectory is used compared to considering ensembles with a much larger number of trajectories (ranging between 150 and 300), and the quantitative effects of the nuclear uncertainties strongly depend on the adopted conditions for the individual trajectory. We use the predicted Th/U ratio to estimate the cosmochronometric age of six metal-poor stars and find the impact of the nuclear uncertainties to be up to 2 Gyr.

The nuclear level density and γ-decay strength of 192 Os have been extracted using particle-γ coincidence data from the 192 Os (α, α γ) 192 Os reaction by means of the Oslo method. The level density is found to be a rather smooth function of excitation energy, approximately following the constant temperature model. The γ-decay strength is compared to photoneutron cross-section data above the neutron separation energy, and to E1 and M1 strengths for nuclei in this mass region derived from primary transitions following neutron capture. Our results are in good agreement with these previous data and draw a consistent picture of the γ-strength function in the range E γ ≈ 1.5 − 6 MeV.Using the measured nuclear level density and γ-decay strength as input to the nuclear-reaction code TALYS, we provide the first experimentally constrained Maxwellian-averaged cross section (MACS) for the 191 Os(n, γ) 192 Os reaction relevant to s-process nucleosynthesis. The systematic uncertainties introduced by the normalization procedure of the level density and γ-strength function were investigated and propagated to the calculated Maxwellian-averaged cross section. The obtained result of the Maxwellian-averaged cross section at k B T = 30 keV, σ n,γ = 1134 ± 375 mb, is in very good agreement with the theoretical estimate provided by the KADo-NiS project, giving experimental support to the adopted KADoNiS value. Good agreement is also found with MACS values obtained from other libraries, such as TENDL-2017, ENDF/B-VII.0, and JEFF.