International Scholarly Research Notices

International Scholarly Research Notices / 2013 / Article

Research Article | Open Access

Volume 2013 |Article ID 170954 |

M. Kiss, Z. Trócsányi, "Phenomenological Description of Neutron Capture Cross Sections at 30 keV", International Scholarly Research Notices, vol. 2013, Article ID 170954, 8 pages, 2013.

Phenomenological Description of Neutron Capture Cross Sections at 30 keV

Academic Editor: H. Dehnen
Received07 Dec 2012
Accepted26 Dec 2012
Published30 Jan 2013


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, [14]). 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 [59].

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.

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.

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.

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.


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

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

Nucleus Nucleus Nucleus

Ac138 15.9 Th141 16.7 Np145 16.5
Pu152 12.3 Am147 22.7 Bk148 12.4
Bk149 6.41 Bk151 7.13 Cf155 20.0
Es155 5.17 Es152 7.48 Es154 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.

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.

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

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.


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.


  1. 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. View at: Publisher Site | Google Scholar
  2. A. G. W. Cameron, “Nuclear astrophysics,” Annual Review of Nuclear and Particle Science, vol. 8, pp. 299–326, 1958. View at: Publisher Site | Google Scholar
  3. 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. View at: Publisher Site | Google Scholar
  4. 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. View at: Publisher Site | Google Scholar
  5. 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. View at: Publisher Site | Google Scholar
  6. 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. View at: Publisher Site | Google Scholar
  7. 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. View at: Publisher Site | Google Scholar
  8. T. Nakagawa, S. Chiba, T. Hayakawa, and T. Kajino, “Maxwellian-averaged neutron-induced reaction cross sections and astrophysical reaction rates for kT=1 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. View at: Publisher Site | Google Scholar
  9. 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. View at: Publisher Site | Google Scholar
  10. C. E. Rolfs and W. S. Rodney, Cauldrons in the Cosmos, The University of Chicago Press, Chicago, Ill, USA, 1988.
  11. 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. View at: Google Scholar
  12. 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. View at: Publisher Site | Google Scholar
  13. F.-K. Thielemann, M. Arnould, and J. W. Truran, “Advances in nuclear astrophysics,” in Proceedings of the 2nd IAP Workshop, Paris, France, July 1986. View at: Google Scholar
  15. H. Beer and F. Käppeler, “Neutron capture cross sections on 138Ba, 140,142Ce, 175,176Lu, and 181Ta at 30 keV: prerequisite for investigation of the 176Lu cosmic clock,” Physical Review C, vol. 21, no. 2, pp. 534–544, 1980. View at: Publisher Site | Google Scholar

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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.