From the slow to the rapid neutron capture process

This contribution starts with memories about Franz Käppeler, as a human being and a scientist, his impact on neutron capture nucleosynthesis via experiments and astrophysical studies, before extending from his focus (the s-process) to the r-process, discussing its mechanism, related astrophysical sites and examining in a concluding section its role during galactic evolution.


Introduction
Before starting the science related to this article, let me say a few words about Franz. He has been a highlight in terms of science, excellence, creative thinking, diligence, and character. We met first at one of the Ringberg Workshops in Nuclear Astrophysics, where he participated together with his colleagues Hermann Beer and Klaus Wisshack. One of the important decisions regarding his science was that he turned from measurements of nuclear cross sections for reactor applications to neutron captures in astrophysics. The idea (in 1988) to utilize the 7 Li(p,n)-reaction to provide a neutron spectrum which corresponds to a Maxwellian neutron energy distribution, equivalent to a thermal environment of kT = 25.3 keV [1], was revolutionary and permitted to measure directly already energy/velocity-averaged neutron capture rates for typical s-process conditions. This and his enormous productivity, filling a number of Atomic and Nuclear Data volumes [2][3][4], which were then readily applied to s-process studies [5][6][7], led to a tremendous impact in that fields, spreading from Karlsruhe around the world, including finally also nToF at CERN [8,9]. Franz was a centerpiece of the field, he brought the second Nuclei in the Cosmos (II) conference to Karlruhe, educated a large family of students who are now at important places world-wide, and collaborated energetically with many colleagues, covering experiments as well as theory. I am proud of 14 joint publications, a nice Nuclear and a e-mail: f-k.thielemann@unibas.ch (corresponding author) Particle Science review with him and Michael Wiescher [10], and we had a joint graduate student [Iris Dillmann, [11,12]]. Shortly after I returned from the US to Europe in 1994, I took the freeway from Frankfurt/Darmstadt to Basel (my then new home) and decided to pick an exit to visit Franz' home, which I had known from previous visits. It was a surprise for all of us, finding a joyful group of visitors, Michael Wiescher and Roberto Gallino were visiting for joint work [13]. Figure 1 (right panel) is testimony of this highly enjoyable encounter.

What features determine the outcome of neutron capture processes?
Franz has mostly worked on the s-process, i.e. slow neutron captures, where -if the resulting isotope is unstable -the betadecay wins in the majority of cases, and one can just think of a sequence of neutron captures connecting nuclei A and A + 1. The identification of a participating isotope is then unique for each A, as one can essentially think of the beta-decay of an instantaneous transition in comparison to neutron captures. Therefore, the speed of the build-up of heavy elements is determined by neutron capture cross sections / rates with the slowest reactions for small cross sections related to the smallest reaction Q-values at magic neutron numbers. If neutron densities n n are enhanced, the rates increase proportional to the neutron density and neutron capture life times (proportional to 1/n n ) can become shorter than the betadecay of unstable nuclei. We display this behavior in Fig. 2 for an abundance composition starting with initial solar abundance for 127 I, 235 U and 238 U in an environment with a neutron density of n n = 6 × 10 18 cm −3 . One can see that within less than 10 −4 s neutron captures hit the limiting isotopes of the nuclear network utilized here.
This leads to the two major neutron capture processes, the s-and the r-process which are both indicated in Fig. 3. (A modification of the s-process at slightly higher neutron den-  [13]]. BTW, Franz had smuggled Rosel in a Fiat 500, hidden below the back seat, from East to West Germany long before the fall of the Berlin wall. He was not only creative, but had an enormous courage Fig. 2 Starting from solar abundances of 127 I, 235 U and 238 U, in a plasma with a neutronen density of n n = 6 × 10 18 cm −3 and a temperature of 10 9 K, it takes typically less than 10 −4 s until neutron captures hit the limiting isotopes included in the nucleosynthesis network sities with n n up to 10 16 cm −3 , the i-process, is discussed in the chapter by Falk Herwig. In that case the unique identification of a mass number A with a specific isotope is lost, as branchings with comparable beta-decay lifetimes come into place). The rapid neutron capture process is not only dominated by fast neutron captures, in most conditions the inverse photo-disintegrations play an essential role as well. For a reaction i(n, γ )m, with i standing for nucleus (Z , A) and m for (Z , A + 1), the relation between the forward rate σ v i;n,γ and the photodisintegration rate λ m;γ,n is given by containing the reduced mass μ in and the neutron capture Qvalue for the reaction, g n = 2×(1/2)+1 for the neutrons and the partition functions G for the nuclei involved. In chemical equilibrium the difference between forward neutron capture and backward photodisintegration flux for nucleus i vanishes, i.e.Ẏ i = −n n σ v i;n,γ Y i + λ m;γ,n Y m = 0. Making also use of the relation between σ v i;n,γ and λ m;γ,n , based on detailed balance, results in For the two adjacent isotopes (Z , A) and (Z , A+1) or (Z , N ) and (Z , N +1), connected by neutron capture, the Q-value is equivalent to the neutron separation energy S n of (Z , A + 1). Equation (2) leads to an abundance distribution in an isotopic chain for element Z with a maximum at a unique neutron separation energy, identical for all different elements Z . This leads to a process path along a contour line for a specific neutron separation energy S n which depends on the neutron density n n and the temperature T . A specific choice with a typical S n ≈ 2MeV is shown in Fig. 3, which indicates the main difference between s-and r-process features. The s-process (with neutron number densities n n of 10 6 to 10 11 cm −3 ) is strongly affected by neutron captures with the smallest rates, leading to the slowest process speed and abundance pile-ups at neutron magic numbers for nuclei at or close to stability. In the r-process (with neutron densities n n > 10 20 cm −3 ) the neutron captures are extremely fast (see Fig. 2) which is also the case for reverse photodisintegrations of very neutron-rich nuclei with temperatures in excess of 10 9 K. Therefore, a chemical or (n, γ ) − (γ , n) equilibrium is obtained (almost) instantaneously for each isotopic chain Z . The process speed is determined by the beta-decay rates between isotopic chains for nuclei at the position of the rprocess path, i.e. those at contour lines of constant neutron separation energies. These contour lines experience kinks at closed neutron shells (far from stability) with longest halflives for the points closest to stability, i.e. at the top of these Fig. 3 (courtesy EMMI, GSI/Different Arts) Nuclear chart with major nucleosynthesis processes, (i) the formation of elements up to Fe, Ni, Zn in fusion reactions during stellar evolution and explosions, (ii) the s-process path along stability up to Pb and Bi, and (iii) the (later to be discussed) r-process path. Stable nuclei are indicated by black squares, experimentally known β + -unstable nuclei by red squares, β − -unstable nuclei by blue squares, alpha-unstable nuclei by yellow squares. The area of β-unstable nuclei which are still stable against neutron or proton emission, as predicted by nuclear mass models, appear in a light green shade. Horizontal or vertical lines stand for neutron or proton shell closures, N = 8, 20, 28, 50, 82, 126, 184 and Z = 8, 20, 26, 28, 50, 82 (, 114), the one in brackets not shown here  [14]: solar s-process abundances (blue), resulting from a superposition of stellar s-process sources [6]. The subtraction from global solar abundances provides the abundance contributions due to the remaining processes which create heavy nuclei up to the actinides, mainly the before-mentioned and later discussed r-process (red) kinks. Therefore, the abundance maxima of the r-process are shifted to the left of the s-process maxima, because the betadecay of isotopes at the top of these kinks leads to isotopes at stability with a smaller A. This is shown in Fig. 4.
We know from the ample amount of Franz Käppeler's work, that the solar s-process contribution is caused by a superposition of several processes, i.e. the weak and strong s-process from massive and low to intermediate mass stars (with variations according to the contributions from the stellar mass interval as well as the stellar metallicity which determines the neutron to Fe-seed ratio of this secondary nucleosynthesis process). Opposite to this it has long been believed that the solar r-process abundances are due to unique process features. The following sections will discuss this in more detail, combined with an attempt to present an overview of the present understanding of the possible astrophysical sites, with the result that also the solar r-process abundances are caused by a superposition of different contributions.

Initial complete chemical equilibrium, NSE
Opposite to the s-process, which works on typical neutroncapture timescales of months to years with neutrons being produced via (α, n)-reactions during stellar evolution, the high density of neutrons (with a half-life of the order 10 min), required for an r-process to occur, indicate an explosive environment. This starts typically with a subset of explosive Si-burning, where temperatures in excess of 5-6×10 9 K are attained, leading initially to a nuclear statistical equilibrium, i.e. a complete set of chemical equilibria, balancing all forward and backward flows in reactions, in particular also for proton or neutron capture reactions p + (Z , A) (Z + 1, A + 1) + γ and n + (Z , A) (Z , A + 1) + γ , corresponding to a relation between the chemical potentials μ p + μ(Z , A) = μ(Z + 1, A + 1) and μ n + μ(Z , A) = μ(Z , A + 1), because the chemical potential of photons vanishes. Such relations between the chemical potentials lead to identical equations for abundance ratios as given for neutron captures in Eq.(2), derived there via expressions for the reaction rates. If this is not only the case for a particular reaction, but across the whole nuclear chart, the complete reaction sequence is in chemical equilibrium, i.e. Z μ p + N μ n = μ(Z , A), termed complete chemical or also nuclear statistical equilibrium NSE [e.g. [15,16]]. For Boltzmann distributions (which apply in general in astrophysical plasmas, with the exception of highly degenerate conditions, where Fermi distributions have to be utilized for the chemical potentials), the abundances of nuclei can be expressed by nuclear properties like the binding energies B i (Z i , A i ), the abundances of free neutrons and protons Y n and Y p , and environment conditions like temperatures T and densities ρ, leading to the abundance of nucleus i [with Z i protons and N i neutrons or A i = Z i + N i nucleons, see e.g. [15]] where G i corresponds to the partition function of nucleus i, as the ground and excited state population is in thermal equilibrium. Reactions moderated by the weak interaction, i.e. β-decays, electron captures, and charged-current neutrino interactions, change the overall proton to nucleon ratio also called electron fraction) and occur on longer time scales than particle captures and photodisintegrations. They are not necessarily in equilibrium and have to be followed explicitly. Thus, as a function of time the NSE will follow the corresponding densities ρ(t), temperatures T (t), and Y e (t), determining the two unknowns Y n and Y p as a solution of the equations for total mass conservation and Y e resulting from weak interactions This permits then the prediction of the individual nuclear abundances Y i from Eq.(3). In general, very high densities favor large nuclei, due to the high power of ρ A i −1 , and very high temperatures favor light nuclei, due to (kT ) −3(A i −1)/2 in Eq. (3). In the intermediate regime exp(B i /kT ) favors tightly bound nuclei with the highest binding energies in the mass range A =50-60 of the Fe-group, but depending on the given Y e . The width of the composition distribution is determined by the temperature.
Under certain conditions, i.e. not sufficiently high temperatures, when not all reactions are fast enough (especially due to small reaction rates caused by too small Q-values, e.g. proton or neutron binding energies across magic proton or neutron numbers, i.e. closed shells), not a full NSE emerges, but only certain areas of the nuclear chart are in equilibrium. These are called quasi-equilibrium or QSE groups. This happens e.g. during early or late phases of explosive burning, before or after conditions for a full NSE have been fulfilled (in the latter case this is referred to as "freeze-out").
3.2 Charged-particle freeze-out from NSE and the resulting n/seed ratio In general, r-process environment pass initially through a full NSE, but depending on the initial entropy (combining density ρ and temperature T ), a so-called normal freeze-out (ending with the NSE composition for the freeze-out temperature) or an alpha-rich freeze-out takes place, when during the decline from high temperatures not all alpha-particles can be incorporated into heavy nuclei on the small timescale before the freeze-out of reactions takes place [see Fig. 11 in [17]]. The requirement for an r-process to proceed afterwards covers neutron/seed ratios of 10 to 150 after freeze-out of charged particle reactions in order to produce all (from the light to the heaviest) r-process nuclei via neutron capture from seed abundances. This translates for a normal freeze-out into Y e = 0.12−0.3. Alternatively, for a moderate Y e > 0.40 an extremely alpha-rich freeze-out is needed [see for details the discussion in [17]]. Under the latter conditions the large mass fraction in 4 He (with N = Z ) would permit ratios of remaining free neutrons to (small) abundances of heavier seed nuclei, which are sufficiently high to attain r-process conditions. This is shown in detail in Fig. 5. After charged-particle (normal or alpha-rich) freeze-out, QSE-groups of neutron captures and photodisintegrations are formed along the isotopic chains of heavy elements during the working of the r-process, as discussed in the previous section in Eq.(2). This is an approximate behavior, in reality full network calculations with all necessary nuclear input are needed. But such an approximation gives clues to a full understanding of the results. Exemptions from this behavior can occur in the late r-process phases, when also neutron capture reactions and photodisintegrations start to freeze-out from equilibrium and final neutron captures can reshape the abundance curves. The main conclusion from this discussion is that the characteristics of an r-process are determined by the combination of entropy S and Y e (and an expansion timescale for a close to adiabatic expansion from the initial densities and temperatures determining S). For moderate Y e 's only very high entropies can support an r-process which proceeds up to the actinides, while a low Y e -like in neutron star matter -supports a full r-process already for the lowest entropies. The r-process results for various astrophysical Fig. 5 neutron-to-seed ratios (shown as contour lines) resulting in expanding hot plasmas from explosive burning as a function of the electron abundance Y e and the entropy (measured in k b per baryon). A strong r-process, producing the actinides with n/seed of 150, requires for moderate Y e 's of about 0.45, entropies beyond 250 [18]. The alter-native is that for vanishing entropies (right panel), i.e. cold matter like in neutron stars, the n/seed curves turn over and become independent of entropy. Then an n/seed ratio of 150 requires a Y e of about 0.15 or less. Image reproduced with permission from [18], copyright by AAS environments, spanning a possible continuum between these extremes, reflect directly these connections.

What determines Y e ?
This leaves us with the causes for the initial values of entropy and Y e . While in some cases the initial Y e is set by the starting abundance composition like in neutron star ejecta (where their small Y e was obtained previously at high densities by capture of electrons with high Fermi energies on protons and nuclei), in many environments neutrinos play a significant role which dial the initial Y e before the expansion from maximum temperatures and densities. This is due to neutrino and antineutrino capture on neutrons and protons ν e + n → p + e − ,ν e + p → n + e + , where for similar neutrino and antineutrino energies the first reaction wins and leads to an increase of Y e . The equilibrium electron fraction for environments with neutrino-driven outflows reads [19,20] Here, = 1.293 MeV is the neutron-proton mass difference, = E 2 / E relates to neutrino energies, and L to their luminosities. In the approximation on the right we have assumed that neutrino energies are large enough so that, to acceptable accuracy, the terms containing can be neglected. The above given general considerations are the ones which then directly enter into r-process nucleosynthesis for specific environments.

Nuclear properties
In addition to the astrophysical environments, determining these starting conditions, the nuclear physics aspects (nuclear mass models far from stability, β-decay half-lives and delayed neutron emission, fission, neutron captures and photodisintegrations determining the final abundances during neutron freeze-out, direct capture contributions, α-decays of the final abundances) enter decisively in producing the abundance pattern of the r-process. Major aspects are related to experimental progress in accessing unstable nuclei far from stability, combined with a growing theoretical understanding of their properties. Without giving here a substantial number of related references, these issues are discussed in extended detail in a recent Reviews of Modern Physics article [17] and from experimental as well as theoretical points of view in the chapters by Ani Aprahamian and by Gabriel Martinez-Pinedo and Karlheinz Langanke in the present volume.

Neutrino-driven core-collapse supernovae
The present volume does not contain a specific chapter on core-collapse supernovae and their explosion mechanism, driven by energy deposition going back to neutrinos streaming out of the collapsed proto-neutron star. We refer here to recent review articles [21,22]. The explodability and an overview on final results of the core collapse of massive stars, ending in regular core-collapse supernovae and neutron stars or in black holes has been given e.g. in [23]. The neutrino and antineutrino capture on neutrons and protons is the key ingredient for energy deposition and the supernova explosion mechanism. Thus, combined with the impact on Y e , neutrinos -as discussed in the previous section -play a central role in determining the conditions of the innermost ejecta of core-collapse supernovae for their nucleosynthesis (i.e. the three parameters electron fraction Y e , entropy S, and expansion time scale τ ) from the dynamics of the explosion. Mass accreted through the (standing) shock during the early phase of the explosion gets dissociated into nucleons and alpha particles, depending on the temperature of the shock. This material can flow down to the proto-neutron star or stay in the neutrino heated region and expands. In spherically symmetric simulations, the accretion takes place only in the first second after bounce while in multidimensional simulations downflows stay active during several seconds, providing a continuous nucleosynthesis source. In general, the matter ejected comes from two regions: the surface of the proto-neutron star and the layers which have experiences the outgoing shock. The outer layers of the proto-neutron star and the matter that has been accreted on its surface can be ejected by neutrinos. When this ejection occurs after accretion flows are terminated and the outflow becomes supersonic, this is known as a neutrino-driven wind, suggested initially as the r-process site [see e.g. [19,24] for reviews and references]. Thus, the neutrinos do not only deposit the necessary energy to unbind and eject material, they can also change Y e and change neutrons into protons when neutrino capture wins over antineutrino capture, preventing this way an r-process to occur.
The nucleosynthesis follows the pathway described in the previous section: early on an NSE is attained until the charged-particle freeze-out leads to a seed distribution and the neutron-to-seed (Y n /Y seed ) or possibly a proton-to-seed (Y p /Y seed ) ratio, depending on Y e being smaller or larger than 0.5. At temperatures around 3×10 9 K three processes are possible: (i) a (strong) r-process [18,20,26,27] if Y n /Y seed > 100 [which does not generally occur in neutrino-driven supernovae, opposite to what was expected initially [28,29]], (ii) a weak r-process [24,25,30,31], if Y n /Y seed ∼ 1, or (iii) a νpprocess [e.g. [25,32]], if Y n /Y seed is very small and Y p > Y n , i.e. for proton-rich conditions. In the latter case a proton-rich freeze-out leads to a proton-capture path beyond 56 Ni up to 64 Ge, where a long beta-decay half-life would prevent further build-up of heavier nuclei. However, antineutrino capture on protons provides some neutrons which permit an (n, p)reaction, faster than the beta-decay, and a further build-up. (ii) and (iii) can be recognized in Fig. 6.

Magneto-rotationally driven supernovae
In addition to neutrino-driven explosions, observations of very energetic supernovae [33], long gamma-ray bursts (GRB) [34], pulsars with extremely high magnetic fields [magnetars [35,36]] indicate the key role of magnetic fields in some explosions, here summarized under the name of so-called magneto-rotational supernovae (MR-SN). There are very energetic explosion known as hypernovae (HN) [37] with nucleosynthesis features different from regular core-collapse supernovae [38] that reach energies of about ∼ 10 52 erg which cannot be explained by the neutrino-driven mechanism. These HNe are associated to the collapse of massive stars with > 30 − 40 M , which probably lead to blackhole formation. In all these cases magnetic fields and rotation play a dominant role in the explosion mechanism, which can end either with a central highly magnetized neutron star (magnetar) or alternatively a central stellar-mass black hole. In this subsection we will concentrate on the first subclass of magneto-rotational supernovae, while in the following subsection we will focus on the second case with black hole formation (collapsars or also GRB supernovae), often in the literature identified with hypernovae.
The magneto-rotational mechanism, proposed in the 1970 s, relies on the impact of rotation and magnetic fields on the core-collapse and explosion phase. MR-SNe were proposed by [39] and later by [40] as an r-process site. First simulations, although ignoring or simplifying the neutrino treatment, found a successful r-process in 2D [41] and 3D [42]. Rapid rotation of the iron core is necessary, which leads to an amplification of the magnetic field by rotational winding, being possibly aided by magneto-rotational instabilites (MRI) [43]. After bounce, the strong magnetic pressure launches jets along the rotational axis [e.g. [31,42,44,45]]. However, some 3D models lead to kink instabilities and the formation of less collimated bipolar jets [46,47]. This is the case if the initial pre-collapse magnetic fields are not sufficiently high and the winding of magnetic fields during the collapse does not support an immediate ejection of neutronrich matter along the rotation axis. In this case neutrino interactions enhance Y e before the MRI-effects lead to enhanced magnetic field strengths up to a value which permits ejection and explosion. Based on this, the initial results by [42], starting with pre-collapse magnetic fields of 10 12 G and leading to a full strong r-process up to the actinides, might have been Fig. 6 Innermost ejecta composition from a 16 M core-collapse supernova explosion, resulting in moderately low Y e -matter from collapse and moderately proton-rich conditions due to neutrino interactions during the explosion [shown for simulations with different equations of state [25]]. This leads to a weak r-process and a νp-process with abundance distributions shown in the insert. Image reproduced with permission from [25], copyright by the authors too optimistic for typical conditions as presented in [31] or [48].
A fully self-consistent treatment of these processes requires high resolution simulations which can resolve magnetorotational instabilities (MRI) and predict reliably the possible amplificaton of magnetic fields during the explosion. Present calculations depend on the quality of high-resolution 3D magneohydrodynamics plus the uncertain and therefore assumed initial conditions [49], that either cause strong jet ejection or can develop kink instabilities of the jets. Based on initial conditions, somewhat parametrized in terms of the impact, either neutrino heating or magnetic pressure is causing the supernova explosion. This is related to the question whether very high magnetic fields exist from the beginning of the collapse or are only later developing due the the MRI instability. The latter delays the polar ejection of matter, exposing it longer to neutrinos from the collapse and reducing the neutron-richness of matter and the strength of the r-process. Thus, the production for heavy neutron capture elements varies strongly, being either Fe and Zn dominated like in regular core-collapse SNe or Eu dominated, indicating a strong r-process. This is shown in Fig. 7.
Our present conclusions would be that magneto-rotational supernovae which lead to neutron stars will on average only cause a relatively weak r-process, but further investigations are needed [49].

Collapsars leading to GRB supernovae/hypernovae
Other potential r-process sites associated with a magnetorotational core collapse are due to the accretion discs that form surrounding a massive black hole, originating from the collapse of a massive star with high rotation (a collapsar). Pioneering nucleosynthesis studies [50] have demonstrated Fig. 7 Nucleosynthesic features of rotating core-collapse SN models (h, i−, i, i+, m) with varying ratios of neutrino luminosity and magnetic field strengths. Model m represents a strong MHD-jet supernova. One can see the transition from a regular core-collapse SN pattern, dominated by 56 Ni, total Fe (after decay), and Zn, to a strong r-process pattern with a high Eu abundance. Image reproduced with permission from [31], copyright by AAS that also here neutrinos can play a critical role. While the possibilities for an r-process in recent studies are not conclusive, yet, we discuss here this option in further detail.
Opposite to the above early studies, an interesting question is related to whether some environments can actually cause an actinide boost, i.e. a very strong r-process with large actinide abundance, while others produce a normal solar-type r-process distribution. Recent studies [51] conclude that actinides are substantially overproduced relative to lanthanides for Y e -values in the range 0.1-0.15, sensitive Fig. 8 Abundances predictions by Wu as a function of the mass number A at 10 9 s after an r-process event for different initial Y e 's, focusing on low Y e values and their influence on the abundances of actinides. Image reproduced with permission from [53], copyright by the authors to the influence of fission cycling. This is consistent with other studies [52,53, see Fig. 8], which find, for a variety of nuclear mass models, that electron fractions in the range of ∼ 0.15 are most favorable to explain "actinide boost" matter. Given that "actinide boost" compositions depend on a dominant fraction of the ejecta to originate from a very narrow Y e -range in order to reproduce their observed abundance pattern, requires that nature produces a restricted range of conditions where such Y e -values occur. While to date this question is not settled, we want to discuss here a possibility, recently suggested in [54], based on the equilibrium electron fractions Y e presented already in Eq. (5).
If a black hole is present or forms, the neutrino irradiation is dramatically quenched and the gas flow around the BH regulates itself into a state of mild electron degeneracy. Due to negative feedback between electron degeneracy and neutrino cooling (higher degeneracy leads to fewer electrons and positrons, therefore reducing the neutrino emission, this leads to a higher temperature and thus to a lowering of the degeneracy), the disk midplane settles inside of the inner ∼ 10 G M BH /c 2 to electron fractions of Y e ≈ 0.1 [55]. Interestingly, this occurs once the accretion rates exceed an "ignition value" that depends on the BH spin [56] and the corresponding accretion rates are those that are needed to power (long or short) GRBs [57]. While discovered in semianalytic models, this self-regularization to low Y e -values in the disk midplane is also found in full-fledged numerical (magneto-)hydrodynamic simulations [see e.g. [58][59][60]].
The simulation of such neutrino-cooled accretion flows is a major challenge since one needs to resolve the small length scales of the magneto-rotational instability [61] and to evolve the black hole torus system for a very large number of dynamical time scales (up to several seconds, while the dynamical time scales are ∼ ms). Therefore, it is not entirely surprising that the exploration of this topic is still in initial stages, that large parts of the relevant parameter space are not explored yet and, where parameters are comparable, the results do not yet agree (at least not concerning the ejecta composition). The currently existing GRMHD explorations [58,59,62,63] agree that a large fraction (∼ 40%) of the initial torus mass becomes unbound, but to date there is no agreement about the resulting Y e and composition of the ejecta. For example, [59] find Y e values around 0.12, those of [58] peak around ∼ 0.14, while [63] find a broad distribution between 0.2 and 0.4. Despite a still missing consensus about the ejecta properties, we think that the black hole torus idea for the source of "actinide boost" material is compelling.
The progenitor systems of actinide boost material could then come from massive accretion disks around black holes, forming substantial tori (either from low mass black holes or large BH spins) and, potentially, also collapsar accretion disks. If black hole torus systems indeed manage to eject matter with properties similar to what they produce robustly in their inner torus regions, and a black hole is needed to launch an (either long or short) GRB (rather than, say, a magnetized neutron star), then it would be the GRB engines that produce the "actinide boost" matter [64].

Neutron star and neutron star black hole mergers
A number of proposals for producing the heaviest (r-process) nuclei in Nature have come forward over the years. These include especially neutron star mergers, going back to early suggestions [65][66][67], and concrete predictions for resulting abundance features [see e.g. [68], for a review before the observation of GW170817]. The follow-up of the gravitational wave event GW170817 [69] revealed strong electromagnetic emission in the aftermath of the merger and showed in particular the expected signatures of an r-process powered kilonova. The decay of its bolometric lightcurve agreed well with the expectations for radioactive heating rates from a broad range of r-process elements [70][71][72][73]. Observations of the 2017 kilonova showed that it was initially blue, indicating the synthesis of elements lighter than barium (Z = 56) [74,75]. These elements have low densities of atomic levels and make the medium less opaque, allowing light to escape (decouple from matter) earlier and without having lost too much energy. The presence of lighter r-process elements was also confirmed by the direct observation of strontium (Z = 38) in the spectra [76]. After a few days, the light turned from blue to red, pointing to the presence of lanthanides and actinides [further reviews are given in [17,73,[77][78][79][80]].
Based on these considerations, three components of neutron star merger ejecta contribute to the overall nucleosynthesis site: (i) dynamical ejecta including compressed and shock heated material from the initial collision as well as possible -cold -tidal spiral arm-type ejecta, (ii) winds driven by neutrinos (emitted from the central hot very massive neutron star -until black hole formation -and the accretion disk) and Fig. 9 Abundance comparison between BNS merger models and solar r-process abundances: (left) for dynamical ejecta (i) based on [81], (right) a combination of tidally dominated dynamical ejecta [82] and abundance yields from neutrino-driven wind ejecta (ii) for different massive NS lifetimes until black hole formation after the merger [83]. In the left panel dark gray and brown lines correspond to cold r-process conditions from slow tidal ejecta and to hot r-process conditions from shock-heated fast ejecta. Averaged lines for all trajectories are marked in red, for the fast ejecta in blue. Image reproduced with permission from [81], copyright by APS (left) and from [83], copyright by AAS (right) Image reproduced with permission from [84], copyright by Oxford University Press potentially also by magnetic fields, and (iii) finally mass outflow from the black hole accretion disk. A common feature of these scenarios is that matter reaches NSE with Y e given by weak reactions or in the cold dynamical ejecta component by beta equilibrium in the cold neutron stars before the merger.
Combining the three types of ejecta (and their nucleosynthesis features) discussed above with kilonova lightcurve and spectra observations, there is no doubt that neutron star mergers are indeed a major r-process source. The blue emission points to the production of a light (lanthanide-free) r-process in dynamical and wind ejecta of Fig. 9, affected by neutrino interactions enhancing the related Y e of (i) and (ii), while the late red emission is the natural outcome for heavy (lanthanides and beyond) r-process ejecta from accretion disk outflows with low Y e (iii, see Fig. 10). The r-process outcome as a function of Y e is discussed in detail in refs. [83][84][85][86][87][88][89][90][91]. In summary, there is strong evidence that this neutron star merger event has produced at least a broad, and maybe the whole, r-process range. However, based on the observed lanthanide fraction X La , [92] find that at least for the neutron star merger GW170817 this does not represent a typical solar r-process pattern.

Galactic chemical evolution of neutron capture elements
We have gone in this article through a number of astrophysical sites responsible for neutron capture elements, especially concentrating on the r-process. There exist many observations of individual objects, which hint at their ejecta compositions, from supernovae and their remnants, from neutron star mergers, and from indirect investigations via meteoritic inclusions, just to name a few. However, there exists also an integrated view, how all these additions changed the composition of the interstellar gas out of which the next genera-  [95][96][97][98][99][100]]. The production of the heavier elements, i.e. the s-process and specifically the r-process, has been the focus of this contribution. The role of the weak or main s-process in massive or low and intermediate mass stars has been discussed by other authors in this volume. The overall relation of the sprocess vs. r-process impact can be viewed as a function of metallicity. The dominant s-process contribution comes from low and intermediate mass stars with long evolution-ary lifetimes, therefore entering late in galactic evolution. A figure of observed La/Eu ratios as a function of metallicity in [101] underlines the general tendency that the abundance ratio is at low metallicities dominated by the r-process contribution, while the combination with the delayed s-process leads finally to solar ratios. This supports the understanding that at the lowest metallicities (for [Fe/H]<−2.5) abundance features represent solar r-process ratios, before during the later evolution also s-contributions enter, leading eventually to the solar abundance ratios.
The possible sites of r-process nucleosynthesis have been presented in the previous section, including neutron star mergers, magneto-rotational jet supernovae and collapsars. Quark-deconfinement (QD) supernovae [102] have recently also been suggested as another weak r-process site. We will now have a look at their impact in galactic evolution, examining especially Eu (a typical r-process element, strongly dominated by the r-process contribution) and Sr (pointing to weak r-process contributions). Figure 11 shows, dependent on the [Eu/Fe]-ratio categorized into limited-r and r-enriched r-I and r-II stars, that the observations of low-metallicity stars proof the existence of a weak or limited r-process, while most stars show rprocess enhanced abundance patterns. This goes together with a variation of e.g. the Sr/Eu ratio, ranging from about 1120 down to 0.5 [110], indicating the different decline of the abundance curve as a function of A. Some of the r-process enriched stars show an "actinide boost", i.e. their Th or U to Eu ratio is supersolar [found in most cases for r-II stars, [51,106,107,109,114]].
In summary: We have discussed a number of suggested rprocess sites, but only one of them is proven by a direct observation of the explosive event. Observations of low metallicity stars show essentially three types of patterns, a weak or limited r-process, a strong solar-type r-process, and an actinide-boosted r-process [for reviews and impact on galactic chemical evolution see e.g. [17,54,80,[115][116][117][118][119]]. There exist strong indications that the weak r-process pattern is related to core-collapse supernovae (probably mostly due to magneto-rotational version?). Whether the latter two types are produced in different sites or result from variations within the same site (e.g. neutron star mergers) is still debated. The question is now how to identify features which can point back to individual sites. In terms of chemical evolution models it is important to utilize methods that can treat the impact of rare events and test whether the sizable spread in abundance ratios in Fig. 11 is due to different origins or just the result of an inhomogeneous environment, going back to rare events in the early galaxy.
A promising approach is to look for correlations among different elements, especially in relation to Fe [54,120]. Recent data from the SAGA and JINA databases [103,121], presented in Fig. 11, permit a correlation between Eu and Fe for [Eu/Fe]<0-0.3, i.e. for stars with lower than average r-process enrichment which also show high [Sr/Fe] values. Interpreted in a straight-forward way this would point to a negligible Fe/Eu ratio (in comparison to solar ratios) in the major r-process sources (like e.g. neutron star mergers and collapsars), while a noticeable co-production of Fe with Eu is possible in less strong r-process sources, e.g. with a weak r-process (probably pointing to rare supernova types). The latter options have been discussed in more detail in [54], possibly linking limited-r stars to (rare) supernovae, r-I stars mainly to neutron star mergers, and r-II stars to collapsars. Alternative views, trying to relate the entire r-process abundances to neutron star mergers alone have been presented recently [114,[122][123][124]. Such a type of interpretation requires further and continuous consistency checks through observations, modelling of the contributing objects, and the treatment of their impact in galactic evolution.
Funding Information Open access funding provided by University of Basel.
Data availability statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: This manuscript has the style of a review and does not provide data predictions or observational results. For predictions see the cited literature, for observations the SAGA and JINA databases [103,121].] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.