ISRN Astronomy and AstrophysicsVolume 2013 (2013), Article ID 170954, 8 pageshttp://dx.doi.org/10.1155/2013/170954

Research Article

Phenomenological Description of Neutron Capture Cross Sections at 30 keV

M. Kiss¹ and Z. Trócsányi²

¹Berze N.J. Gimnázium, Kossuth 33, Gyöngyös 3200, Hungary²Institute of Physics, University of Debrecen, P.O. Box 105, Debrecen 4010, Hungary

Received 7 December 2012; Accepted 26 December 2012

Academic Editors: H. Dehnen and R. Hawkes

Copyright © 2013 M. Kiss and Z. Trócsányi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Studying the published values of Maxwellian-averaged neutron capture cross sections, we found simple phenomenological rules obeyed by the cross sections as a function of proton and neutron numbers. We use these rules to make predictions for cross sections of neutron capture on nuclei with proton number above 83, where very few MACS data are available. These predictions may be useful in certain models of nucleosynthesis of heavy nuclei in stars.

1. Introduction

Theoretical descriptions of nucleosynthesis in stars rely heavily on the knowledge of capture cross sections of slow neutrons on nuclei. The classical model of nucleosynthesis in weak neutron flux is based on slow neutron capture (the process) that occurs along a path in the stability valley of nuclei (see, for instance, [1–4]). The -process evolution codes take into account the most important processes (those with largest cross sections) along the stability valley. The necessary information on the neutron capture cross sections and -decay life times, needed to describe qualitatively the abundances of the -process elements, is rather well known from laboratory experiments [5–9].

The -process model is capable to explain the observed abundance of heavy elements fairly well [10]. The difference of observation and prediction is largely attributed to another process that occurs in stellar environment with high neutron flux, typically in supernovae. In such circumstances, the neutron capture is very likely, and neutron-rich nuclei far from the stability valley build up very quickly due to repeated capture of neutrons. The nuclei produced in such a way are so unstable and short lived that experimental information about their capture cross sections and decay life times is not generally available.

In a recent work we proposed a unified model of nucleosynthesis of heavy elements in stars [11] that interpolates between the process and process smoothly. The approach takes into account all possible types of production and depletion mechanisms and solves the whole system of differential equations numerically. The result of such an approach is that (instead of the -process path) the evolution of the synthesis proceeds along a band in the valley of stable nuclei. The width of this band—and consequently the final abundances of nuclei—depends on the neutron flux and the capture cross sections on individual nuclei characterized by both their proton and neutron numbers, , which constitute an essential input to the model calculations. Therefore, it is important to learn about these cross sections as much as possible.

In this paper, we study the general features of MACS collected in recent compilations [8, 9]. We used both measured and calculated cross sections. According to the Hauser-Feshbach (statistical) model, these cross sections can be computed, provided there are sufficiently many energy levels in the nucleus, which is usually satisfied by medium and heavy nuclei [12]. At 30 keV the calculated and measured cross sections agree within a factor of two [13], which is sufficient for our target precision. We have taken values of MACS from [14].

In Section 2 we show some phenomenological observations. In the following section we use those to make some order of magnitude predictions for the capture cross sections for proton numbers , where only very few MACS data are available. Section 4 contains our conclusions.

2. Observations

Neutron capture cross sections depend on the energy or equivalently the speed of neutrons. At moderate energies [10]; therefore, the averaged reaction rate is approximately constant within a wide range of neutron energies (10–100 keV), relevant for nucleosynthesis of heavy elements in stars. As a result one can choose a convenient value for the temperature, traditionally keV.

There are cases when the capture cross section is due to several narrow resonances, or values at higher temperatures are necessary (for instance in the case of the astrophysical process). Then one has to measure the cross section on a wider energy range (the range 1–300 keV is relevant for slow neutron capture) and compute the average with the Boltzmann energy distribution. In some measurement setup, it was possible to measure MACS directly [15]. These measurements showed that the neutron spectrum can be well approximated by a Maxwellian fit with temperature close to keV [3].

MACS at 30 keV have been calculated for many nuclei and made available in public depositories. Currently MACS values based on experimentally measured cross sections are available for stable and for some long-lived targets. For unstable ones cross sections obtained from theoretical reaction models are used for calculating MACS. A comprehensive and complete review of available MACS values has been presented recently in [9].

Studying the available MACS, we can make several observations: (i) although MACS of many nuclei are available, there are still many missing, or rather uncertain values, especially for nuclei with (see Figure 1 where all MACS, currently available, are shown); (ii) the cross sections vary over very large range of values (about four orders of magnitude); (iii) for any fixed neutron number the cross section is maximal for a corresponding value of the proton number and decreases rapidly as increases (see Figure 2). The last point implies that in the plain for each there is a unique value where the capture cross section attains its maximum value. The quantitative understanding of the existence of such a maximum is not our goal in the present paper.

Figure 1: All presently available MACS (at 30 keV) on nuclei as a function of the proton and neutron numbers.

Figure 2: Dependence on the proton number of the MACS on nuclei with fixed neutron number N = 30, 31, 60, 61, 90, and 91 (indicated by vertical lines on Figure 1).

If we plot the function, a rather simple picture emerges: it appears that a simple, almost linear function can describe the MACS values, especially for small . This feature becomes even more salient if we divide the nuclei into four groups according to the even/odd number of protons and neutrons: (i) for even, even, (ii) for odd, odd, (iii) for even, odd, and (iv) for odd, even, as shown in Figure 3.

(a)

(b)

(c)

(d)

Figure 3: The function for even-even, odd-odd, even-odd, and odd-even nuclei. The crosses are the experimental values, and the solid line represents the fit to the function in (1). The dashed line runs through the bottom of the stability valley.

In Figure 3 crosses mark the values of where the -capture cross section is maximal for a fixed value of the neutron number as taken from [8]. The solid lines represent fits of simple functions to these points in the form of with , , and being fitted parameters, and ee, oo, eo, or oe. We determined the values of these parameters in two steps. First, we minimized the function that is, nuclei belonging to all four groups are taken into account. As the values of can only be integers; therefore, we assumed that the uncertainty of determining is for any . The upper limit of the summation in each group was chosen as the largest value for which can be identified. With such a choice we find , 56, 53, and 55 maxima in the groups of even-even, odd-odd, even-odd, and odd-even nuclei, respectively, . This fit gives with correlation index that is, the coefficient of determination is close to one, .

In the second step we minimize the functions separately for each group ( ee, oo, eo, oe). (Again, we assume for all ). These fits give with uncertainty being 0.16 and coefficient of determination above 0.99 in all cases. In principle the parameters and could also depend on , but such dependence did not show any improvement in the fits; therefore, we kept the fitted formula as universal as possible.

We also exhibit the line of the stability valley in Figure 3 as a function of (instead of the usual ) with parameters We see clearly that the highest -capture cross sections lie above the stability valley, and the separation grows with .

We can also observe regularity in the dependence of the cross section at fixed (see Figure 2). We can extrapolate this regularity as well as the values to the region in the nuclide chart where very few MACS values are available for -capture cross sections on nuclei (nuclei with proton number above 83, see Figure 1).

The first observation is a simple trend in the behaviour of the function . Putting on a double logarithmic plot as shown in Figure 4(a) (left panel), we find that the general trend is well described by a fourth-order power function, This general trend is slightly modulated with some oscillatory behaviour, with minima around magic numbers, as seen on Figure 4(b), where the ratios of the published cross sections [9] to are shown.

(a)

(b)

Figure 4: (a) Largest neutron capture cross sections as a function of the neutron number. (b) Ratio of the published [9] largest cross sections to given in (9).

The second observation is that if we normalize the cross sections for a fixed neutron number with the largest cross section , the profile of the dependence on the proton number is rather similar for all neutron numbers. This similarity is best seen if the position of the largest cross section is shifted by to zero; therefore, we define these normalized and shifted cross section values by for all values of , where MACS values are available. Then we define the average -dependence by with squared standard deviation where is the number of available MACS values for fixed . This average is shown in Figure 5. As seen from Figure 5(b) this function is well approximated with an almost exponential function in both positive and negative directions, but with different exponents. More precisely, we fit the logarithm of the average with quadratic functions of the form with subscript of the coefficients referring to three regions in : (i) for , (ii) for , and (iii) for . For and 3 we fix . This form ensures the constraint . We also require the continuity of the fitted function at . We measure the goodness of the fit by the weighted sum of squares summed over values of in the three regions separately. The result of these fits is presented in Table 1 and shown in Figure 5.

Table 1: Result of the fit to the average function .

/d.o.f

1
0.0044
1.135
17.95
5.29/5

2
−0.0025
0.2658
0
6.15/9

3
−0.0058
−0.3948
0
7.14/4

(a)

(b)

Figure 5: Average of the normalized neutron capture cross sections as a function of . The error bars represent the standard deviation .

Each function differs from the average in two ways: (i) typically the larger the wider (as seen on Figure 1), (ii) in addition there are seemingly random fluctuations. The origin of the latter could be either a small physical effect or simply error of the measurement: there are published values for cross sections that differ by a factor of two. While it is difficult to consider the effect of the latter, the first effect can be taken into account by a simple appropriate scaling of the width of the average to those of the functions , which we discuss in the next section.

3. Predictions

The phenomenological observations made in the previous section can be used to make predictions for the order of magnitude of neutron capture cross sections in regions of the nuclide chart where MACS values based on measured data are not available. We make these predictions in two steps. First we validate our procedure by comparing our predictions to published [9] cross sections. Then we use our procedure to make predictions.

3.1. Procedure

Our procedure relies on three pieces of information concluded from the analysis of the shape of ridge of Maxwellian-averaged neutron capture cross sections: (1)position of as a function of the neutron number (location of the ridge top on the nuclide chart) obeys the simple function (1), (2)values of (height of the ridge for given value of ) obey the simple function (9), (3)characteristic behaviour of the average function (slope of the ridge) is given by Figure 5.

In order to predict the cross section values for fixed neutron number, we proceed along the following steps.(1)Given , find the position of from (1), which gives two maxima, one for even proton numbers () and one for odd proton numbers ().(2)Given (either or ), position the maximum location of the average function to .(3)Scale the height and width of the function to the available MACS values by performing a two-parameter fit: (i) the scale factor of the height and (ii) the scale factor of the width.

The third step is hampered by the discrepancies in the published cross section values, which can sometimes be quite significant as shown in Table 2 for some heavy nuclei. Discrepancies exist among MACS values for lighter nuclei, but generally within a factor of two [9].

Table 2: Ratios of largest and smallest published neutron capture cross sections () for nuclei beyond bismuth with [9].

Nucleus

Nucleus

Nucleus

Ac₁₃₈
15.9
Th₁₄₁
16.7
Np₁₄₅
16.5

Pu₁₅₂
12.3
Am₁₄₇
22.7
Bk₁₄₈
12.4

Bk₁₄₉
6.41
Bk₁₅₁
7.13
Cf₁₅₅
20.0

Es₁₅₅
5.17
Es₁₅₂
7.48
Es₁₅₄
46.8

3.2. Validation

We can compare the values of the predicted cross sections to those calculated MACS [9] over the regions of the nuclide chart where values are abundantly available (). In Figure 6 we show again the cross sections of Figure 2 together with the predicted values following from our procedure described in the previous subsection. Considering the simple nature of our procedure, the agreement is striking for all neutron numbers. Of course, the predictions rarely coincide exactly with the published values, but the order of magnitude is usually correct, especially where the cross sections are large, which is the most important region for nucleosynthesis. Similar agreement can be observed over the large region of the nuclide chart where published MACS are available.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: Dependence on the proton number of MACS (at 30 keV) on nuclei with fixed neutron number = 30, 31, 60, 61, 90, and 91: comparison of the predictions of the phenomenological model to published MACS values.

3.3. Predictions of Unknown Cross Sections

Our procedure can be used to make prediction for cross sections in regions of the nuclide chart where some experimental information are available, such as . In this region the general trend can be fitted to the published MACS to complete the ridge. With such a procedure we obtain cross section values shown in Table 3. We can now use those predictions to complete the picture exhibited on Figure 1. The result of such completion is shown in Figure 7.

Table 3: Some predictions for neutron capture cross sections (in mb) as a function of the proton and neutron numbers for nuclei beyond bismuth. Based on the validation results, we estimate the relative uncertainty of the largest values along isotones below 10%.

127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142

84
6
1
8
2
15
3
28
3
216
53
331
99
121
90
157
82

85
50
16
72
28
152
74
577
99
772
534
229
530
171
385
29
215

86
10
2
12
3
26
6
59
5
330
3
512
156
205
150
259
146

87
76
23
108
48
268
135
1024
163
1276
875
437
792
302
693
74
376

88
16
3
20
5
46
13
118
8
601
227
653
448
343
6
421
258

89
114
34
161
83
463
242
1786
266
2083
1560
2020
1366
523
1229
184
650

90
25
5
32
8
79
26
233
15
751
252
1400
429
1400
433
1550
484

91
171
49
237
139
786
428
3063
428
3359
2269
600
1770
695
2140
1213
2250

92
38
8
50
14
135
51
450
25
1118
412
1790
427
492
770
425
1550

93
253
70
346
229
1314
743
5168
681
5347
3590
2717
2514
1506
3692
600
1020

94
59
13
78
23
225
98
849
41
1648
666
2667
861
1496
1036
1693
750

95
370
100
500
374
2162
1267
8577
1069
8406
5613
4816
3635
2499
6256
2345
3093

96
90
20
120
37
371
184
1568
68
2407
1063
3938
1289
2383
1631
2635
2191

97
209
59
288
600
3500
2124
14000
1657
13048
8675
8379
5212
4087
10442
5132
5065

98
48
30
63
60
600
336
2831
111
3485
1675
5761
1913
3748
2538
4055
3611

99
115
34
162
295
1701
3500
6718
2536
20000
13250
14312
7411
6584
17168
10876
8185

100
25
16
32
30
291
600
5000
178
5000
2605
8353
2813
5820
3904
6172
5873

Figure 7: Ridge of MACS (at 30 KeV) on nuclei as a function of the proton and neutron numbers.

4. Conclusions

We studied the dependence of the published MACS values on the proton and neutron numbers. We found a simple characteristic behaviour that we call the shape of the ridge of MACS in the nuclide chart. This shape can be described by the position and height of the ridge and the decrease of the slope. Although there are local deviations from this general trend, we were able to quantify these characteristics and made predictions for cross sections in regions of the nuclide chart where only few MACS values are available. Such predictions are needed for computer programs aimed at simulating the formation of all (i.e., not only along the -process path) heavy nuclei in stars.

Acknowledgments

This paper was supported by the TÁMOP 4.2.1./B-09/1/KONV-2010-0007 Project. The authors are grateful to I. Angeli for the useful discussions.

References

E. M. Burbidge, G. R. Burbidge, W. A. Fowler, and F. Hoyle, “Synthesis of the elements in stars,” Reviews of Modern Physics, vol. 29, no. 4, pp. 547–650, 1957.

A. G. W. Cameron, “Nuclear astrophysics,” Annual Review of Nuclear and Particle Science, vol. 8, pp. 299–326, 1958.

F. Käppeler, H. Beer, and K. Wisshak, “S-process nucleosynthesis-nuclear physics and the classical model,” Reports on Progress in Physics, vol. 52, no. 8, article 945, 1989.

R. Gallino, C. Arlandini, M. Busso, et al., “Evolution and nucleosynthesis in low-mass asymptotic giant branch stars. II. Neutron capture and the s-process,” The Astrophysical Journal, vol. 497, p. 388, 1998.

T. Rauscher and F. K. Thielemann, “Astrophysical reaction rates from statistical model calculations,” Atomic Data and Nuclear Data Tables, vol. 75, no. 1-2, p. 1, 2000.

Z. Y. Bao, H. Beer, F. Käppeler, F. Voss, K. Wisshak, and T. Rauscher, “Neutron cross sections for nucleosynthesis studies,” Atomic Data and Nuclear Data Tables, vol. 76, no. 1, pp. 70–154, 2000.

I. Dillmann, M. Heil, F. Käppeler, R. Plag, T. Rauscher, and F. K. Thielemann, “KADoNiS—The Karlsruhe Astrophysical Database of Nucleosynthesis in Stars,” AIP Conference Proceedings, vol. 819, pp. 123–127.

T. Nakagawa, S. Chiba, T. Hayakawa, and T. Kajino, “Maxwellian-averaged neutron-induced reaction cross sections and astrophysical reaction rates for keV to 1 MeV calculated from microscopic neutron cross section library JENDL-3.3,” Atomic Data and Nuclear Data Tables, vol. 91, no. 2, pp. 77–186, 2005.

B. Pritychenko, S. F. Mughaghab, and A. A. Sonzogni, “Calculations of Maxwellian-averaged cross sections and astrophysical reaction rates using the ENDF/B-VII.0, JEFF-3.1, JENDL-3.3, and ENDF/B-VI.8 evaluated nuclear reaction data libraries,” Atomic Data and Nuclear Data Tables, vol. 96, no. 6, pp. 645–748, 2010.

C. E. Rolfs and W. S. Rodney, Cauldrons in the Cosmos, The University of Chicago Press, Chicago, Ill, USA, 1988.

M. Kiss and Z. Trócsányi, “A unified model for nucleosynthesis of heavy elements in stars,” Journal of Physics, vol. 4, Article ID 012024, 2010.

K. Langanke, F. K. Thielemann, and M. Wiescher, “Nuclear astrophysicsand nuclei far from stability,” Lecture Notes in Physics, vol. 651, Article ID 383467, pp. 383–467, 2004.

F.-K. Thielemann, M. Arnould, and J. W. Truran, “Advances in nuclear astrophysics,” in Proceedings of the 2nd IAP Workshop, Paris, France, July 1986.

http://adg.llnl.gov/Research/RRSN/semr/30kev/rath00_7.4.30kev_calc

H. Beer and F. Käppeler, “Neutron capture cross sections on ¹³⁸Ba, 140,142Ce, 175,176Lu, and ¹⁸¹Ta at 30 keV: prerequisite for investigation of the ¹⁷⁶Lu cosmic clock,” Physical Review C, vol. 21, no. 2, pp. 534–544, 1980.