Neutron capture reactions relevant to s-process and p-process in the domain of the N = 50 shell closure Saumi Dutta,∗ G. Gangopadhyay,† and Abhijit Bhattacharyya‡ Department of Physics, University of Calcutta 92, Acharya Prafulla Chandra Road, Kolkata-700 009, India Radiative thermal neutron capture cross sections for nuclei participating in s-process and p- process nucleosynthesis in and around N = 50 closed neutron shell have been calculated in statis- tical semi-microscopic Hauser-Feshbach approach for the energy range of astrophysical interest. A folded optical-model potential is constructed utilizing the standard DDM3Y real nucleon-nucleon interaction. The folding of the interaction with target radial matter densities, obtained from the relativistic-mean-field approach, is done in coordinate space using the spherical approximation. The standard nuclear reaction code TALYS1.8 is used for cross-section calculation. The cross sections are compared with experimental results and reasonable agreements are found for almost all cases. Maxwellian-averaged cross sections (MACS) for the nuclei are presented at a single thermal energy of 30 keV relevant to s-process. We have also presented the MACS values over a range of energy from 5 to 100 keV for neutron magic nuclei with (N = 50). PACS numbers: I. INTRODUCTION have been devoted to this respect in order to measure the thermal neutron capture cross sections relevant to s-process temperature. However, reactions on some im- Elements heavier than iron are produced via two prin- portant nuclei are still not available due to their unavail- ciple processes, namely, the slow neutron capture pro- ability in the terrestrial laboratory. K¨appeler et al. [4] cess (s-process) and the rapid neutron capture process showed the present status of the uncertainty of stellar (r-process), differing in the respective neutron capture (n,γ) cross sections and commented that improvements timescales with respect to the β-decay half-lives. There are certainly necessary especially in the mass region be- is a minor contribution from another process, namely, p- low A=120 and above A=180. Some reactions are of sig- process, producing a subset of proton-rich isotopes. The nificant importance, basically those with closed neutron detailed study of the heavy element nucleosynthesis was shells acting as bottlenecks to the s-process reaction flow. done in the fundamental work of Burbidge et al. [1] and They have very low cross-sections and in some scenarios, also of Cameron [2]. Recently K¨appeler et al. [3] pre- the s-process reaction flow cannot overcome these bottle- sented a review of progress on the studies of s-process necks. The neutron capture cross sections of the branch- nucleosynthesis with advanced nuclear physics inputs, point nuclei in the s-process path also have to be known observational data, and stellar models. with better accuracy. The branching in the s-process While the majority of the theory of the s-process is path appears when a competition between neutron cap- well-developed, uncertainty still remains in constraining ture and β-decay takes place due to the appearance of the neutron capture rates of nuclei involved in nucle- a long-lived nucleus with a larger probability of neutron osynthesis chain. The capture cross sections are highly capture than β-decay. Further propagation depends on scattered and uncertain in the energy range appropriate the neutron density and temperature and hence branch- for astrophysical applications. These uncertainties lead arXiv:1604.06757v1 [nucl-th] 22 Apr 2016 ing analysis can be used as an efficient monitor of stel- to significant errors in determining exact abundances of lar neutron density and temperature. These branchings elements involved in different processes. Many works also have an overall influence in changing the nucleosyn- thesis path of s-process and hence on the final yields of the nuclei produced. Thus, more accurate cross sections are nowadays also needed to resolve the discrepancies in ∗Electronic address: [email protected] †Electronic address: [email protected] overproduction problems of some elements. ‡Electronic address: [email protected] The s-process is subdivided into weak, main, and 2 strong components. Recently, another component, Rayet et al. found that the suppression of (n,γ) channel called lighter element primary process (LEPP) is also leads to change the p-nuclei overproduction factors, in proposed. The neutron exposure for weak s-process is some cases by order as large as 10 or 100 [8]. Hence, ac- too low to achieve flow equilibrium and hence a small curate determination of cross sections is highly needed uncertainty in capture cross section may lead to signifi- for them. Particularly the (n,γ) cross sections on p- cant changes in abundance for a number of heavier ele- nuclei are difficult to measure as they are not found in ments. This so-called propagation problem requires new significant amounts for time-of-flight (TOF) or activa- data for the weak component in between Fe to Sr. Our tion measurements. Experimental data are very scarce present study deals with nuclei near N = 50 shell clo- in this respect and by far most of the rates are to be sure where mainly the weak s-process component creates inferred from statistical model calculations. elements between Fe to Sr-Y-Zr region. In a previous In the domain of N = 50 shell proximity, some nuclei work, we studied the thermal (n,γ) cross sections near exist which are predominantly produced via s-process N = 82 shell closure in between isotopes of xenon to only. These so called s-only isotopes (86Sr, 88Sr, and samarium [5] which are produced by the main compo- 96Mo) are of special importance as they can provide nec- nent of s-process nucleosynthesis. The astrophysical ori- essary clues to the s-process branchings. Hence, accu- gins of the weak and main components are completely racy is also required in the cross sections of these s-only different. The weak component occurs during the con- isotopes. vective He core and C-shell burning in massive stars driv- Some other important aspects, for which more accu- ing the material with masses 56 < A < 90, in contrast rate and enhanced nuclear data are required, are those to the main component, occurring mostly in asymptotic of unstable elements having isotopic anomalies. For ex- giant branch (AGB) stars. However, there remains some ample, in the region of N = 50 shell closure, Kr and confusion about the marginal border between these two Zr isotopes show anomalies in meteoritic isotopic dis- components and most of the nuclei considered in the tribution. There remain discrepancies in between mea- present study reside in the vicinity of this border re- sured and model-based isotopic ratios and a spread in gion. They are considered to be produced in both main the ratio with grain size of the material is a persistent and weak components of stellar s-process nucleosynthe- problem.The isotopes 93Zr and 99Tc are of particular sis. For example, 50% of the solar abundances of the interest as they are being used in the nuclear transmuta- 86,87 s-only pair Sr are formed in the main component tion of long-lived fission products (LLFP). The capture and remaining 50% are from the weak component of s- cross sections of the stable zirconium isotopes as well process [4]. In the concerned mass region, there are six as 93Zr are very important for the transmutation study stable nuclei on the nucleosynthesis path with magic neu- on 93Zr. Similar to 93Zr, the isotope 99Tc with half- tron number N = 50, namely, 86Kr, 87Rb, 88Sr, 89Y, and 5 90 92 life of 2.11 10 years and a large cumulative yield in Zr, and Mo. the thermal× neutron fission of 235U, is also one of the The present study also involves some nuclei produced most important LLFPs for transmutation technology to in astrophysical p-process near N = 50 closed shell (84Sr, convert the radioactive wastes into stable or short-lived 92,94Mo, 96,98Ru). The p-nuclei, those are 10 to 100 nuclei. The unstable isotopes are not available for direct times less abundant than s- or r-nuclei in the same mass measurement. Even the experimental data on stable iso- range are produced in the high-temperature environment topes are inadequate and hence, theoretical estimates are via the sequences of photo-disintegration and β+ decays strongly required for them. driving the path through the extreme proton-rich side. It is also possible to predict the values of stellar mean However, current theoretical models are incapable of re- neutron exposure relevant to s-process from the measure- producing the p-nuclei abundances in good agreements ment of (n,γ) cross sections. For example, Bauer et al. (except in certain mass range) with observation and this [9] recommended a set of neutron density, temperature, inadequacy is directly related to the absence of proper and mean neutron exposure from the cross section mea- nuclear inputs with reduced uncertainties. Thus, (n,γ) surements, both for main and weak components. There reactions have significant impact in p-process abundance are a few other works which investigated the same quan- determination as they hinder the flux via competing with tities from radiative cross sections [3]. The (n,γ) cross (γ,n) reactions. Some studies [6, 7] have revealed a high sections can also be used in the study of nuclear cos- degree of sensitivity towards the (n,γ) rate variations in mochronology to determine the age of the Galaxy [10– p-process during the network calculations. For example, 12]. Though the neutron capture cross sections, in gen- 3 eral, have 1/v dependence, this can significantly differ tion. Next, in section III, we have discussed the results when the p-wave capture is superimposed on pure s- obtained from our theoretical approach.