# Phenomenological Description of Neutron Capture Cross Sections at 30 keV

Authors: M. Kiss; Z. Trócsányi
Journal: ISRN Astronomy and Astrophysics
 (2013)

---

## 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.


---

## Body

## 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 (thes process) that occurs along a path in the stability valley of nuclei (see, for instance, [1–4]). The s-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 s-process elements, is rather well known from laboratory experiments [5–9].Thes-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 s process and r 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 s-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, σ(Z,N), 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 Section2 we show some phenomenological observations. In the following section we use those to make some order of magnitude predictions for the capture cross sections σ(Z,N) for proton numbers Z>83, 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σ∝1/v [10]; therefore, the averaged reaction rate 〈σv〉 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 kBT=30keV.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 astrophysicalr 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 kBT=30keV [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 withZ>83 (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 N the cross section is maximal for a corresponding value of the proton number Zmax and decreases rapidly as |Z-Zmax| increases (see Figure 2). The last point implies that in the Z-N plain for each N there is a unique value Zmax(N) 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 numberZ 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 theZmax(N) function, a rather simple picture emerges: it appears that a simple, almost linear function can describe the MACS values, especially for small N. 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) Zmax(ee) for Z even, N even, (ii) Zmax(oo)  for Z odd, N odd, (iii) Zmax(eo) for Z even, N odd, and (iv) Zmax(oe) for Z odd, N even, as shown in Figure 3.The functionZmax(N) 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.
(a)
(b)
(c)
(d)In Figure3 crosses mark the values of Zmax where the n-capture cross section is maximal for a fixed value of the neutron number N as taken from [8]. The solid lines represent fits of simple functions to these points in the form of
(1)f(N;ax,b,c)=N+ax1+bNc,
with ax, b, and c being fitted parameters, and x= ee, oo, eo, or oe. We determined the values of these parameters in two steps. First, we minimized the function
(2)χ2(a,b,c)=∑i=1214‍(Zmax(Ni)-f(N;a,b,c)σi)2,
that is, nuclei belonging to all four groups are taken into account. As the values of Zmax(Ni) can only be integers; therefore, we assumed that the uncertainty of determining Zmax(Ni) is σi=1 for any i. The upper limit of the summation in each group was chosen as the largest value for which Zmax can be identified. With such a choice we find Nmax(x)=50, 56, 53, and 55 maxima in the groups of even-even, odd-odd, even-odd, and odd-even nuclei, respectively, (50+56+53+55=214). This fit gives
(3)a=0.60±0.08,b=0.013±0.0001,c=0.666±0.002,
with correlation index
(4)i=1-χ2(a,b,c)∑i=1214(Zmax(Ni)-Z¯)2=1-590.2116029=0.997,
that is, the coefficient of determination is close to one, i2=0.994(Z¯=(1/214)∑i=1214‍Zmax(Ni)=45.5).In the second step we minimize the functions(5)χ2(ax)=∑i=1Nmax(x)‍(Zmax(x)(Ni)-f(N;ax,0.013,0.666)σi)2
separately for each group (x= ee, oo, eo, oe). (Again, we assume σi=1 for all i). These fits give
(6)aee=2.47,aoo=-0.38,aeo=0.45,aoe=0.07,
with uncertainty being 0.16 and coefficient of determination above 0.99 in all cases. In principle the parameters b and c could also depend on x, 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 Figure3 as a function of N (instead of the usual A=Z+N)
(7)Zstab=N+as1+bsNcs,
with parameters
(8)as=0.682,bs=0.027,cs=0.614.
We see clearly that the highest n-capture cross sections lie above the stability valley, and the separation grows with N.We can also observe regularity in theZ dependence of the cross section at fixed N (see Figure 2). We can extrapolate this regularity as well as the Zmax values to the region in the nuclide chart where very few MACS values are available for n-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 σmax(N)≡σ(Zmax(N)). Putting σmax(N) 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,
(9)σmax(N)=(N10)4mb.
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 σmax(N) are shown.(a) Largest neutron capture cross sections as a function of the neutron number. (b) Ratio of the published [9] largest cross sections to σmax given in (9).
(a)
(b)The second observation is that if we normalize the cross sections σ(Z,N) for a fixed neutron number N with the largest cross section σmax(N), 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 -Zmax to zero; therefore, we define these normalized and shifted cross section values by
(10)ρN(z)=σ(z+Zmax,N)σmax(N)≡σ(Z,N)σ(Zmax(N)),z=Z-Zmax,
for all values of N, where MACS values are available. Then we define the average z-dependence by
(11)ρ(z)=1Nz∑N=1Nz‍ρN(z),
with squared standard deviation
(12)σ(z)2=1Nz(Nz-1)∑N=1Nz‍[ρN(z)-ρ(z)]2,
where Nz is the number of available MACS values for fixed z. 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 αiz2+βiz+γi with subscript of the coefficients referring to three regions in z: (i) i=1 for z<-26, (ii) i=2 for -26≤z<0, and (iii) i=3 for z>0. For i=2 and 3 we fix γi=0. This form ensures the constraint ρ(0)=1. We also require the continuity of the fitted function at z=-26. We measure the goodness of the fit by the weighted sum of squares
(13)χ2≃∑z‍[lnρ(z)-(αiz2+βiz+γi)]2((σ(z))/ρ(z))2,
summed over values of z 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ρ(z).
i
α
i
β
i
γ
i
χ
2/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/4Average of the normalized neutron capture cross sections as a function ofz=Z-Zmax. The error bars represent the standard deviation σ(z).
(a)
(b)Each functionρN(z) differs from the average in two ways: (i) typically the larger N the wider ρN(z) (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 σ(Z,N) 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 ρN(z), 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 ofridge of Maxwellian-averaged neutron capture cross sections: (1)
position ofZmax as a function of the neutron number (location of the ridge top on the nuclide chart) obeys the simple function (1),(2)
values ofσmax(N) (height of the ridge for given value of Zmax(N)) obey the simple function (9),(3)
characteristic behaviour of the average functionρ(z) (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)
GivenN, find the position of Zmax from (1), which gives two maxima, one for even proton numbers (Zmax(e)) and one for odd proton numbers (Zmax(o)).(2)
GivenZmax (either Zmax(e) or Zmax(o)), position the maximum location of the average function ρ(z) to Zmax.(3)
Scale the height and width of the functionρ(z) 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 Table2 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 (R=σmax/σmin) for nuclei beyond bismuth with R>5 [9].
Nucleus
R
Nucleus
R
Nucleus
R

89
227Ac138
15.9

90
231Th141
16.7

93
238Np145
16.5

94
246Pu152
12.3

95
242Am147
22.7

97
245Bk148
12.4

97
246Bk149
6.41

97
248Bk151
7.13

98
253Cf155
20.0

99
254Es155
5.17

99
251Es152
7.48

99
253Es154
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 (Z≤82). 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.Dependence on the proton numberZ of MACS (at 30 keV) on nuclei with fixed neutron number N = 30, 31, 60, 61, 90, and 91: comparison of the predictions of the phenomenological model to published MACS values.
(a)
(b)
(c)
(d)
(e)
(f)


### 3.3. Predictions of Unknown Cross Sections

Our procedure can be used to make prediction for cross sections in regions of the nuclide chart wheresome experimental information are available, such as Z>83. 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%.
Z
N
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
5873Figure 7
Ridge of MACS (at 30 KeV) on nuclei as a function of the proton and neutron numbers.


## 3.1. Procedure

Our procedure relies on three pieces of information concluded from the analysis of the shape ofridge of Maxwellian-averaged neutron capture cross sections: (1)
position ofZmax as a function of the neutron number (location of the ridge top on the nuclide chart) obeys the simple function (1),(2)
values ofσmax(N) (height of the ridge for given value of Zmax(N)) obey the simple function (9),(3)
characteristic behaviour of the average functionρ(z) (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)
GivenN, find the position of Zmax from (1), which gives two maxima, one for even proton numbers (Zmax(e)) and one for odd proton numbers (Zmax(o)).(2)
GivenZmax (either Zmax(e) or Zmax(o)), position the maximum location of the average function ρ(z) to Zmax.(3)
Scale the height and width of the functionρ(z) 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 Table2 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 (R=σmax/σmin) for nuclei beyond bismuth with R>5 [9].
Nucleus
R
Nucleus
R
Nucleus
R

89
227Ac138
15.9

90
231Th141
16.7

93
238Np145
16.5

94
246Pu152
12.3

95
242Am147
22.7

97
245Bk148
12.4

97
246Bk149
6.41

97
248Bk151
7.13

98
253Cf155
20.0

99
254Es155
5.17

99
251Es152
7.48

99
253Es154
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 (Z≤82). 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.Dependence on the proton numberZ of MACS (at 30 keV) on nuclei with fixed neutron number N = 30, 31, 60, 61, 90, and 91: comparison of the predictions of the phenomenological model to published MACS values.
(a)
(b)
(c)
(d)
(e)
(f)


## 3.3. Predictions of Unknown Cross Sections

Our procedure can be used to make prediction for cross sections in regions of the nuclide chart wheresome experimental information are available, such as Z>83. 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%.
Z
N
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
5873Figure 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 thes-process path) heavy nuclei in stars.



---

*Source: 170954-2013-01-30.xml*