Development of a low-background neutron detector array Enhanced Publication

Development of a low-background neutron detector array 1 2 first-author https://orcid.org/0000-0001-5638-6406 3 corresp linwp1204@scu.edu.cn linwp1204@scu.edu.cn linwp1204@scu.edu.cn https://orcid.org/0000-0002-1971-4025 1 2 4 corresp gaobsh@impcas.ac.cn gaobsh@impcas.ac.cn gaobsh@impcas.ac.cn 1 4 3 1 2 1 2 1 2 4 https://orcid.org/0000-0002-0029-4999 1 2 5 https://orcid.org/0000-0002-4546-2781 6 6 https://orcid.org/0000-0003-2738-537X 1 2 6 1 2 1 2 CAS Key Laboratory of High Precision Nuclear Spectroscopy, Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China CAS Key Laboratory of High Precision Nuclear Spectroscopy, Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China School of Nuclear Science and Technology, University of Chinese Academy of Sciences, Beijing 100049, China School of Nuclear Science and Technology, University of Chinese Academy of Sciences, Beijing 100049, China Key Laboratory of Radiation Physics and Technology of the Ministry of Education, Institute of Nuclear Science and Technology, Sichuan University, Chengdu 610064, China Key Laboratory of Radiation Physics and Technology of the Ministry of Education, Institute of Nuclear Science and Technology, Sichuan University, Chengdu 610064, China Joint Department for Nuclear Physics, Lanzhou University and Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China Joint Department for Nuclear Physics, Lanzhou University and Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China Joint Department for Nuclear Physics, Lanzhou University and Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China Joint Department for Nuclear Physics, Lanzhou University and Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China China Institute of Atomic Energy, Beijing 102413, China China Institute of Atomic Energy, Beijing 102413, China A low-background neutron detector array was developed to measure the cross section of the 13 C( α ,n) 16 O reaction, which is the neutron source for the s -process in AGB stars, in the Gamow window ( E c.m. = 190 ± 40 keV) at the China Jinping Underground Laboratory (CJPL). The detector array consists of 24 3 He proportional counters embedded in a polyethylene cube. Owing to the deep underground location and a borated polyethylene shield around the detector array, a low background of 4.5(2)/h was achieved. The 51 V(p, n) 51 Cr reaction was used to determine the neutron detection efficiency of the array for neutrons with energies E n < 1 MeV. Geant4 simulations are shown to effectively reproduce the experimental results. They were used to extrapolate the detection efficiency to higher energies for neutrons emitted in the 13 C( α ,n) 16 O reaction. The theoretical angular distributions of the 13 C( α ,n) 16 O reaction were shown to be important in the estimation of the uncertainties of the detection efficiency. Underground laboratory Neutron detector Low background 3 He Counter Introduction 1 The 13 C( α ,n) 16 O reaction is the dominant neutron source for the slow neutron capture process (s-process), which synthesizes roughly half of the elements heavier than iron in the Universe [ 1 1 ]. The main site of the s-process is the “He intershell” of the asymptotic giant branch (AGB) stars where the temperature is 0.1 GK. This corresponds to a Gamow window of 190 ± 40 keV for 13 C( α ,n) 16 O reaction, which is far below the Coulomb barrier. Theoretical calculations predicted a cross section as low as ∼ 10 -14 barn at 190 keV [ 2 2 , 3 3 ]. The 13 C( α ,n) 16 O reaction is also the main neutron source for the intermediate neutron capture process (i-process), with a typical temperature of 0.2 GK [ 4 4 , 5 5 ]. A reliable direct measurement is still missing within the corresponding energy range of 0.2–0.54 MeV. To push the direct measurements of the cross section towards the Gamow Window, high-intensity beams and the low-background neutron detector array are indispensable. Present measurements performed at ground laboratories are limited by the cosmic-ray-induced background. With a typical background counting rate of a few hundred per hour (e.g. [ 6 6 , 7 7 ]), those measurements only reached a lower limit of E α ∼ 400 keV with a cross section of ∼ 10 -10 barn. The China JinPing underground Laboratory (CJPL) [ 8 8 ] with approximately 2400 m rock overburden (6700 m water equivalent) shields against most of the cosmic rays. The ultra-low background provides new opportunities for neutrino physics [ 9 9 ], dark matter [ 10 10 , 11 11 ], and direct cross section measurements of the key reactions in nuclear astrophysics [ 12 12 , 13 13 ]. The Jinping Underground Nuclear Astrophysics experimental facility (JUNA) [ 12 12 , 13 13 ] has been initiated to measure important stellar nuclear reaction rates, taking advantage of the low-background environments in CJPL. The present JUNA project includes 12 C( α , γ ) 16 O, 13 C( α ,n) 16 O, 19 F(p, γ ) 20 Ne and 25 Mg(p, γ ) 26 Al reactions [ 12 12 , 13 13 ]. A high-intensity accelerator capable of delivering proton and α beams up to 10 mA in the energy range of 0.05-0.4 MeV has been built for JUNA [ 12 12 ]. This has made it possible to measure the 13 C( α ,n) 16 O cross section at the region of the s-process Gamow window when combined with a high-efficiency and low-background neutron detection system. This paper describes the development and characterization of a low-background neutron detector array for the cross section measurement of the 13 C( α , n) 16 O reaction. Design of the detector array 1 The detector array consists of 24 3 He-filled proportional counters (manufactured by GE-Reuter Stokes), embedded in a 50 × 50 × 50 cm 3 polyethylene cubes. Each proportional counter was filled with 3 He and argon gases at pressures of 4 bar and 2.4 bar, respectively. The active length of each counter was 30 cm, with a diameter of 2.54 cm. A schematic representation of the detector array is shown in Fig. 1 Fig. 1 . The 24 proportional counters are evenly distributed in two concentric rings with radii R 1 = 8.5 cm and R 2 = 13 cm, respectively. The radius of the two rings were optimized using Geant4 simulations, following the same procedure as that in [ 14 14 ], to obtain the maximum detection efficiency for 2.0-MeV isotropic neutrons. A central bore hole with a diameter of 10 cm was created to accommodate the beam line and target. Neutrons from the 13 C( α , n) 16 O reaction with energies of 2–3 MeV are first moderated by the polyethylene cube and then captured by the proportional counters. To shield against environmental neutrons, the polyethylene cube is wrapped with a 5 cm-thick 7% borated polyethylene layer. The signal of each detector was processed using a charge-sensitive preamplifier (model CAEN A1422). The output waveforms were digitized and recorded using a XIA Pixie-16 100 MHz card. Figure 2(a) Figure 2(a) shows a typical energy spectrum measured by the detector array using a 252 Cf source. The “Neutron & alpha” cut shown in Fig. 2(b) Fig. 2(b) was used in the plot to select the neutron events (see the following sections for details). The peak at 0.76 MeV corresponds to full-energy events, where the triton and proton produced from the 3 He + n → p + t reaction deposit all their kinetic energies in the sensitive region of the counter. If the proton or triton hits the wall of the detector before depositing all their kinetic energies in the sensitive region, the detector records only part of the reaction Q value. Those are wall-effect events that correspond to the flat region with lower energies starting from 0.18 MeV in the spectrum. The argon gas filled in the proportional counters has little effect on the neutron detection efficiency but can significantly decrease the track lengths of the produced proton and triton, and therefore reduce the wall effect. By defining a region of interest around the full energy peak between 0.64–0.8 MeV (ROI) as indicated in Fig. 2(a) Fig. 2(a) , we find that 81.0±0.4% of the detected events are located in it. Background measurements 1 Stainless steel is chosen as the wall material of the proportional counters because it has a significantly lower α emission rate compared to aluminum [ 15 15 , 16 16 ]. The inset in Fig. 2(a) Fig. 2(a) shows an energy spectrum up to 7 MeV measured at the CJPL. The events located beyond the neutron full-energy peak are mainly due to the α radioactivity from the wall material of the counters. The α spectrum has a relatively flat distribution below ∼ 3 MeV and exhibits a peak at 5 MeV, a typical energy of α decays from U/Th nuclides. Below the 0.3 MeV energy region, numerous β/γ events exist, overlapping with the neutron events. Following the procedure described in Ref. [ 15 15 ], the β/γ events can be well separated from neutrons in the raise time versus the energy two-dimensional plot as shown in Fig. 2(b) Fig. 2(b) taken at the CJPL. To characterize the background of the detector array, ground- and underground-measurements were performed at the China Institute of Atomic Energy and the A1 hall of the CJPL, respectively. Neutron events were selected using the “neutron & alpha" cut from the rise time versus energy two-dimensional plot as shown in Fig. 2(b) Fig. 2(b) . Such events also include contributions from the α background. To disentangle the α background from the “neutron & alpha" events, a flat α energy spectrum is assumed [ 15 15 , 16 16 ] in the energy range 0.18–2.0 MeV. The total background inside the "neutron and alpha" cut is shown in Fig. 2(b) Fig. 2(b) measured at CJPL for the whole array as 4.5(2)/h, of which 1.94(5)/h is from the α background. The above α background is evaluated using the α counts in the energy range 1.0–2.0 MeV, which was assumed to have the same α emission rate as 0.18–0.8 MeV. This assumption was verified using a 4 He proportional counter, as discussed below. The total background was more than two orders of magnitude lower than the result of 1238(11)/hour obtained from ground-measurement. A 4 He counter, which has exactly the same parameters as the 3 He counters but filled with 4 He gas instead of 3 He, manufactured by the same supplier was used to investigate the α background in more detail. The 4 He counter is insensitive to neutrons but expected to have similar α background shape as the 3 He counters. The α background rate, with a 28.9-hour measurement time in CJPL, was found to be 4.24 ± 0.48 and 4.22 ± 0.38 MeV -1 h -1 , in 0.18–0.8 and 1.0–2.0 MeV energy regions, respectively. This indicates that the two energy regions have the same α emission rates. However, we found an α emission rate of 2.63(30)/hour using the same neutron cut as the 3 He counters in the raise time versus the energy plot. This is approximately a factor of 30 higher than the 3 He counters. It should be noted that the α emission rate of 1.94(5)/h quoted above includes the background of the 24 3 He counters while the 2.63(30)/h is only for one 4 He counter. This implies that the α background could differ significantly even for the same type of wall material. This has to be carefully evaluated at the manufacturing stage to achieve a low α background. To further reduce the α background, one can choose events within the ROI only [see Fig. 2(a) Fig. 2(a) ]. Consequently, the α background was reduced by 74% at the cost of losing only approximately 19% of the neutron detection efficiency. This corresponds to a background of 2.5(1)/h, of which 0.50(1)/h is from the α background. A better approach to the α -background suppression is to use a plastic scintillator as the moderator, where most of the α background can be eliminated via coincident measurements [ 17 17 ]. The α background is the intrinsic properties of the proportional counters used in the array, which does not change with environments. However, the neutron background may change with different environments and therefore should be characterized for each use. Efficiency calibration 1 51 V(p, n) 51 Cr Experiment and simulation 2 The detection efficiency of the array was calibrated at the nuclear physics experiment (NPE) terminal of the 3-MV tandetron accelerator [ 18 18 ] at Sichuan university. The detection efficiency hereafter includes both full-energy and wall-effect events. Quasi-mono-energetic neutrons were produced using the 51 V(p, n) 51 Cr reaction ( Q = -1534.8 keV) at incident energies between 1.7 and 2.6 MeV with a step of 0.15 MeV. The 51 V(p, n) 51 Cr reaction is widely used in the calibration of the neutron detectors [ 14 14 , 16 16 , 19 19 - 22 22 19-22 ] owing to the slow variation in the neutron intensity and energy with the angle, and the well-known target preparation and utilization. Figure 3 Figure 3 shows the emitted neutron energy as a function of the incident proton beam energy. When the incident energy is above 2.33 MeV, neutrons from the transitions feeding the first excited state of 51 Cr ( E x = 749 keV) were mixed with those feeding the ground state. However, as pointed out in Ref. [ 16 16 ], the contribution from the transitions to the first excited state is negligible for incident energies up to 2.6 MeV [ 19 19 ], which was also confirmed in the present work (see discussions in the following sections). The proton beams were focused on a target with a diameter of less than 5 mm. (full width at half maximum). The beam spots were monitored using the fluorescence target with a beam position monitor located 1.5 m upstream the target for each beam energy. Vanadium targets with thicknesses of 110舁micro g/cm 2 were used in the measurements. They were produced by evaporating natural vanadium on 1-mm thick tantalum disks with a diameter of 30 mm. The target thickness corresponds to an energy loss of 8 ∼ 11 keV for proton beams with energies of 1.7–2.6 MeV. The beam intensity varied from 4舁μA to 120 nA as the beam energy increased from 1.7 to 2.6 MeV, making the counting rate of the detector array remain at approximately 10 4 /s with a dead time of less than 1% for all beam energies. A direct water cooling of the reaction target was used to reduce the sputtering and target loss. The possible contribution of the beam induced background was evaluated with a blank tantalum target at the beam energy of 2.0 MeV. The background contribution was found to be less than 1%. The total number of emitted neutrons was determined based on the activation method as described in Ref. [ 16 16 ]. The number of radioactive products 51 Cr equals to the number of emitted neutrons, during the reaction. 51 Cr decays via the electron capture with a half-life of T 1/2 = 27.7025(24) days, and has a branching ratio of B = 9.91(1)% to decay to the first excited state of its daughter nuclei 51 V, which is followed by the emission of a 320-keV γ ray. Offline measurements of the γ ray were conducted using a GEM-series HPGe detector whose relative efficiency was 30%. The absolute efficiency of the HPGe detector was measured at a distance of 20 cm using 137 Cs [1.534(19)×10 5 Bq] and 152 Eu [5.72(6)×10 4 Bq] γ -ray sources. The distance of 20 cm is large enough to avoid pileups of the cascading γ rays from the sources. A 51 Cr γ source was produced by a 51 V(p, n) 51 Cr reaction. Its activity was measured at a position of 20 cm and then used to determine the efficiency of the 320-keV γ line at a position of 10.2 cm ( η 320 ), where all the irradiated targets were placed for off-line measurements. η 320 was determined to be 0.498(7)% at 10.2 cm. Because the beam current was stable within the irradiation time, the number of emitted neutrons was determined by the offline measurement of the 320 keV γ ray with the activation formula [ 16 16 ] <math display="block" id="M1"><mrow><msub><mi>N</mi><mtext>R</mtext></msub><mo>=</mo><mfrac><mrow><msub><mi>N</mi><mi>γ</mi></msub></mrow><mrow><mi>B</mi><mo>⋅</mo><msub><mi>η</mi><mrow><mn>320</mn></mrow></msub></mrow></mfrac><mo>⋅</mo><mfrac><mrow><msup><mi>e</mi><mrow><mi>λ</mi><msub><mi>t</mi><mtext>w</mtext></msub></mrow></msup></mrow><mrow><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><msub><mi>t</mi><mtext>c</mtext></msub></mrow></msup></mrow></mfrac><mo>⋅</mo><mfrac><mrow><mi>λ</mi><mo>⋅</mo><msub><mi>t</mi><mtext>i</mtext></msub></mrow><mrow><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><msub><mi>t</mi><mtext>i</mtext></msub></mrow></msup></mrow></mfrac><mn>,</mn></mrow></math> where N γ , t i , t c , and t w are the number of the detected 320-keV γ rays, the activation time, the counting time, and the waiting time elapsed between the end of irradiation and the start of the counting, respectively. λ = (ln2)/ T 1/2 is the decay constant of 51 Cr. After the dead time correction, the detection efficiency of the neutron detector is calculated as <math display="block" id="M2"><mrow><msub><mi>η</mi><mtext>n</mtext></msub><mo>=</mo><mfrac><mrow><msub><mi>N</mi><mtext>n</mtext></msub></mrow><mrow><msub><mi>N</mi><mtext>R</mtext></msub></mrow></mfrac><mn>,</mn></mrow></math> where N n is the detected number of neutrons by the detector array. Because the energy of the neutrons emitted in the 51 V(p, n) 51 Cr reaction is only up to ∼ 1 MeV, it is still far below that of the 13 C( α , n) 16 O reaction. The Monte Carlo simulation code, Geant4 [ 23 23 , 24 24 ] of version 10.4.6 was used to determine the neutron detection efficiency at neutron energies above 1 MeV. The detailed detector setup used in the simulation is shown in Fig. 4 Fig. 4 . The beam pipe, target backing, and water cooling loops were included in the simulation, which resembles the physical setups. An overestimation of the neutron detection efficiency was obtained from the simulation compared to the experimental results. Similar overestimations have also been observed in several other setups [ 14 14 , 16 16 ]. This overestimation may be because some of the neutrons are absorbed by small contaminants in the moderating polyethylene. Instead of using a normalization factor, we added a small amount of boron into the moderating polyethylene in the simulation to consider the neutron absorption effects. Figure 5(a) Figure 5(a) shows the reduced χ 2 , obtained by comparing the simulated and measured detection efficiencies, as a function of the boron mass fraction. The minimum of χ 2 was found at the boron mass fraction of 0.054% by fitting the curve using a parabola. The simulated total, inner ring, and outer ring detection efficiencies are shown in Fig. 5(b) Fig. 5(b) . An estimation of the first excited state contributions was also carried out using the first excited state neutron component estimated from a statistical model calculation, calibrated using experimental data [ 19 19 ] and TALYS [ 25 25 ]. The simulated total efficiencies with the first excited state neutron contributions are also shown by the blue short-dashed curve in Fig. 5(b) Fig. 5(b) . The maximum deviation between the simulation with and without the contribution of the first excited state was 0.9%, which is negligible compared to the experimental uncertainty of 3.7%. The boron contaminant affects not only the total efficiency but also the ratio R in/out of the detection efficiencies of the inner and outer ring detectors. With the boron mass fraction obtained above, the measured ratio R in/out was also well reproduced using the two simulations (see Fig. 5(c) Fig. 5(c) ). Except for the first point, the simulation results are 1%-3% lower than those of the experiment at E p < 2.4 MeV. Isotropic angular distributions were used in the simulations. The discrepancy in R in/out between the simulations and experiment is probably due to the non-isotropic distribution of the emitted neutrons. For E p = 2.45 and 2.6 MeV, the discrepancy between the simulated (with ground state neutrons) and the measured R in/out ratios increase to 1.6% and 5%, respectively. This is due to the opening of the decay channel to the first excited state in 51 Cr (see Fig. 3 Fig. 3 ). Considering the first excited state neutrons estimated from TALYS, the discrepancy is reduced to 0.6% and 2.3% for E p = 2.45 and 2.6 MeV, respectively, as shown in Fig. 5(c) Fig. 5(c) . The contribution from the decay channel to the first excited state is negligible regarding the total detection efficiency although its effects on the R in/out is larger. The total, inner ring and outer ring detection efficiencies for monoenergetic neutrons up to 4.5 MeV with an isotropic angular distribution were also simulated, as shown in Fig. 5(d) Fig. 5(d) . The dependence of the detection efficiency on the source position was measured by placing the detector array at different positions along the beam line ( Z axis) using the 51 V(p, n) 51 Cr reaction at E p = 2 MeV. Position Z = 0 corresponds to the target being at the center of the array. Moving the detector array forward in the beam direction corresponds to positive Z values and backward corresponds to negative Z values. The relative efficiencies at different positions were normalized using integrated incident beam currents on the target. The results for the total, inner-ring, and outer-ring detectors are shown in Fig. 6 Fig. 6 , together with the Geant4 simulations. In Fig. 6 Fig. 6 , the measured total efficiency at Z = 0 was normalized to the simulated value. Overall good agreements were found between the measurements and simulations. Considering the average deviation between the experimental data and simulations in Fig. 5(b) Fig. 5(b) and Fig. 6 Fig. 6 , the difference between the Geant4 simulations and experimental results is 2.8%, which reflects one of the systematic uncertainties of our simulation and is quoted as the uncertainty of the Geant4 simulations. Extrapolating the detection efficiency for the study of the 13 C( α , n) 16 O reaction at stellar energies 2 It should be noted that the simulated efficiency in Fig. 5 Fig. 5 cannot be directly applied to the 13 C( α , n) 16 O reaction in which the emitted neutrons are neither mono-energetic nor isotropic. The asymmetry in the efficiency curves, shown in Fig. 6 Fig. 6 indicates that the simulation predicts a slightly higher detection efficiency for the neutrons emitted at backward angles than that at forward angles. Because the angular distribution of the 13 C( α , n) 16 O reaction is not measured at low energies close to the Gamow window, theoretically predicted angular distributions were used in our simulations. Legendre polynomials up to the third order were used in calculating the angular distributions. Figure 7 Figure 7 shows the Legendre polynomial coefficient a 1 used in the angular distributions of the 13 C( α , n) 16 O reaction as a function of the beam energies in (a) and the representative angular distributions obtained at beam energies of 0.2, 0.7, 1.0534, and 1.2 MeV in (b). Adopting the angular distributions in Refs. [ 26 26 , 27 27 ], the angular distribution effect on the efficiencies is corrected in our Geant4 simulations. The detection efficiencies were simulated using both isotropic and calculated angular distributions in the incident α energies between 0.2–3.15 MeV. In the simulation, the narrow resonances are ignored. Figure 8(a) Figure 8(a) shows the simulated efficiencies using isotropic and calculated angular distributions. For α energies between 0.3-0.8 MeV, the difference is 0.19%–2.28%, which is relatively small. However, above 0.8 MeV, the difference increases with a maximum of 4%. To evaluate the effect of the sharp resonance, we also simulated the detection efficiency at the energy where the a 1 coefficient reached its maximum at approximately E = 1.0534 MeV, as shown in Fig. 7(a) Fig. 7(a) . The resulting efficiency was 11% less than that with an isotropic angular distribution. Therefore, the angular distribution effect needs to be carefully evaluated for experiments aimed at achieving high precision. The effects of the angular distributions were further evaluated by comparing the simulated and measured ratios of the inner and outer ring detection efficiencies R in/out . As shown in Fig. 8(b) Fig. 8(b) , the R in/out is insensitive to the angular distribution and the simulations in both cases agree reasonably with the measured values from an underground measurement. The difference between the measured and simulated R in/out values was ∼ 4%. An alternative method for extrapolating the total detection efficiency was carried out as follows. The inner-ring detection efficiency was extrapolated to the neutron energy range relevant to the 13 C( α , n) 16 O reaction. The total detection efficiency was then obtained using the measured R in/out values from an underground measurement. The difference in the total detection efficiency between the two extrapolation methods is ∼ 2.5%, which is considered as another systematic uncertainty. Considering the uncertainties of the Geant4 simulations (2.8%), extrapolation (2.5%), neutron angular distribution (2.3%), and detection efficiency of the HPGe detector (1.5 %), the overall uncertainty of the detection efficiency for the 13 C( α , n) 16 O reaction in the energy range E α = 300–800 keV is determined as 5%. Excluding the narrow resonances, the maximum uncertainty of the neutron angular distributions is 4% if the energy range is extended to 2.4 MeV. According to the ENDF angular distribution, the maximum uncertainty of the neutron angular distributions at energies greater than 2.4 MeV is 12%. Summary 1 A high-efficiency and low-background neutron detector array consisting of 24 3 He proportional counters embedded in a polyethylene moderator was developed for the cross-section measurement of the 13 C( α , n) 16 O reaction at the China Jinping underground laboratory. Because of the deep underground location and 5-cm thick borated polyethylene shield, a low background of 4.5(2)/hour was achieved, of which 1.94(5)/hour was from the internal α radioactivity. The 51 V(p, n) 51 Cr reaction was used to calibrate the neutron detection efficiency of the array for neutrons with energies of E n < 1 MeV. For E n > 1 MeV, Monte Carlo simulation code Geant4 was used to extrapolate the detection efficiency. Specifically, an energy dependent detection efficiency, which can be directly applied to the 13 C( α , n) 16 O reaction, was obtained from the simulation. The effects of the angular distribution of the 13 C( α ,n) 16 O reaction on the detection efficiency, which were overlooked in previous works, were evaluated in the present work and were shown to be non-negligible. References 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 . Astrophysical J. 497 , 388 ( 1998 ). doi: 10.1086/305437 http://doi.org/10.1086/305437 Evolution and nucleosynthesis in low-mass asymptotic giant branch stars. II. neutron capture and the s-process B. Guo , Z. H. Li , M. Lugaro et al. , New determination of the 13 C( α ,n) 16 O reaction rate and its influence on the s-process nucleosynthesis in AGB stars . Astrophysical J. 756 : 193 - 195 ( 2012 ). doi: 10.1088/0004-637x/756/2/193 http://doi.org/10.1088/0004-637x/756/2/193 New determination of the 13C(α,n)16O reaction rate and its influence on the s-process nucleosynthesis in AGB stars R.J. deBoer , C.R. Brune , M. Febrarro et al. , Sensitivity of the 13 C( α,n ) 16 O S factor to the uncertainty in the level parameters of the near-threshold state . Phys. Rev. C 101 , 045802 ( 2020 ). doi: 10.1103/PhysRevC.101.045802 http://doi.org/10.1103/PhysRevC.101.045802 Sensitivity of the 13C(α,n)16O S factor to the uncertainty in the level parameters of the near-threshold state F. Herwig , M. Pignatari , P.R. Woodward et al. , Convective-reactive proton-12C combustion in sakurai’s object(V4334 SAGITTARII) and implications for the evolution and yields from the first generations of stars . Astrophysical J. 727 , 89 ( 2011 ). doi: 10.1088/0004-637X/727/2/89 http://doi.org/10.1088/0004-637X/727/2/89 Convective-reactive proton-12C combustion in sakurai’s object(V4334 SAGITTARII) and implications for the evolution and yields from the first generations of stars P. Banerjee , Y.-Z. Qian , A. Heger , New neutron-capture site in massive Pop III and Pop II stars as a source for heavy elements in the early galaxy . Astrophysical J. 865 120 ( 2018 ). doi: 10.3847/1538-4357/aadb8c http://doi.org/10.3847/1538-4357/aadb8c New neutron-capture site in massive Pop III and Pop II stars as a source for heavy elements in the early galaxy H.W. Drotleff , A. Denker , H. Knee et al. , Reaction rates of the s-process neutron sources 22 Ne (alpha, n) 25 Mg and 13 C( alpha, n) 16 O . Astrophysical J. 414 735 ( 1993 ). doi: 10.1086/173119 http://doi.org/10.1086/173119 Reaction rates of the s-process neutron sources 22Ne (alpha, n) 25Mg and 13C( alpha, n) 16O S. Harissopulos , H.W. Becker , J.W. Hammer et al. , Cross section of the 13 C( α,n ) 16 O reaction: A background for the measurement of geo-neutrinos . Phys. Rev. C 72 062801 ( 2005 ). doi: 10.1103/PhysRevC.72.062801 http://doi.org/10.1103/PhysRevC.72.062801 Cross section of the 13C(α,n)16O reaction: A background for the measurement of geo-neutrinos J.P. Cheng , K.J. Kang , J.M. Li et al. , The China Jinping Underground Laboratory and its early science . Ann. Rev. Nucl. Part. Sci. 67 231 - 251 ( 2017 ). doi: 10.1146/annurev-nucl-102115-044842 http://doi.org/10.1146/annurev-nucl-102115-044842 The China Jinping Underground Laboratory and its early science Z.Y. Li , Y.M. Zhang , G.F. Cao et al. , Event vertex and time reconstruction in large-volume liquid scintillator detectors . Nucl. Sci. Tech. 32 , 49 ( 2021 ). doi: 10.1007/s41365-021-00885-z http://doi.org/10.1007/s41365-021-00885-z Event vertex and time reconstruction in large-volume liquid scintillator detectors K.L. Giboni , P. Juyal , E. Aprile et al. , A LN 2 -based cooling system for a next-generation liquid xenon dark matter detector . Nucl. Sci. Tech. 31 76 ( 2020 ). doi: 10.1007/s41365-020-00786-7 http://doi.org/10.1007/s41365-020-00786-7 A LN2-based cooling system for a next-generation liquid xenon dark matter detector P. Juyal , K.L. Giboni , X.D. Ji et al. , On proportional scintillation in very large liquid xenon detectors . Nucl. Sci. Tech. 31 , 93 ( 2020 ). doi: 10.1007/s41365-020-00797-4 http://doi.org/10.1007/s41365-020-00797-4 On proportional scintillation in very large liquid xenon detectors W.P. Liu , Z.H. Li , J.J. He et al. , Progress of Jinping Underground laboratory for Nuclear Astrophysics (JUNA) . Science China Physics, Mechanics & Astronomy , 59 , 4 ( 2016 ). doi: 10.1007/s11433-016-5785-9 http://doi.org/10.1007/s11433-016-5785-9 Progress of Jinping Underground laboratory for Nuclear Astrophysics (JUNA) L.Y. Zhang , J. Su , J.J. He et al. , Direct measurement of the astrophysical 19 F( p,αγ ) 16 O reaction in the Deepest Operational Underground Laboratory . Phys. Rev. Lett. 127 , 152702 ( 2021 ). doi: 10.1103/PhysRevLett.127.152702 http://doi.org/10.1103/PhysRevLett.127.152702 Direct measurement of the astrophysical 19F(p,αγ)16O reaction in the Deepest Operational Underground Laboratory S. Falahat , A. Best , M. Couder et al. , A 3 He neutron detector for the measurement of ( α ,n) reactions . Nucl. Instrum. Meth. A , 700 , 53 - 58 ( 2013 ). doi: 10.1016/j.nima.2012.10.036 http://doi.org/10.1016/j.nima.2012.10.036 A 3He neutron detector for the measurement of (α,n) reactions T.J. Langford , C.D. Bass , E.J. Beise et al. , Event identification in 3He proportional counters using risetime discrimination . Nucl. Instrum. Meth. A 717 , 51 - 57 ( 2013 ). doi: 10.1016/j.nima.2013.03.062 http://doi.org/10.1016/j.nima.2013.03.062 Event identification in 3He proportional counters using risetime discrimination L. Csedreki , G.F. Ciani , J. Balibrea-Correa et al. , Characterization of the LUNA neutron detector array for the measurement of the 13C( α , n)16O reaction . Nucl. Instrum. Meth. A 994 , 165081 ( 2021 ). doi: 10.1016/j.nima.2021.165081 http://doi.org/10.1016/j.nima.2021.165081 Characterization of the LUNA neutron detector array for the measurement of the 13C(α, n)16O reaction H.X. Huang , B. Gao , Y.T. Li et al. , Low event rate neutron detector array using the coincidence between plastic scintillator and Helium-3 proportional counters . Nucl. Instrum. Meth. A 1003 , 165323 ( 2021 ). doi: 10.1016/j.nima.2021.165323 http://doi.org/10.1016/j.nima.2021.165323 Low event rate neutron detector array using the coincidence between plastic scintillator and Helium-3 proportional counters J.F. Han , Z. An , G.Q. Zheng et al. , An ion beam facility based on a 3 MV tandetron accelerator in Sichuan University, China . Nucl. Instrum. Meth. B 418 , 68 - 73 ( 2021 ). doi: 10.1016/j.nimb.2018.01.002 http://doi.org/10.1016/j.nimb.2018.01.002 An ion beam facility based on a 3 MV tandetron accelerator in Sichuan University, China G. Deconninck , J. Royen , La reaction 51 V(p,n) 51 Cr comme source de neutrons monoenergetiques . Nucl. Instrum. Meth. 75 266 - 270 ( 1969 ). doi: 10.1016/0029-554X(69)90607-7 http://doi.org/10.1016/0029-554X(69)90607-7 La reaction 51V(p,n) 51Cr comme source de neutrons monoenergetiques J. Pereira , P. Hosmer , G. Lorusso et al. , The neutron long counter NERO for studies of β -delayed neutron emission in the r-process . Nucl. Instrum. Meth. A 618 , 275 - 283 ( 2010 ). doi: 10.1016/j.nima.2010.02.262 http://doi.org/10.1016/j.nima.2010.02.262 The neutron long counter NERO for studies of β-delayed neutron emission in the r-process E. Ramström , T. Wiedling , The excitation function of the 13C( α , n)16O reaction and its astrophysical application . Nucl. Phys. A 272 259 - 268 ( 1976 ). doi: 10.1016/0375-9474(76)90330-4 http://doi.org/10.1016/0375-9474(76)90330-4 The excitation function of the 13C(α, n)16O reaction and its astrophysical application E. Lund , P. Hoff , K. Aleklett et al. , Delayed neutron emission probabilities of gallium, bromine, rubidium, indium, antimony, iodine, and cesium precursors . Zeitschrift für Physik A Atoms and Nuclei , 294 ( 1980 ). doi: 10.1007/BF01438160 http://doi.org/10.1007/BF01438160 Delayed neutron emission probabilities of gallium, bromine, rubidium, indium, antimony, iodine, and cesium precursors S. Agostinelli , J. Allison , K. Amako et al. , Geant4—a simulation toolkit . Nucl. Instrum. Meth. A 506 , 250 - 303 ( 2003 ). doi: 10.1016/S0168-9002(03)01368-8 http://doi.org/10.1016/S0168-9002(03)01368-8 Geant4—a simulation toolkit J. Allison , K. Amako , J. Apostolakis et al. , Recent developments in Geant4 . Nucl. Instrum. Meth. A 835 , 186 - 225 ( 2016 ). doi: 10.1016/j.nima.2016.06.125 http://doi.org/10.1016/j.nima.2016.06.125 Recent developments in Geant4 A.J. Koning and D. Rochman , Modern nuclear data evaluation with the TALYS code system . Nuclear Data Sheets 113 2841 - 2934 ( 2012 ). doi: 10.1016/j.nds.2012.11.002 http://doi.org/10.1016/j.nds.2012.11.002 Modern nuclear data evaluation with the TALYS code system D.A. Brown , M.B. Chadwick , R. Capote et al. , ENDF/B-VIII.0: The 8th major release of the nuclear reaction data library with CIELO-project cross sections, new standards and thermal scattering data . Nuclear Data Sheets 148 1 - 142 ( 2018 ). doi: 10.1016/j.nds.2018.02.001 http://doi.org/10.1016/j.nds.2018.02.001 ENDF/B-VIII.0: The 8th major release of the nuclear reaction data library with CIELO-project cross sections, new standards and thermal scattering data R.B. Walton , J.D. Clement , F. Boreli , Interaction of neutrons with oxygen and a study of the C 13 ( α, n )O 16 reaction . Phys. Rev. 107 , 1065 - 1075 ( 1957 ). doi: 10.1103/PhysRev.107.1065 http://doi.org/10.1103/PhysRev.107.1065 Interaction of neutrons with oxygen and a study of the C13(α, n)O16 reaction National Nuclear Data Center . https://www.nndc.bnl.gov/nudat3/getdatasetClassic.jsp?nucleus=17O&unc=nds https://www.nndc.bnl.gov/nudat3/getdatasetClassic.jsp?nucleus=17O&unc=nds All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Yu-Tian Li, Tao-Yu Jiao, Bing-Shui Gao and Wei-Ping Lin. The first draft of the manuscript was written by Yu-Tian Li, Bing-Shui Gao, Wei-Ping Lin, and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. The authors thank the staff of the 3 MV tandetron accelerator facility of Sichuan University for their support during the experiment. The authors also thank Carl Brune for estimating the neutron angular distribution corrections. 10.1007/s41365-022-01030-0 This work was supported by the National Natural Science Foundation of China (Nos. 11490564 and 11805138). 2021-11-22 2022-02-27 2022-03-07 2022-04-15 2022-11-30 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2022 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-022-01030-0 http://dx.doi.org/10.1007/s41365-022-01030-0 . 2022 Nuclear Science and Techniques 04 33