Sensitivity studies for the weak r process : neutron capture rates

Rapid neutron capture nucleosynthesis involves thousands of nuclear species far from stability, whose nuclear properties need to be understood in order to accurately predict nucleosynthetic outcomes. Recently sensitivity studies have provided a deeper understanding of how the r process proceeds and have identified pieces of nuclear data of interest for further experimental or theoretical study. A key result of these studies has been to point out the importance of individual neutron capture rates in setting the final r-process abundance pattern for a ‘main’ (A ∼ 130 peak and above) r process. Here we examine neutron capture in the context of a ‘weak’ r process that forms primarily the A ∼ 80 r-process abundance peak. We identify the astrophysical conditions required to produce this peak region through weak r-processing and point out the neutron capture rates that most strongly influence the final abundance pattern.


I. INTRODUCTION
In the 1950s it was recognized that the solar system abundances of nuclei heavier than iron could be divided roughly in half based on the nucleosynthesis processes that create them. 1,2 low neutron capture process, or s-process, nuclei lie along the middle of the valley of stability, and rapid neutron capture process, or r-process, nuclei are found on the neutron-rich side of stability.A third process, the p process, is responsible for the significantly less abundant nuclei on the proton-rich side of stability.Since that time much progress has been made, e.g.Ref. 3, and the basic mechanisms of and astrophysical sites for the creation of the s-process 4 and heavy p-process 5,6 nuclei are on a firm footing.The largest remaining open questions involve the site of the synthesis of the heaviest elements (A > 130) in the r process, and the uncertain origins of the elements above the iron peak (80 < A < 130).
The solar r-process abundance pattern traditionally consists of three main abundance peaks around A ∼ 80, 130, and 195, associated with the N = 50, 82, and 126 closed neutron shells.The locations of the peaks in A show that they must have formed in neutron-rich conditions, with the higher-mass peaks requiring greater neutron-richness than the lighter-mass peak.Thus the astrophysical conditions required to produce a good match to the A ∼ 80 solar abundance pattern are distinct from those required to make the heavier peaks. 7Observational data also points to separate origins for these nuclei, with the first hint coming from meteoritic data. 8More recently, spectroscopic studies of low metallicity halo stars have been performed that show for the most part a solar r-process a Electronic mail: surmanr@union.edupattern for elements from barium to lead ('main' r process) but large scatter below barium. 9It may be that these light elements owe their origins to a range of nucleosynthesis processes, some of which may even be proton-rich; the question is difficult to answer without isotopic information for the halo star abundances (see, e.g., Ref. 10 and references therein).Given this uncertainty, the production mechanism for these elements has often been dubbed the light element primary process, or LEPP.
Studies of A ∼ 80 nucleosynthesis are additionally complicated by large uncertainties in nuclear physics data of the short-lived, exotic nuclei that assumably participate in the processes.These nuclei are difficult to measure in the laboratory owing to the small production cross sections and short lifetimes.Nuclei that participate in the r-process typically have half lives that are too short to allow them to be made into a target (i.e. less than at least a month).As neutron targets are not available, neutron capture experiments cannot be performed on these nuclei.Nuclear data relevant to r-process calculations (J π assignments, excitation energies, spectroscopic factors) can however be extracted from transfer reactions, such as (d, p).2][13][14][15][16][17] These measurements were all performed at the Holifield radioactive ion beam facility, 18 which has since ceased operations as a user facility.It is expected that research in this area will resume once the appropriate capabilities are developed at current and next generation facilities.Even then, many interesting nuclear species will still lie outside the reach of anticipated experimental investigations and so their properties need to be calculated using nuclear models.
0][21][22][23] A sensitivity study begins with the choice of a baseline astrophysical simulation.Then an individual nuclear parameter, such as mass, β-decay rate or neutron capture rate, is varied, the simulation is repeated, and the final abundance patterns are compared.This process is repeated for each parameter of interest in the network.While sensitivity studies for masses, β-decay rates, and neutron capture rates have been performed for a main r process, 80 < A < 100 nucleosynthesis has only been addressed on the proton-rich side 24 and for a weak s process. 25The importance of neutron capture rates for a main r process has also been discussed in Refs.26-30.Here we extend the neutron capture rate sensitivity studies of Ref. 20 and 21 to nuclei below the A ∼ 130 r-process peak.While we recognize that these nuclei may have multiple nucleosynthetic origins, here we will address only a potential weak r process-a rapid neutron capture process that forms a solar-type A ∼ 80 r-process peak and potentially nuclei up to the A ∼ 130 peak.Though neutron capture rates may play a role in a proton-rich LEPP, we will not discuss it here.In Section II, we identify astrophysical conditions appropriate for a weak r process.In Section III, we describe our sensitivity studies and review the mechanisms by which individual neutron capture rates in the region of a major r-process peak can influence the final abundance pattern.The compiled results identifying the neutron capture rates with the greatest impact on the final isotopic abundances from a weak r process are presented in Section IV.

II. ASTROPHYSICAL CONDITIONS
We study a weak r process in the context of neutron-rich supernova and collapsar accretion disk winds. 31,32 ind trajectories are characterized by a measure of heating or entropy per baryon, s/k, the dynamic timescale, τ , that describes the drop in temperature and density, e.g., ρ ∼ e −t/τ , and the neutron-richness of the ejected material or electron fraction, Y e .In order to perform a weak r-process sensitivity study, we first need to identify astrophysical conditions that produce these nuclei.We ran many simulations, varying the wind parameters in discrete units of 1 for s/k, 5 ms for τ , and 0.05 for Y e , as in Ref. 33.To assess how well an abundance pattern matches the solar A = 80 data we calculate the deviation, where abundances, Y(A) = X(A)/A, are computed from the atomic mass, A, and the mass fraction along an isobaric chain, X(A) = A=Z+N X(Z, N) (note mass fractions for any simulation will sum to 1, A X(A) = 1).Each abundance curve is scaled to solar by a factor α which is computed by solving: A slice of this astrophysical parameter space at constant Y e = 0.4 is shown in Fig. 1.The best matches to the solar pattern (small deviations) occur between s/k = 30 and s/k = 40 with timescales longer than τ = 0.15 s.Small deviations also occur outside the range of the parameter space shown in Fig. 1 at constant Y e = 0.4, however we rule these conditions out because they overproduce abundances of the heavier elements.
A sample best fit abundance pattern from the Y e = 0.4 grid of astrophysical parameter space is shown in Fig. 2 (case 3 from Table I, short dashed green line) compared to the solar pattern.Note the A ∼ 80 peak is reproduced but element synthesis does not proceed beyond it.Heavier nuclei can be made at a slightly lower Y e , but this has the additional effect of shifting the peak location to higher A, weakening the match to the solar pattern.An example of this effect is shown in Fig. 2 (case 2 from Table I, long dashed blue line).It is a general feature that becomes even more pronounced at higher entropies, which can lead to the so-called A ∼ 90 overproduction problem in slightly neutron-rich supernova winds. 10,34 he best match to the solar pattern from 70 < A < 110 seems to be found at lower entropies and lower Y e 's; Fig. 2 shows an example of this as well (case 1 from Table I, solid red line).
In addition to general parameterized trajectories, we also consider trajectories from the supernova neutrino-driven wind simulations used in Ref. 10 and from the black hole-neutron star merger FIG. 2. Shows abundance patterns for three sample weak r-process trajectories, with wind conditions listed in Table I for case 1 (red solid line), case 2 (long dashed blue line), and case 3 (short dashed green line).The scaled solar r-process residuals (crosses) from Ref. 35

III. SENSITIVITY STUDY
The sensitivity studies took the baseline weak r-process simulations described above as starting points.For each astrophysical trajectory chosen, we determined the baseline abundance pattern with a nuclear network calculation.The nuclear network code employed includes all relevant two-and three-body charged particle reactions, weak reactions, neutron captures, and photodissociations, with reaction rates from JINA REACLIB v1.0. 37Note all of the REACLIB neutron capture rates relevant here are theoretical values, from NON-SMOKER. 38We followed the composition from free nucleons and alpha particles in nuclear statistical equilibrium (NSE) through the assembly of seed nuclei by charged particle reactions and the subsequent weak r process.The reaction network is described in detail in Ref. 39.
Once the baseline abundance pattern was determined, the neutron capture rates of approximately 300 nuclei in the A ∼ 80 region and above were individually varied by a constant factor.Theoretical compilations of neutron capture rates can differ from each other by orders of magnitude, 19,20,27   the factor was taken to be 100.With each individual variation, we reran the simulation, and compared the resulting abundance pattern to the baseline pattern using the sensitivity measure F as in Ref. 20: where X baseline (A) are the final mass fractions of the baseline abundance pattern and X(A) are the final mass fractions for the simulation with the rate change.The outcome of each sensitivity study is thus a set of sensitivity measures F for each of the 300 nuclei whose capture rates were varied.The resulting sensitivity measures from a single sensitivity study are shown in the bottom panel of Fig. 3 for an astrophysical trajectory with entropy s/k = 10, timescale τ = 0.2 s, and electron fraction Y e = 0.3 (case 1 from Table I).These conditions are similar to those that may be obtained in the winds from the accretion disk formed in a compact object merger. 36,40 he final abundance pattern for this baseline simulation is a good match to the solar A ∼ 80 region and above and as shown in Fig. 2, plotted with the red line.This case is a 'true' weak r process, in that we see the establishment of (n, γ )-(γ , n) equilibrium, and the behavior of the nuclear flow near and through the N = 50 closed shell is similar to that in the N = 82 and N = 126 regions in a main r process, with the r-process path as shown in the top panel of Fig. 3. Thus, we expect the role of neutron capture rates to be similar to that identified for the N = 82 region in a main r process, as described in detail in Refs.19 and 20.Indeed, the pattern of sensitivity measures for this example follow the N = 82 region results, with the greatest sensitivity measures F found for nuclei above and to the left of the closed shell, along the β-decay pathways of the closed shell nuclei, as seen in Fig. 3.
Refs. 19 and 20 identified two mechanisms by which neutron capture rates in the A ∼ 130 region influence the final r-process abundance pattern: an early-freezeout photodissociation effect and a All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported license.See: http://creativecommons.org/licenses/by/3.0/Downloaded to IP: 71.193.127.205On: Sun, 02 Nov 2014 20:14:00 late-freezeout neutron capture effect.Both operate toward the end of the r-process, after the neutronto-seed ratio drops below one and (n, γ )-(γ , n) equilibrium begins to fail.If the conditions are still hot at this point, the r-process path can begin to move toward stability by photodissociations in addition to β decays.Individual photodissociation rates can therefore become important, particularly for long-lived, highly populated even-N nuclei.An increase of the photodissociation rate of a highly populated nucleus, for example, can enhance the shift of material to lower A, providing a burst of neutrons that can then be captured by nuclei throughout the network.The late-freezeout neutron capture effect operates later, when the path moves toward stability primarily via β decays and there is competition for the capture of the final few free neutrons.The most important capture rates at this time are for the nuclei along the β-decay pathways of the closed shell nuclei, since these nuclei are so abundant that an increase to the neutron capture rate, for example, can effectively steal neutrons away from the rest of the pattern.
These same two effects operate in this weak r-process trajectory, as we show here for two example nuclei in the A ∼ 80 region in Figs. 3 and 4. The nuclear flow from the equilibrium path nucleus 82 Zn to 80 Zn as the path begins moving toward stability at the onset of freezeout is controlled by the photodissociation rate of 82 Zn, which depends on the capture rate of 81 Zn by detailed balance: where T is the temperature, S n is the one-neutron separation energy, and σ v is the thermallyaveraged neutron capture rate of the adjacent isotope.An increase in the capture rate of 81 Zn by a factor of 100 thus also increases the photodissociation rate of 82 Zn by the same factor.Here this shifts material from 82 Zn to 80 Zn (via 81 Zn, which immediately photodissociates to 80 Zn) at the beginning of freezeout and results in an associated release of neutrons back into the environment.This is illustrated in the top panel of Fig. 4, which shows the change in the rate at which neutrons are captured in the 70 < A < 90 region in the simulation with the rate change compared to the baseline simulation.It is negative since the photodissociations result in a net release of neutrons.These neutrons are captured throughout the rest of the network, as also shown in Fig. 4, which subsequently produces changes throughout the abundance pattern.The second example is of a latefreezeout neutron capture effect, for another isotope of zinc, 78 Zn.This nucleus is populated in late freezeout from the β-decay chain 78 Ni → 78 Cu → 78 Zn.When its capture rate is increased, 78 Zn captures more of the few remaining free neutrons than in the baseline simulation, which leaves fewer neutrons for capture elsewhere in the network.This is shown in the bottom panel of Fig. 4. The element synthesis in many of the trajectories studied does not proceed by a clean weak r process, however.At moderately higher entropies and initial electron fractions, the nucleosynthetic flow through the N = 50 closed shell primarily proceeds via (n, p) and (n, α) reactions. 10,34 hus the flow has moved beyond the N = 50 closed shell by the time charged particle reactions terminate and neutron capture takes over.This is the case for the simulation with s/k = 30, τ = 0.1 s, and Y e = 0.35 (case 2 from Table I), with final abundances plotted with the long dashed blue line in Fig. 2 and abundances at a temperature of T = 1.75 GK shown in the top panel of Fig. 5.There is sufficient r-processing of the material in this example for the early freezeout photodissociation effect and the late freezeout neutron capture effect to operate.However, since the flow has largely moved past N = 50, the nuclei with the highest sensitivity measures F are also found at higher Z, N. As shown in the bottom panel of Fig. 5, the highest sensitivity measures are found primarily for nuclei which lie along the β-decay pathways of the most abundant nuclei along the r-process path with 30 < Z < 34.
The final set of trajectories considered have even higher entropies and electron fractions and produce primarily A ∼ 90 − 100 nuclei.In these cases there is essentially no r-processing of the material, as all of the neutrons have been expended in (n, p) and (n, α) reactions by the time the temperature drops to T ∼ 2 − 3 GK.The simulation with s/k = 35, τ = 0.3 s, and Y e = 0.4 (case 3 from Table I), with abundances given by the short dashed green line in Fig. 2, is an example of this.The neutron abundance for this simulation is shown as a function of decreasing temperature in Fig. 6, as compared to the two example simulations described above that do show weak r-processing.Without any neutrons available for late time neutron captures, we see essentially no neutron capture rate sensitivity measures F greater than one for these trajectories.

IV. COMBINED RESULTS
The neutron capture rate sensitivity study procedure was repeated for the approximately ninety astrophysical trajectories chosen for this investigation.Of these, 55 show some significant sensitivity to neutron capture rates.The lower entropy 10 < s/k < 30, lower electron fraction Y e < 0.35 trajectories that produce a clean weak r-process show neutron capture rate sensitivity measures similar to Fig. 3. Moderate to higher entropy 30 < s/k < 100 and electron fraction 0.35 < Y e < 0.48 trajectories that show partial r-processing show neutron capture rate sensitivity measures similar to Fig. 5.The results from all 55 sensitivity studies, Fig. 7 and Table II, show the maxiumum sensitivity measures F obtained in the full set of studies.

V. CONCLUSION
We have examined the role of neutron capture of A = 80 nuclei in the context of a weak r process from winds that could occur in supernova or collapsar accretion disks.In the interest of studying neutron capture in these environments we have identified astrophysical conditions which match features of the A = 80 solar abundances.The best matches to solar are low entropy 10 < s/k < 30, low electron fraction Y e < 0.35 outflows while moderate to higher entropy 30 < s/k < 100, larger electron fraction 0.35 < Y e < 0.48 outflows show partial r-processing for these nuclei.
The distribution of influential neutron capture rates in the A = 80 region depends critically on the r-process path which is set by these conditions.Low entropy simulations tend to pile up more material at the N = 50 closed shell and so exhibit important neutron capture rates at lower atomic mass compared to high entropy simulations, see Figures 3 and 5 respectively.
In agreement with previous studies of neutron capture rates in a main r process, we show that neutron capture rates of individual nuclei in a weak r process are significant only once (n, γ )-(γ , n) equilibrium breaks down and the neutron-to-seed ratio drops below unity.We also find both the early time freezeout photodissociation mechanism and the late freezeout neutron capture mechanism are active in this region.
A complete description of the A = 80 abundances we find in nature is fraught with compilations.These light elements may owe their origins to a range of nucleosynthesis processes including the weak r-process models studied here.The resolution of isotopic abundances in low metallicity halo stars will help settle this issue however these observations alone will not be enough.The reduction of large nuclear physics uncertainties, such as those found in neutron capture rates, is critical for accurately predicting final abundance patterns. 41Direct (n, γ ) measurements are not possible with the very short-lived nuclear species pointed out in this work.Indirect techniques such as (d, p) or Coulomb dissociation measurements can be employed, though further development is required to connect this type of experimental data to the thermal cross sections used in r-process calculations.Thus improvements in theoretical models of neutron capture are needed along with the advances in experiment to reduce these key nuclear physics uncertainties.

FIG. 1 .
FIG.1.Shows the calculated deviation, D, from solar for a slice of astrophysical parameter space with Y e = 0.4.Small deviations (blue) denote conditions which best reproduce the A = 80 solar region.The timescale is given in units of milliseconds.

FIG. 3 .
FIG. 3. Sensitivity measures F (bottom panel) resulting from a sensitivity study in which each neutron capture rate in the A ∼ 80 region is varied by a factor of 100.The baseline simulation is case 1 from Table I and Fig. 2. The top panel shows the abundances at the point in the simulation where the neutron-to-seed ratio is one (at the onset of freezeout).The two boxed nuclei, 78 Zn and 81 Zn, are singled out for further analysis.

FIG. 4 .
FIG.4.The resulting change in the rate at which neutrons are captured in the A ∼ 80 region (solid line) and elsewhere (dashed line) when the neutron capture rate of 81 Zn (top panel) or 78 Zn (bottom panel) is increased by a factor of 100, compared to a baseline simulation with no neutron capture rate modifications.

TABLE I .
are included for comparison.Wind conditions for three example weak r-process trajectories.
accretion disk outflows of Ref.36.A total of approximately ninety distinct astrophysical trajectories are chosen as baseline simulations for our neutron capture rate sensitivity studies.Included in this set are trajectories that result in an excellent match to the solar weak r-process pattern as well as those that make mostly A ∼ 90 − 100 nuclei; main r-process trajectories are excluded.

TABLE II .
Nuclei with maximum neutron capture rate sensitivity measures F > 10 from the combined results of fifty-five neutron capture rate sensitivity studies run under a range of distinct astrophysical conditions, from Fig.7.