1. Introduction
There are 17 short-lived radioactive (SLR, with half lives of the order of 0.1 to 100 Myr) nuclei known to have been (or potentially have been) present in the early Solar System (ESS) [
1]. Three of them,
I,
Pu, and
Cm, have the specific property to be produced in the Galaxy almost exclusively by the process of
neutron captures (the
r process). Among those three, live
Pu from the present-time interstellar medium has also been detected in young sediments of the ocean floor [
2,
3]. Furthermore,
Pu and
Cm are actinides located beyond Pb and Bi at mass numbers around 208–210, the end point of the
neutron-capture (
s) process [
4]. Therefore, they are exclusively of
r-process origin. Being located beyond the classical third
r-process peak at Pt and Au, actinides are typically produced if the number of neutrons per seed is relatively large. Instead,
I belongs to the classical second
r-process peak. It has only a very minor (a few percent) contribution from the
s process because the unstable isotope that precedes it on the
s-process path,
I, has a half life of 25 min only and decays faster than the typical time required for capturing a neutron. These three
r-process isotopes have most likely the same
r-process origin (as indicated by elemental abundance observed in halo stars, [
5]). They can be studied individually or together to provide evidence on the history of the material that made up the Solar System [
6] and to set constraints on the
r-process astrophysical site and its nuclear input, which are both extremely uncertain [
7,
8,
9,
10]. In particular, Côté et al. [
8] (hereafter Paper I) constrained the last
r-process source to have contributed to the solar material by comparing the
I/
Cm ratio observed in primitive meteorites to nucleosynthesis calculations based on neutron star (NS-NS) merger, black hole–neutron star (NS-BH) merger, and magneto-rotational supernova simulations. Here, we extend that study to
Pu, to investigate if it is possible to find an explanation for the presence of this SLR isotope in the ESS compatible with the explanation for the presence of the other SLR isotopes heavier than iron, also well known to have been present in the ESS.
Table 1 summarises the main properties and information available on the four isotopic ratios under consideration here:
I/
I,
Pu/
U,
Cm/
U, and
I/
Cm. We only analyse isotopic ratios because the most direct evidence that comes from the analysis of meteoritic material on ESS values is not absolute abundances, but abundance values relative to each other. Absolute abundances suffer from many uncertainties, e.g., chemical separation in the nebula, in the meteorite parent body, and/or during chemical analysis, as well as dilution from the original stellar source. The ratios of interest are those of each estimated SLR abundance relative to a long-lived or a stable isotope. These ratios are directly measured in primitive meteorites and their components (the first three rows of
Table 1), or derived from the ratios directly measured (last row), as in the case of the
I/
Cm ratio. This last ratio provides us with a further observational constraint because
I and
Cm have very similar half lives [
11]
1. This allowed to remove several theoretical uncertainties in Paper I, providing a direct window into the astrophysical conditions of the
r-process site that produced the
I and
Cm in the ESS. Note that, instead, it is not possible to extract any further meaningful constraints from the
I/
Pu and
Cm/
Pu ratios because their half lives are very different from each other [
11].
Out of the four ratios reported in
Table 1,
Pu/
U has not been considered yet within a global analysis of origin of the SLR nuclei heavier than iron in the ESS. This is for two main reasons: first, its half life of 80 Myr is very different from that of the other two isotopes of roughly 15 Myr, therefore, the modelling of its abundance in the interstellar medium (ISM) is likely to present a different behaviour (see discussion in
Section 3). Second, its ESS abundance is less certain than those of the other two isotopes. The ESS
I/
I ratio has an uncertainty of roughly 2% at 2
and many studies agree on its value, suggesting that systematic uncertainties are not significant [
12]. The
Cm/
U was established with an uncertainty of roughly 6% at 2
thanks to the discovery of the special meteoritic inclusion, named Curious Marie, rich in U [
13]. More data on different samples is still needed to completely establish this value.
In the case of the ESS
Pu abundance (i.e., the
Pu/
U ratio), instead, not only is the uncertainty for the value reported in
Table 1 roughly 30%, but also there are potential systematic uncertainties in the determination of the ESS value. The ESS
Pu abundance can be estimated by xenon isotope studies of meteorites, since
Xe and the heavy
Xe are stable isotopes produced by the spontaneous fission of
Pu. Moreover, solids are extremely poor in noble gases and the radiogenic and fissiogenic xenon signatures may become significant over time and, hence, can be quantified at high precision. Studies have been focusing on gas-poor meteoritic materials: mineral separates [
14]; CAIs [
15,
16]; differentiated meteorites with simple cooling histories, such as angrites [
17]) and eucrites [
18]; and high-metamorphic-grade ordinary chondrites [
19]. Currently, there are two “best estimates” of ESS using different approaches. Lugmair and Marti [
17] normalized
Pu to
Nd, an
r-process-only isotope of Nd, because they found an achondrite (Angra dos Reis) where they could prove that the geochemical analogue of Pu is Nd, and potential modification of the Pu/Nd ratio with respect to the Solar System abundances can be ruled out. They reported
Pu/
U ratios ≃ 0.0043 at the adjusted time of Solar-System formation [
20]. The value reported in
Table 1 (0.007) is a different estimate by [
19], who used a different approach. As the fissiogenic signature is dominated by
Pu-derived xenon in meteorites, they irradiated an exceptionally gas-poor ordinary chondrite (St Severin) with thermal neutrons to induce the fission of
U and derived the
Pu/
U ratio from the component analysis of xenon isotope measurements alone. This value is almost twice as high as the value provided by the Angra dos Reis study, and it is in better agreement with the more recent analysis of Xe in ancient terrestrial zircons from Western Australia [
21]. In summary, the major challenge is to find a meteorite sample that is representative of the Solar System, and for which geochemical processes that could potentially modify the relative abundances of Pu to U or rare earth elements with respect to the chondritic composition is well-understood, and the effect can be corrected for. Here, we will consider for the ESS
Pu/
U ratio the value reported in
Table 1. If the “true” value was eventually found to be lower, for example, by a factor of two, all the times calculated and reported in our analysis below would have to be increased by 80 Myr.
The aim of this paper is to investigate possible self-consistent solutions for the origin of the abundances of all the SLR nuclei heavier than iron observed to have been present in the ESS, including
Pu. These observed abundances are represented by the four
r-process ratios reported in
Table 1, as well as the SLR isotopes produced by
neutron captures (the
s process, specifically
Pd,
Cs, and
Hf, as discussed in [
22]). We start by discussing predictions from state-of-the-art models of the
r process for the three SLR nuclei of interest and their reference isotopes (
Section 2). Then, in
Section 3, we consider the temporal evolution of the
Pu/
U ratio in the ISM and discuss if there are solutions for its ESS value that are consistent with the abundances of the other SLR nuclei heavier than iron. Finally, in
Section 4, we present our summary, conclusions, and suggestions for future work.
2. Nucleosynthesis Calculations
We consider the large set of
r-process abundances published with Paper I and calculated with the nuclear network code WINNET
2 [
23,
24] and the nucleosynthesis network PRISM
3 [
25,
26]. All the abundances reported and used in this work, are taken at 1 Myr after the nucleosynthetic event, i.e., they are not decayed completely, given that we are interested in SLR nuclei.
Table 2 lists all the WINNET models considered here and the relationship between the labels used in Paper I and the shorter labels used here. The sites and the nuclear physics sets, with all their relevant references, are described in detail in Paper I and Ref. [
27]. Here, we remind briefly that the nomenclature of the nuclear input is as follows: [D,J,Jm] denotes the mass model (D for Duflo Zuker, J for JINA reaclib, Jm for JINA with Marketin theoretical
decays). The “h” indicates that the nuclear heating subroutine was turned on [
28]), modifying the temperature evolution of the trajectory. Finally, [f1,f2,f4] represent three different fission fragment models. There are, in total, 3 (top) × 3 (bottom) nuclear labels, i.e., nine sets of nuclear inputs (right side of the table), and seven astrophysical sites (left side of the table), therefore, a total of 63 WINNET models. The tabulated abundances of the six isotopes of interest here (together with the Eu isotopes and
Th for future reference) can be found in
Supplementary Table S1.
In the case of the PRISM models, the nomenclature is identical to that used in Paper I (
Table S3), and references therein. In this case, four sites are considered (the dynamical ejecta of two NS–NS mergers and two NS–BH models), and combined with ten different mass models, of which four are also investigated using alternative
decays (the “Mkt” label, here, corresponding to the “D3C*” label in Paper I). The total is, therefore, 4 × 14 = 56 PRISM models. The tabulated ratios of interest here can be found in
Supplementary Table S2.
The four ratios of interest from all the models are plotted in
Figure 1 and
Figure 2. As expected, ratios of isotopes of similar mass (
Figure 2) are much less dependent on the model than the
I/
Cm ratio (
Figure 1). Variations in those three ratios are typically of factors ∼2.5 to 3 in the WINNET models. Additionally found by some of these models is the
I/
I ratio of 1.35 derived from the
r-process abundance of the stable
Xe in the Solar System
4. When considering the PRISM models, which explored a larger set of nuclear inputs, variations are somewhat larger, especially in the case of
Cm/
U (up to a factor of 10).
All the models show variations in the
I/
Cm ratio of up to three orders of magnitudes (
Figure 1, corresponding to Figure S2 of Paper I, but with the PRISM models also included). Of the WINNET models, 13 of them (the nine FMdef models, the three FMs6 Jmhf models, and the FMs6 Jhf4 model) could match the observed
I/
Cm in the ESS. The PRISM models all represent dynamical ejecta and, therefore, provide similar results to the corresponding WINNET models (Bs125, R1010, R1450). As noted above, the PRISM models explore a larger set of nuclear inputs and only one out of those 14 choices (TF_Mkt) provides a solution for the
I/
Cm ratio in all the four sites.
A quick estimate indicates that self-consistent values of the time elapsed from the last
r-process event (
Section 3.1) would result from models with similar values of the
I/
I and
Cm/
U ratios. This is because both the ESS value and the
K values for
I/
I are roughly twice those of
Cm/
U. Therefore, the two differences cancel each other out in the calculation of the ISM ratio needed to derive the time interval of the decay by comparison to the ESS ratio
5. While there are no models with the same values of the
I/
I and
Cm/
U ratio, when we consider the uncertainties of the
and ESS values, many solutions can be found for a much larger range of relative ratios, as shown in
Section 3.1. This is because the time interval of the decay is a function of the natural logarithm of the abundance ratio; therefore, variations in the relative ratios up to a factor of 5 (or even 10) result in a difference, by subtraction, between the time intervals of 1.6 (2.3), which correspond to a percent variation by 30% (50%) only, i.e., well within the uncertainties. For
Pu/
U, instead, it is more difficult to make a quick estimate because of the very different
. In the next section, we evaluate quantitatively, using the WINNET set, the values of the time elapsed between production and incorporation into the first solids in the ESS to verify which models can match the three constraints simultaneously.
3. Galactic Evolution and Origin of the SLRs in the ESS
When considering the ESS data, we need to process the stellar abundances for their recycling within the ISM material from which the Sun formed. Such recycling implies a certain time delay, which is crucial to consider when analysing radioactive isotopes that decay within a given timescale. Côté et al. [
29,
30] and Yagüe López et al. [
11] provided a methodology and tools to address the evolution in SLR nuclei in the ISM of the Galaxy, and we base our analysis on such works.
First, we need to take into account the uncertainties related to galactic chemical evolution (GCE) itself over the whole lifetime of the Galaxy. These result in a factor
K, by which any ratio predicted by nucleosynthesis calculations involving a stable or long-lived reference isotope needs to be multiplied. This factor takes into account the history of the Galaxy and how it influences the evolution, and therefore the abundance, at the galactic time of the formation of the Sun, of the stable or long-lived isotope that is used as reference for the ESS ratio. The values of
K we calculated from the full GCE models [
29] are reported in
Table 1 for the three isotopic ratios considered here. Three values are provided: the middle value is the best-fit case and the other two reflect the GCE uncertainties, which provide a minimum and maximum value of the SLR to stable or long-lived isotope ratios. Summarizing
Table 1 of [
29], the GCE parameters that mostly affect the value of
K are those related to the first and second infall episodes (
and
) and the star formation efficiency (
). The observational constraints whose uncertainties affect
K the most are the current inflow rate and mass of gas. Due to the feedback between all these quantities, there is not a simple relation with the value of
K. For example, the
values are found for the highest values of
,
, and inflow rate, together with the lowest values of
and mass of gas. The reasons for this behaviour are explained in detail in [
29].
We found that, if the reference isotope is stable, as in the case of I, the best-fit value of K is 2.3. When the reference isotope is unstable and long-lived (such as U), the value of K decreases with the half life of the nucleus because the abundance is affected by a shorter time scale within the full history of the Galaxy. For example, in the case of U, with a half life of 1 Gyr, i.e., roughly ten times shorter than the age of the Galaxy, the K factor decreases by roughly a factor of two. In the case of I/Cm, there are no values of K to be applied (in other words, K is always equal to 1) because these two isotopes are both short-lived and insensitive to the past history of our Galaxy. This is one the several advantages of using such a ratio, as discussed in detail in Paper I.
The other potential problem is that injection of SLR nuclei into the ISM by the stellar objects that produce them is not continuous, because stellar ejection events happen in correspondence to very specific discrete events, e.g., supernova explosions or neutron–star mergers. For stable nuclei, this effect is not significant because their ESS abundances are primarily defined by the total number of events that enriched the pre-solar nebula, rather than by the exact times at which the events occurred. However, for SLR nuclei, this effect can completely control their abundances in the ISM since they freely decay between events. One way to account for this is to consider the average of the interval
between additions to a given parcel of ISM gas from events of a given type, and compare it to the mean-life
of the SLR nuclei produced by this type of events. Therefore, the
ratio is the crucial parameter to consider, or equivalently
, where
is the time interval between the births of the event progenitors
6. We do not know a priori the value of
for any SLRs and their sources because it depends on uncertain effects such as diffusive transport in the ISM, supernova energetic in carrying material in the ISM, spatial distribution of the events, and distance of the events from the pre-solar ISM parcel of gas (see, e.g., [
31] and Wehmeyer et al. in prep). Our approach has therefore been to first develop a general framework and then test its implications and derive its predictions for different values of
.
Côté et al. [
30] found that, if
where
is the time interval between the births of the event progenitors, then we can treat the injection of SLR from such an event as continuous (hereafter Regime I). We just need to add an uncertainty resulting from the statistical spread of the SLR abundance. If, instead,
, the most likely scenario is that the ESS abundance of the given SLR came from one last event only, without any memory of the previous events (hereafter Regime III). For values of
in-between 0.3 and 2, the SLR abundance carries the memory of a few events (hereafter Regime II). Finally, we note that considering SLR ratios such as the
I/
Cm ratio in the last row of
Table 1 significantly reduces the uncertainties resulting from the discrete nature of stellar ejections, especially for cases when the half lives are comparable, as discussed in general in Ref. [
11], and in detail for the
r-process SLR in Paper I.
In the following, we use the WINNET abundances to derive more information on the early Solar System and the history of presolar matter from the three r-process isotopes considered here in different possible scenarios. We remind that the main difference between I and Cm, on the one hand, and Pu, on the other hand, is that the half life of the latter is roughly five times longer than those of the former two. Therefore, the criterion needs to be applied differently, even if all the three isotopes are exclusively produced by r-process events.
3.1. One (Regime III) or Few (Regime II) Events and Time Elapsed from Last Event
In the case of the two
r-process SLR
I and
Cm, as discussed and presented in detail in Paper I, we can justify statistically the assumption that their abundances in the ESS originated from one last event only, (Regime III ), which occurred roughly 100–200 Myr before the formation of the first solids in the Solar System. The criterion
under which Regime III is valid also for
Pu is that
, or equivalently
in the equation, is greater than 345 Myr. Therefore, possible solutions for this scenario are those for which
is around or larger than this value. To evaluate the ISM
Pu/
U ratio under the assumption that
I,
Cm, and
Pu in the ESS originated from one event, we then use Eq. S2 of Paper I, as performed in that paper for
I/
I and
Cm/
U, and the values of
K reported in
Table 1. Some examples of the calculation of the time from the last event are shown in
Figure 3. There, self-consistent solutions are represented by the overlapping areas of the three different colored bands, each representing one of the three SLR isotopes and their uncertainties. The trend with the
of the time elapsed calculated using
Pu is steeper than those calculated using the other two isotopes. This is due to its much longer
value and the fact that the time elapsed is a linear function of
.
Figure 3 also shows some examples of possible solutions for Regime II, which corresponds to
, i.e.,
= 68
7–345 Myr. In this case,
Pu originated from a few discrete events and the lower the value of
, the larger the number of events. The last event would have contributed only a fraction,
(assuming a constant production factor), of the ESS abundance of
Pu. Therefore, at the lower limit of Regime II,
Myr, the last event contributed 45% of the ESS abundance of
Pu.
Overall, the WINNET set comprises 63 sets of models, and the GCE model provides three values of
K for a total of 189 possibilities. For the first three ratios of
Table 1, we found that 92% of the models can provide overlapping solutions: 62 of those with
, 60 of those with
, and 52 of those with
. Therefore, solutions are common, partly thanks to the degree of freedom provided by the relatively free parameter
. Times of the last event are in the range 100–200 Myr as derived in Paper I, this is expected given that the analysis presented here is just an extension of that presented there, to check if the
Pu/
U ratio could also be explained. A new result is that these elapsed times are lower in Regime II than in Regime III, due to the steeper trend with
of those calculated using
Pu.
If we include the requirement that the
I/
Cm ratio should be between 254 and 622
8, the number of solutions becomes much more restricted. In fact, the
I/
Cm ratio is a much more stringent constraint because the two isotopes are very far apart in mass and therefore located in very different regions of the nuclide chart. Such relative abundances are more sensitive to the general features of the process and its astrophysical site (such as the amount of free neutrons) as well as the uncertainties in the nuclear model, than the ratios of isotopes that are closer to each other in mass (see
Figure 2). As shown in Figure S2 of Paper I, the simulations that produce the best matches to the observations are those dominated by moderately neutron-rich ejecta (in the specific case of the WINNET models, these correspond to the nine FMdef and the three FMs6 Jmhf models).
Out of these models, we find that six out of the nine FMdef models (the three Dhf plus the three Jhf cases) and one of the FMs6 model (Jmhf4) can also account for Pu/U when using either of the three values of K, for 21 solutions in total. The main difference between using and versus using is that the former two values provide solutions for values typical of Regime III, while the latter corresponds to solutions within Regime II. (The FMs6 Jmhf1 and Jmhf2 cases produce I/Cm ratios of 236 and 242, respectively, just outside the required range). In summary, more than half (21) of all the 36 possible models that match the I/Cm (12 models × 3 values of K = 36) provide a global solution for all the four isotopic ratios.
Finally, we note that, if the ESS
Pu/
U ratio was lower than the value used here, the green shaded area in
Figure 3 would shift upwards, for example, by 80 Myr if the ESS ratio was twice as low, due to a longer decay time needed to match the lower ESS value. This would remove most of the Regime III solutions and shift the Regime II solutions to lower values of
.
3.2. Steady-State Equilibrium (Regime I) and Isolation Time
In Regime I,
is
, corresponding to
Myr, and
Pu evolves in steady-state equilibrium. In this case,
I and
Cm would be in Regime II; the time elapsed from the last event would decrease with
(as discussed in
Section 3.1), and reach roughly 80–130 Myr. The steady-state regime for
I and
Cm would require, instead, roughly
Myr (for this value of
, at the limit of Regime II, the last event contributed roughly 40% of their ESS abundances). This can be excluded with reasonable confidence because it is the typical value obtained for core-collapse supernovae, which are much more frequent than the currently accepted
r-process sources.
If
Pu was in steady-state equilibrium, then we can use Equation (11) of [
1] (where
) and consider the productions ratios from the
r-process models as a continuous wave of enrichment. In this case, the time interval needed to decay the ISM ratio to its corresponding ESS ratio is an isolation time rather than a time from the last event. This time interval can then be compared to the isolation time obtained from the
s-process isotopes,
Pd,
Cs, and
Hf, under the assumption of same regime, i.e.,
Myr for the
s-process events in the Galaxy, which correspond to asymptotic giant branch (AGB) stars of initial mass
M
[
22]. For the three values of
K to be used when studying SLR/stable isotope ratios (i.e., 1.6, 2.3, 5.7, as reported in
Table 1 for
I/
I), the isolation times reported by [
22] for the
s-process SLR isotopes are 9–12, 10–16, and 18–26 Myr, respectively. We also need to consider the statistical uncertainty due to stochasticity and discussed in [
30]. We can use here the uncertainties reported in Table 3 of [
30], for the specific case
and
= 31.6 Myr, which are close to the maximum uncertainty that would correspond to the case of
Pu in this regime. The error is almost symmetric and corresponds to variations in the ISM ratio of +1.16 and −0.84. These translates into error bars to be applied to each isolation time of +17 and −19 Myr.
In the case of
, no solutions are present because all the isolation times derived from
Pu are in the range 61–213 Myr, much higher than the range derived for the
s-process SLR nuclei of 18–26 Myr. This is controlled by the large value of
combined with the long half life of
Pu. In the case of
and
, instead, there are 18 and 14 solutions possible, respectively, which have an overlap with the ranges of isolation time derived from the
s-process SLR nuclei. These solutions are all obtained from the models that produce
Pu/
U abundance ratios in the range 0.19–0.34. As shown in
Figure 2, these correspond mostly to the WINNET models run with the Jmhf nuclear inputs (out of the total 32 solutions, 23, , i.e., 72%, are Jmhf solutions) and the six NS–NS merger PRISM models with SLY4, TF_Mkt, and UNEDF0. For the other nuclear input choices, instead, only specific astrophysical sites results in
Pu/
U abundance ratios in the required range.
Out of the 12 models that match the three ratios that involve
I and
Cm, seven of them also provide solutions for the isolation time from
Pu/
U compatible with the
s-process SLR isotopes. However, as mentioned above for the value of
considered here, these two SLR isotopes may have more than one event contributing to their ESS abundances; therefore, such constraints become less strong (see also [
31]).
We should also consider the case where
Pu is in steady-state, but the
s-process SLRs came from one last event, which requires roughly
Myr for
s-process event in the Galaxy. In this case, the last
s-process event was identified to have occurred at 25 Myr before the formation of the first solids [
22], therefore, the isolation time from
Pu/
U is simply constrained to be smaller than this value. Also in this case solutions do not exist with
, while there are 13 and 7 more solutions for
and
, respectively. Most of these solutions overlap as they correspond to the same
r-process models but the different value of
K, therefore, they correspond to the same range of
Pu/
U abundance ratios and nuclear models as reported above. The isolation time from
Pu in this case can vary more more freely and there are a few models that given the uncertainties provide values down to zero, which is not a useful constraint. Finally, we note that if the
s-process SLRs originated from a few events (i.e.,
Myr), then the time from the last event would increase and a few more models could produce an isolation time lower than this value. A more detailed statistical analysis would be needed in this case.
Finally, we note that if the ESS value of the Pu/U ratio was twice as high as the value considered here, we would need to add Myr to every isolation time, which would make it impossible to find a solution consistent with the origin of the s-process SLR isotopes.
4. Summary and Conclusions
We presented and analysed the relative production of the short-lived and long-lived
r-process isotopes
I,
U,
U,
Pu, and
Cm and the stable
I in a large set of 119
r-process models from two different sets calculated with the WINNET and PRISM frameworks. We then investigated if it is possible to find solutions for the origin of the ESS abundance of
Pu that provide production at the source and time intervals (either from the last event or from the time of the isolation, depending on the
regime) compatible to those of the other
r-process and
s-process SLR isotopes. A summary of the different possibilities, solutions, and derived time intervals are shown in
Table 3. In brief:
In
Section 3.1 (top section of
Table 3), we considered Regimes II and III for
Pu, corresponding to
Myr and Regime III for
I and
Cm. More than half of the WINNET models that were already shown to reproduce the three ratios that involve
I and
Cm in Paper I, also provide a self-consistent solution for
Pu. These models all correspond to the NS–NS merger disk cases dominated by moderately neutron-rich ejecta.
In
Section 3.2 (bottom section of
Table 3), we considered Regime I for
Pu, i.e.,
Myr, where this SLR reaches a steady-state value in the ISM. It is also possible to find a significant number of
r-process models (mostly corresponding to the Jmhf nuclear input) that provide solutions for the ESS
Pu abundance compatible to those of the SLR isotopes produced also by the
s process:
Pd and
Hf (and the current ESS upper limit of
Cs). However, no solutions exist in Regime I for
Pu if the ESS value of
Pu was twice as high as the value used here or if the Milky Way model was represented by
.
We cannot determine if the solution to the origin of
Pu in the ESS is 1. or 2. above, and which implications on the timescales are valid, since we still do not know how far material from
r-process sources can travel, and, therefore, how many parcels of the ISM are affected by each of these events and the value of
. Although we note that recent hydrodynamical models aimed at calculating how far material travels after being ejected by a hypernova predict relatively short distances [
32], which would support large
values for the r-process events. In any case, we have established that within Point 1., WINNET solutions within the NS–NS disk ejecta favour the Dhf and Jhf nuclear models. Within Point 2., all solutions exclude the case of a Milky Way Galaxy with
, therefore restricting the isolation time to 9-16 Myr (if
Pd and
Hf are also in Regime I), still supporting the hypothesis that the Sun was born in a massive, long-living molecular cloud. We can also conclude that a much lower value of the
Pu/
U ratio in the ESS than that reported in
Table 1 would be impossible to reconcile within Regime I, and would therefore support Regimes II and III. New, future experiments and analysis are needed to confirm the ESS
Pu/
U ratio.