Constraining r-process nucleosynthesis with multi-objective Galactic chemical evolution models

The astrophysical site(s) of the r-process are uncertain, with candidates such as neutron star mergers and magneto-rotational supernovae predicting different event rates, delay times, and heavy-element yields. Galactic chemical evolution models constrain these properties by comparing model predictions with observed abundances. We explore, in a systematic and data-driven way, the astrophysical conditions under which r-process enrichment can reproduce the observed trends of multiple neutron-capture elements in the Milky Way. Rather than assuming a fixed site, we adopt a flexible, parametric approach to test whether a common set of r-process parameters can explain the chemical evolution of several heavy elements. We compute a grid of one-infall, homogeneous models varying: Eu yield per event, r-process event rate, enrichment delay time, and progenitor mass range. For each of the $\sim 1.5 \times 10^5$ models, we predict [X/Fe] vs. [Fe/H] trends by scaling Eu yields with the solar r-process pattern. A multi-objective optimisation based on Pareto fronts identifies models that best reproduce the abundance trends. Best-fitting models favour short delay times ($\leq 30\ \rm Myr$), low-mass progenitors ($\sim 20-25\ \rm M_\odot$), and an effective Eu injection of $\sim 2 \times 10^{-7}\ \rm M_\odot$ per event. Stars more massive than $\sim 80\ \rm M_\odot$ are too rare to dominate the enrichment. While heavy elements can be reproduced, lighter ones show stronger conflicts with Eu, reflecting that the solar r-process scaling relation becomes less valid toward lighter elements. No single class of r-process events, under solar-scaled yields, can explain light and heavy neutron-capture elements; at least two components are required: a main r-process consistent with solar and r-rich stars, and a weaker component producing enhanced light r-process elements, similar to that observed in r-poor stars.