Yield ratio of neutrons to protons in ¹²C(d,n)¹³N and ¹²C(d,p)¹³C from 0.6 MeV to 3 MeV

Yield ratio of neutrons to protons in 12 C(d,n) 13 N and 12 C(d,p) 13 C from 0.6 MeV to 3 MeV 1 2 first-author 3 1 corresp ygma@sinap.ac.cn ygma@sinap.ac.cn ygma@sinap.ac.cn 4 1 corresp zhangguoqiang@zjlab.org.cn zhangguoqiang@zjlab.org.cn zhangguoqiang@zjlab.org.cn 3 1 5 6 7 3 1 1 1 1 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China University of the Chinese Academy of Sciences, Beijing 100080, China University of the Chinese Academy of Sciences, Beijing 100080, China Key Laboratory of Nuclear Physics and Ion-Beam Application (MOE), Institute of Modern Physics, Fudan University, Shanghai 200433, China Key Laboratory of Nuclear Physics and Ion-Beam Application (MOE), Institute of Modern Physics, Fudan University, Shanghai 200433, China Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201210, China Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201210, China College of Physics and Electronics information, Inner Mongolia University for Nationalities, Tongliao 028000, China College of Physics and Electronics information, Inner Mongolia University for Nationalities, Tongliao 028000, China Cyclotron Institute, Texas A&M University, College Station, Texas 77843, USA Cyclotron Institute, Texas A&M University, College Station, Texas 77843, USA Laboratori Nazionali del Sud, INFN, via Santa Sofia, 62, 95123 Catania, Italy Laboratori Nazionali del Sud, INFN, via Santa Sofia, 62, 95123 Catania, Italy The neutron yield in the 12 C(d,n) 13 N reaction and the proton yield in the 12 C(d,p) 13 C reaction have been measured using deuteron beams of energies 0.6–-3 MeV. The deuteron beam is delivered from a 4 MeV electrostatic accelerator and bombarded on a thick carbon target. The neutrons are detected at 0°, 24°, and 48° and the protons at 135° in the lab frame. Further, the ratio of the neutron yield to the proton yield was calculated. This can be used to effectively recognize the resonances. The resonances are found at 1.4 MeV, 1.7 MeV, and 2.5 MeV in the 12 C(d,p) 13 C reaction, and at 1.6 MeV and 2.7 MeV in the 12 C(d,n) 13 N reaction. The proposed method provides a way to reduce systematic uncertainty and helps confirm more resonances in compound nuclei. Proton neutron ratio 12 C(d,n) 12 C(d,p) Trojan Horse Method (THM) 1 Introduction 1 The reaction 12 C(p, γ ) 13 N has been gaining considerable attention of researchers for both its usefulness as an important reaction to study the properties of light nuclei and the significant role it is expected to play in the cores of main sequence stars and the shells of giant stars [ 1 1 ]. Studying the 12 C(p, γ ) 13 N reaction is important to understand the well-known carbon-nitrogen-oxygen cycle (CNO cycle). Reactions occur in the CNO cycle [ 2 2 ] in the following chain: 12 C(p, γ ) 13 N( β + ) 13 C(p, γ ) 14 N (p, γ ) 15 O ( β + ) 15 N(p, α ) 12 C. The result of the 12 C(d,n) 13 N reaction can be used to study the 12 C(p, γ ) 13 N reaction using the Trojan Horse Method (THM) [ 3 3 ]. Compared with 12 C(p, γ ) 13 N, the cross section of 12 C(d,n) 13 N is three times higher with incident energy in the order of MeV, making it a better alternative to determine the cross section for 12 C(p, γ ) 13 N. The 12 C(d,n) 13 N reaction can be employed to analyze various materials; in addition, also it can be used as a source of polarized neutrons [ 4 4 ]. The 12 C(d,p) 13 C reaction can be employed in metallurgy and Ion Beam Analysis (IBA) [ 5 5 ]. For example, the properties of steel depend on its carbon content and spatial distribution. The cross section and the resonances of nuclear reactions are highly beneficial in understanding the structure of materials at micro level [ 6 6 , 7 7 ]. Moreover, deuterated polyethylene (CD 2 ) n is typically selected as the target for nuclear reactions induced by intensive lasers [ 8 8 , 9 9 ]. At very high plasma temperatures, in addition to the d(d,n) 3 He reaction, the 12 C(d,n) 13 N reaction could occur. Both these reactions emit neutrons, which may result in a fake enhancement of d(d,n) 3 He, thereby complicating the measurement. The cross section measurement for 12 C(d,n) 13 N will help clarify the real cross section of d(d,n) 3 He [ 10 10 - 12 12 10-12 ] in laser induced plasma. The cross sections of the 12 C(d,n) 13 N reaction and the 12 C(d,p) 13 C reaction have been measured in various aspects. The cross section of the 12 C(d,n) 13 N reaction was measured by the neutron detection method using a 4 π detector [ 4 4 ]. The cross section of the 12 C(d,n) 13 N reaction was measured by counting the γ -rays from 13 N produced in the reaction [ 14 14 ]. Further, the cross section of the 12 C(d,p) 13 C reaction was measured by counting the γ -rays [ 5 5 ]. Although much effort has been made, a large diversity of resonances and cross sections exists at approximately 1.4 MeV among the data. In this study, we conduct an experiment to study the ratio between the neutron yield in 12 C(d,n) 13 N and the proton yield in 12 C(d,p) 13 C. This ratio can not only provide information on resonance in the two reactions, but also eliminate the systematic errors caused by the intensity of uncertain incident beam, and it can be applied in thick target reaction systems. 2 Experiment 1 The experiment was conducted at the Shanghai Institute of Applied Physics, Chinese Academy of Sciences. In this experiment, a deuteron beam of energy 0.6–-3 MeV was used to bombard a thick graphite target. Further, the following two reactions were observed in the experiment: 12 C(d,n) 13 N and 12 C(d,p) 13 C. The ratio of the yields of neutrons to protons is extracted. 2.1 Layout of the experiment 2 Figure 1 Figure 1 shows the schematic of the devices used in the experiment. The target chamber was a cylinder of diameter 60.0 cm. The beam tube transports the deuterium beam from the accelerator to the chamber. Further, three neutron detectors were placed outside the chamber, and a carbon target, Faraday’s cup, and proton detector were placed inside the chamber. The deuteron beam used in the experiment was produced from a 4 MeV electrostatic accelerator. The maximum voltage of the accelerator was set as 3 MeV in this experiment, and the fluctuation of the voltage was less than 2 keV. The beam intensity was about several nA, and three ion beam analysis terminals were present. The experiment was conducted at the 15° terminal; this terminal provides conventional ion beam analysis. Furthermore, a vacuum of 6.2 × 10 -5 Pa was established. The graphite target was installed in the center of the chamber, directly facing the deuteron beam. The thickness of the target was 1 mm, which was much larger than the deuteron stopping range of about several μm, in carbon. The 12 C constitutes 98.93% of the natural graphite used in the experiment. The Faraday-cup, which was installed touching the target, was used to count the incident charged particles. However, because the conductivity of the thick graphite target was poor, sparking was unavoidable during the experiment, it was difficult to estimate the intensity of incident beam. Our method of calculating and analyzing the ratio of neutron yield to proton yield provided a means to study the resonances in the two reactions, even without accurate information about beam intensity. 2.2 Detectors 2 In the experiment, four detectors were used to count the yields: three EJ-301 detectors for neutrons and one Au-Si surface barrier detector for protons. The EJ-301 detector was a liquid scintillator detector that is suitable for MeV neutron detection with a good neutron-gamma discrimination. Three EJ-301 detectors were placed at a distance of 112.9 cm from the target at 0°, 24°, and 48° from the beam axis in the experiment. Both the diameter and the thickness of the EJ-301 detectors were 12.7 cm. Further, the solid angles of the three detectors were 1.210 -2 sr. A time-of-flight (TOF) system based on a Cf-252 neutron source was installed to calibrate the EJ-301 detectors. The efficiency of the neutron detector can be calculated by the ratio of the counts of the neutron detector per second and the theoretical radiation rate. The energies and the emission probabilities of the neutrons from the Cf-252 source follow the Maxwell-Boltzmann distribution. The threshold is applied by the ADC of the neutron signal, i.e., the height of the signal pulse, which is at 300 channel to filter background noise. The efficiency of the EJ-301 detector at 24° varied from different incident energies shown in Fig. 2 Fig. 2 ; the final neutron yields were then corrected by this efficiency. The Au-Si surface barrier detector was placed at 135 mm from the target at 135° from the beam axis. The diameter of the Au-Si surface barrier detector was 7.5 mm. The solid angle of the Au-Si surface barrier detector was 2.42 × 10 -3 sr. The efficiency of the Au-Si surface barrier detector was practically 100% for the MeV proton in this experiment. 2.3 Methods 2 <math display="block" id="M1"><mrow><mi>Y</mi><mo>=</mo><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>T</mi></msubsup></mrow></mstyle><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>D</mi></msubsup></mrow></mstyle><mi>I</mi><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>t</mi><mo stretchy="false">)</mo><msub><mi>N</mi><mi>v</mi></msub><mi>σ</mi><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>x</mi><mtext>d</mtext><mi>t</mi><mn>,</mn></mrow></math> where D is the thickness, I ( E , t ) is the incident current, and T is the time duration of the projectile in the target. Further, N v is the atom number per unit volume of the carbon target, And σ ( E ) is the cross section of the reaction at energy E . The ratio of the neutron yield to the proton one in the two reactions is expressed as follows: <math display="block" id="M2"><mrow><mfrac><mrow><msub><mi>Y</mi><mi>n</mi></msub></mrow><mrow><msub><mi>Y</mi><mi>p</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>t</mi><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>D</mi></msubsup></mrow></mstyle><mi>I</mi><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><msub><mi>N</mi><mi>v</mi></msub><msub><mi>σ</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>x</mi></mrow><mrow><mi>t</mi><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>D</mi></msubsup></mrow></mstyle><mi>I</mi><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><msub><mi>N</mi><mi>v</mi></msub><msub><mi>σ</mi><mi>p</mi></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>x</mi></mrow></mfrac><mn>.</mn></mrow></math> Theoretically, a thick target is equivalent to the sum of several continuous thin targets. For example, the yield with an incident energy of 3 MeV can be described as follows. <math display="block" id="M3"><mrow><msub><mi>Y</mi><mrow><mtext>total</mtext></mrow></msub><mo>=</mo><mi>Δ</mi><msub><mi>Y</mi><mrow><mn>0.6</mn></mrow></msub><mo>+</mo><mi>Δ</mi><msub><mi>Y</mi><mrow><mn>0.7</mn></mrow></msub><mo>+</mo><mo>⋯</mo><mo>+</mo><mi>Δ</mi><msub><mi>Y</mi><mrow><mn>3.0</mn></mrow></msub></mrow></math> The thickness of each thin target is the range that the deuteron passes through in the target from energy E to energy E -0.1 MeV. △ Y E represents the yield of the thin target. The total yield is the sum of yields of all the thin target slices. For the two reactions in our work, when the incident energy is below 0.6 MeV, the cross section is less than 10 mb. Therefore, the slices in which the deuteron energy is less than 0.6 MeV will not be considered. For a very thin target, the formula for yield and cross section of the nuclear reaction is given by, <math display="block" id="M4"><mrow><mi>Δ</mi><mi>Y</mi><mo>=</mo><mi>I</mi><mi>Δ</mi><mi>t</mi><msub><mi>N</mi><mi>v</mi></msub><mi>σ</mi><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mi>Δ</mi><mi>x</mi><mn>.</mn></mrow></math> The yield of a thin target is proportional to its cross section. The difference between the ratio of neutron yield to proton yield in the two reactions is given by, <math display="block" id="M5"><mrow><mrow><mo>{</mo><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>Y</mi><mrow><msub><mi>n</mi><mrow><mi>E</mi><mo>+</mo><mn>0.1</mn></mrow></msub></mrow></msub><mo>=</mo><msub><mi>Y</mi><mrow><msub><mi>n</mi><mi>E</mi></msub></mrow></msub><mo>+</mo><mi>Δ</mi><msub><mi>Y</mi><mrow><msub><mi>n</mi><mrow><mi>E</mi><mo>+</mo><mn>0.1</mn></mrow></msub></mrow></msub></mrow></mtd><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mfrac><mrow><msub><mi>Y</mi><mrow><msub><mi>n</mi><mrow><mi>E</mi><mo>+</mo><mn>0.1</mn></mrow></msub></mrow></msub></mrow><mrow><msub><mi>Y</mi><mrow><msub><mi>p</mi><mrow><mi>E</mi><mo>+</mo><mn>0.1</mn></mrow></msub></mrow></msub></mrow></mfrac><mo>−</mo><mfrac><mrow><msub><mi>Y</mi><mrow><msub><mi>n</mi><mi>E</mi></msub></mrow></msub></mrow><mrow><msub><mi>Y</mi><mrow><msub><mi>p</mi><mi>E</mi></msub></mrow></msub></mrow></mfrac><mo>≈</mo><mfrac><mrow><mfrac><mrow><mi>Δ</mi><msub><mi>Y</mi><mrow><msub><mi>n</mi><mrow><mi>E</mi><mo>+</mo><mn>0.1</mn></mrow></msub></mrow></msub></mrow><mrow><msub><mi>Y</mi><mrow><msub><mi>n</mi><mi>E</mi></msub></mrow></msub></mrow></mfrac><mo>−</mo><mfrac><mrow><mi>Δ</mi><msub><mi>Y</mi><mrow><msub><mi>p</mi><mrow><mi>E</mi><mo>+</mo><mn>0.1</mn></mrow></msub></mrow></msub></mrow><mrow><msub><mi>Y</mi><mrow><msub><mi>p</mi><mi>E</mi></msub></mrow></msub></mrow></mfrac></mrow><mrow><mfrac><mrow><msub><mi>Y</mi><mrow><msub><mi>p</mi><mi>E</mi></msub></mrow></msub></mrow><mrow><msub><mi>Y</mi><mrow><msub><mi>n</mi><mi>E</mi></msub></mrow></msub></mrow></mfrac></mrow></mfrac></mrow></mtd><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow/></mtd></mtr></mtable></mrow></mrow><mn>.</mn></mrow></math> Although our ratio values are not continuous, the incident beam energy changes by 0.1 MeV, which is small compared to the resonance width larger than 0.2 MeV in this work. Therefore, Eq. 5 is acceptable. The resonances in the reactions can be identified by an extremely high cross section at a certain energy. If resonance occurs in the reaction at an incident energy E, the yield of the thin target <math id="M6"><mrow><mi>Δ</mi><msub><mi>Y</mi><mrow><msub><mi>N</mi><mi>E</mi></msub></mrow></msub></mrow></math> or <math id="M7"><mrow><mi>Δ</mi><msub><mi>Y</mi><mrow><msub><mi>P</mi><mi>E</mi></msub></mrow></msub></mrow></math> changes drastically, which may result in a maximum or a minimum value in the ratio. Because the △ Y is proportional to the cross section, the resonances of the two reactions can be identified by the ratios. 3 Results 1 Table 1 Table 1 lists the deuteron beam current intensity, total charge, beam time, and incident energy used in this work. The beam time of each run is 5 min, except for the first three runs. In addition, a 2 h background measurement was performed on the EJ-301 detectors without the deuteron beam just before the experiment, which ensured that the background measurements were performed in the same condition. 1 1 Incident energy (MeV) 1 1 Time (s) 1 1 Current (nA) 1 1 Charge (nC) 1 1 Incident energy (MeV) 1 1 Time (s) 1 1 Current (nA) 1 1 Charge (nC) 1 1 0.600 1 1 430 1 1 1.50 1 1 676 1 1 1.900 1 1 300 1 1 0.12 1 1 36 1 1 0.700 1 1 312 1 1 1.30 1 1 387 1 1 2.000 1 1 300 1 1 0.15 1 1 77 1 1 0.800 1 1 296 1 1 1.00 1 1 300 1 1 2.100 1 1 300 1 1 0.10 1 1 59 1 1 0.900 1 1 300 1 1 1.00 1 1 295 1 1 2.200 1 1 300 1 1 0.10 1 1 59 1 1 1.000 1 1 300 1 1 1.00 1 1 297 1 1 2.300 1 1 300 1 1 0.10 1 1 64 1 1 1.100 1 1 300 1 1 0.70 1 1 198 1 1 2.400 1 1 300 1 1 0.05 1 1 43 1 1 1.200 1 1 300 1 1 0.50 1 1 175 1 1 2.500 1 1 300 1 1 0.07 1 1 51 1 1 1.300 1 1 300 1 1 0.40 1 1 120 1 1 2.600 1 1 300 1 1 0.10 1 1 92 1 1 1.400 1 1 300 1 1 0.40 1 1 120 1 1 2.700 1 1 300 1 1 0.08 1 1 54 1 1 1.500 1 1 300 1 1 0.20 1 1 125 1 1 2.800 1 1 300 1 1 0.07 1 1 51 1 1 1.600 1 1 300 1 1 0.20 1 1 70 1 1 2.900 1 1 300 1 1 0.07 1 1 51 1 1 1.700 1 1 300 1 1 0.20 1 1 70 1 1 3.000 1 1 300 1 1 0.07 1 1 51 1 1 1.800 1 1 300 1 1 0.20 1 1 48 1 1   1 1   1 1   1 1   Figure 3 Figure 3 shows an example of the signals of neutron detector at 48°, when the incident energy is 2 MeV. The pulse shape discrimination (PSD) is calculated to distinguish neutrons and γ -ray. The PSD is given by, <math display="block" id="M8"><mrow><mi>P</mi><mi>S</mi><mi>D</mi><mo>=</mo><mfrac><mrow><msub><mi>Q</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>Q</mi><mtext>s</mtext></msub></mrow><mrow><msub><mi>Q</mi><mtext>l</mtext></msub></mrow></mfrac><mn>,</mn></mrow></math> where Q l is the total integral of a whole pulse, Q s is the partial integral from the star of the signal to the point when the signal falls to half the maximum. The different fall times of neutron signals and γ -ray signals in the EJ-301 detector make this method feasible. < > Neutron signals have a longer attenuation time than γ -ray signals in the neutron detector. The peak in Fig. 3 Fig. 3 with higher psd represents the neutron signals. The threshold of neutron detection is that the ADC of the neutron signals must be larger than 300 channels. %The neutrons yields are calculated by substract the background signals numbers from the neutron signals number of which psd above threshold. Figure 4 Figure 4 shows an example of the proton signals and the results simulated by SIMNRA at 135° with an incident energy of 2 MeV. The SIMNRA[ 19 19 ] is a software suitable for computer simulation of ion beam analysis methods with MeV ions. The Q value of this reaction is approximately 2.7 MeV, which helps identify and calibrate protons from the detector. The signals below 1600 channel are the Rutherford backscattered (RBS) deuterons recoiled from the target, and the signals above 2000 channel are protons. The energy distribution of proton is caused by the change in energy of deuteron beam in the thick target. The peak at channel 2774 represents the signals of protons when the incident deuteron energy is 2 MeV. Further, the signals between channel 2774 to channel 3555 belong to the protons produced by an incident energy of less than 2 MeV because of ionization in the thick target. The energy shift is caused by the center of mass motion. The Au-Si detector was placed at 135° in a direction opposite to the beam in the center of mass frame. Under these conditions, the higher the incident energy of the deuteron, the faster the center of mass moves, resulting in low speeds of outgoing protons in the backward direction. Figure 5 Figure 5 shows the yields of neutrons and protons when the incident energy is between 0.6-3 MeV. The yield is corrected by excluding the background and being converted into 300 s. The errors in yields depend on the statistical magnitude. Because every count in each measurement is larger than 560, the errors are less than 4.2%. Because of the sparking of the Faraday-cup and the imprecise incident beam intensity, the extreme value in the plot does not provide any reasonable conclusion about the cross section of the reactions. Our method enables exclusion of error in the incident beam and finds the physical information in the two reactions. 4 Discussion 1 The target used in this work is natural graphite, which contains 12 C and 13 C. The proton and the neutron channels from the reaction between 13 C and deuteron must be considered. The cross section of the 13 C(d,n) 14 N reaction is about two to three times larger than that of the 12 C(d,n) 13 N reaction, as indicated by the reference [ 14 14 , 15 15 ]. The abundance of 12 C is 98.89%, whereas that of 13 C is 1.109%. Although the 13 C has a larger cross section in this experiment than that of the 12 C, the large difference in abundance makes the neutron yield in the 13 C target, which is within 3.3% of the total neutron yield, insignificant. In addition, other possible reaction channels in the experiment, such as 12 C(d,n+p) 12 C, 12 C(d,n+d) 11 C, and 12 C(d,n+ α ) 9 B must be considered. The data from the TENDL-2017 Nuclear data library [ 16 16 ] indicate that the cross-section of these reactions from 0.6 MeV to 3 MeV are too small to be considered. For the energy range considered in this work, the dominant reactions in this experiment are 12 C(d,n) 13 N and 12 C(d,p) 13 C. In Fig. 6 Fig. 6 , the ratio is different in three neutron detection directions because of angular distribution. For low resonance states, the outgoing particles could interfere with each other, as well as with the coulomb force, resulting in anisotropy. Thus, the ratio in every direction could provide the physical information of nuclear reactions. The ratio between neutron and proton yields presents some insights into the resonance in 12 C(d,n) 13 N and 12 C(d,p) 13 C. The maximum and the minimum points in the curves in Fig. 6 Fig. 6 exhibit some features that indicate the complexity of nuclear reactions. The minimum points in our curves are attributed to either a sudden jump in proton yield or a sudden drop in neutron yield, as discussed in the Methods section. Thus, if the cross section of the neutron does not decrease at the resonance points, the cross section of the proton must reach a maximum value, which helps understand the resonance state for the proton channel. Fig. 7 Fig. 7 shows the cross sections of neutron and proton from Ref. [ 5 5 , 14 14 , 17 17 ]. The plot illustrates whether the cross of section of neutron decreases at the possible resonance energies. The possible resonances in Fig. 6 Fig. 6 can be checked using these data. The proton detection reported in Ref. [ 5 5 ] was performed at 170°, which is different from our detection. However, because the angular distribution does not influence the resonance recognition, Ref. [ 5 5 ] still works in the following analysis. For the 12 C(d,p) 13 C reaction, two resonances are found at 1.7 MeV and 2.5 MeV in all directions, and one resonance is found at 1.4 MeV in the 0° direction. In Fig. 7 Fig. 7 , the cross section of neutron does not decrease at 1.4 MeV, 1.7 MeV, or 2.5 MeV, which guarantees the resonance points for the proton channel. The resonance point at 1.7 MeV originates from the shell in the 14 N with an energy of 11.76 MeV [ 5 5 , 18 18 ]. The one at 2.5 MeV originates from the shell in the 14 N with an energy of 12.4 MeV [ 5 5 ]. Similarly, the one at 1.4 MeV is consistent with Ref. [ 5 5 ], and originates from the shell in the 14 N, whose energy is 11.51 MeV, and spin and parity are 3 + . Similarly, for the 12 C(d,n) 13 N reaction, resonances are found at 2.7 MeV at 0° and 24°. Further, two other resonances are found at 1.6 MeV and 1.9 MeV in all directions. In Fig. 7 Fig. 7 , the cross section of proton does not increase at 1.6 MeV and 2.7 MeV, which confirms a 1.6 MeV resonance point in 12 C(d,n) 13 N. On the contrary, the cross section at 1.9 MeV decreases significantly, indicating that this minimum value may not be a resonance. The resonance at 2.7 MeV is consistent with that of Ref. [ 4 4 , 14 14 ] originating from the excitation state of 14 N with an energy of 12.58 MeV. Further, the resonance at 1.6 MeV is consistent with that of Ref. [ 13 13 ]. The resonance at 1.6 MeV in the 12 C(d,n) 13 N reaction originates from the shell in the 14 N, whose excitation energy, and spin and parity are 11.66 MeV and 2 + or 3 + , respectively. Table 2 Table 2 lists the resonances in 12 C(d,n) 13 N and 12 C(d,p) 13 C reported by Ref. [ 4 4 , 5 5 , 13 13 , 14 14 , 17 17 , 18 18 ] and us. 5 1 12 C(d,n) 13 N 3 1 12 C(d,p) 13 C 1 1 Ref. [ 14 14 ] (MeV) 1 1 Ref. [ 4 4 ] (MeV) 1 1 Ref. [ 17 17 ] (MeV) 1 1 Ref. [ 13 13 ] (MeV) 1 1 This work (MeV) 1 1 Ref. [ 5 5 ] (MeV) 1 1 Ref. [ 18 18 ] (MeV) 1 1 This work (MeV) 1 1   1 1 1 1 0.85 1 1 1 1 1 1 0.89 1 1 1 1   1 1 0.945 1 1 1 1 1 1 0.91 1 1 1 1 1 1 0.94 1 1 1 1   1 1 1 1 1 1 1.16 1 1 1 1 1.16 1 1 1.18 1 1   1 1 1.279 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1   1 1 1 1 1 1 1.3 1 1 1 1   1 1 1.31 1 1   1 1   1 1 1 1 1 1 1 1 1 1 1.41 1 1 1.45 1 1 1.400 1 1 1.586 1 1 1 1 1.5 1 1 1 1   1 1   1 1   1 1   1 1   1 1 1 1 1 1 1.62 1 1 1.600 1 1 1.69 1 1 1 1   1 1 1.727 1 1 1 1 1 1 1.76 1 1   1 1 1 1 1.74 1 1 1.700 1 1 2.279 1 1 1 1 1 1 1 1   1 1   1 1   1 1   1 1   1 1 2.32 1 1 1 1 1 1   1 1   1 1   1 1   1 1   1 1 1 1 1 1 1 1 1 1 2.45 1 1 1 1 2.500 1 1 2.684 1 1 2.6 1 1 1 1 1 1   1 1   1 1   1 1   1 1   1 1 2.7 1 1 1 1 1 1 2.700 1 1   1 1   1 1   1 1   1 1 2.83 1 1 1 1 1 1   1 1   1 1   1 1   Figure 6 Figure 6 indicates that the ratio changes rapidly in the lower energy part; however, with the increase in energy, it becomes relatively more stable, making it difficult to recognize the resonance points. Although the yield change rate reflects the cross section, the higher the energy of the incident particles in the thick target, the higher the energy change of deuterons in the entire bombarding process; this leads to a more significant difference between the cross section and the ratio of yield. Under this scenario, only the great resonance points are found in our curves. This is another reason for the fewer number of resonances than in other studies. 5 Summary 1 In this study, by performing and studying the bombardment of deuteron on a thick carbon target, the resonances in the two reactions, 12 C(d,n) 13 N and 12 C(d,p) 13 C, were observed. The resonances when the incident deuteron energies were 1.4 MeV, 1.7 MeV, and 2.5 MeV in the 12 C(d,p) 13 C reaction, and the incident deuteron energies were 1.6 MeV and 2.7 MeV in the 12 C(d,n) 13 N reaction, were confirmed. We propose the ratio of neutron yield to proton yield as a new means to study the resonances in the 12 C(d,n) 13 N and the 12 C(d,p) 13 C reactions. This approach greatly reduces the possible systematic errors arising from both the fluctuation of the incident beam and the complex process of projectile in the thick target. However, some improvements are necessary in the future, because (a)although both neutron resonance and proton resonance points are found, their numbers are relatively small; (b) the feasibility of this method strongly depends on the incident energy. Moreover, the resonance could be difficult to recognize if the incident energy is very high or the incident energy range is very wide. The corresponding experiments can be improved by the following ways in the future: first, a thin target (approximately 1 μm thick) is recommended to replace a thick target, by which errors arising from the change in deuteron energy in the thick target can be eliminated. Second, the event-by-event method can be applied to eliminate the background signals. This method can be realized by placing a detector to detect the incident beam and performing a coincidence measurement among the incident particles and the outgoing particles. However, as a first step, our method ensures a low systematic error and provides some insights for more precise cross section measurements. References J.L. Vogl . Radiative capture of protons by 12C and 13C below 700 keV . California Institute of Technology , 1963 . Radiative capture of protons by 12C and 13C below 700 keV R.E. Brown . Experimental study of the 17O(p,γ)14N reaction and a calculation of the rate of this reaction in the CNO Cycle in stars . Physical Review 125 , 347 ( 1962 ). doi: 10.1103/PhysRev.125.347 http://doi.org/10.1103/PhysRev.125.347 Experimental study of the 17O(p,γ)14N reaction and a calculation of the rate of this reaction in the CNO Cycle in stars S. Typel , G. Baur . Theory of the TrojanHorse method . Annals of Physics 305 , 228 ( 2003 ). doi: 10.1016/S0003-4916(03)00060-5 http://doi.org/10.1016/S0003-4916(03)00060-5 Theory of the TrojanHorse method R.J. Jaszczak , R.L. Macklin , J.H. Gibbons . 12C(d,n)13N total cross section from 1.2 to 4.5 MeV . Physical Review 181 , 1428 ( 1969 ). doi: 10.1103/PhysRev.181.1428 http://doi.org/10.1103/PhysRev.181.1428 12C(d,n)13N total cross section from 1.2 to 4.5 MeV F. Papillon , P. Walter . Analytical use of the multiple gamma-rays from the 12C(d,p)13C nuclear reaction . Nucl. Instrum. Meth. B 132 , 468 ( 1997 ). doi: 10.1016/S0168-583X(97)00438-2 http://doi.org/10.1016/S0168-583X(97)00438-2 Analytical use of the multiple gamma-rays from the 12C(d,p)13C nuclear reaction X.D. Tang , S.B. Ma , F. Xiao et al ., An efficient method for mapping the 12C+12C molecular resonances at low energies . Nucl. Sci. Tech . 30 , 126 ( 2019 ). doi: 10.1007/s41365-019-0652-9 http://doi.org/10.1007/s41365-019-0652-9 An efficient method for mapping the 12C+12C molecular resonances at low energies M. Peng , G.Z. He , Q.W. Zhang et al ., Study of neutron production and moderation for sulfur neutron capture therapy . Nucl. Sci. Tech ., 30 : 2 . doi: 10.1007/s41365-018-0529-3 http://doi.org/10.1007/s41365-018-0529-3 Study of neutron production and moderation for sulfur neutron capture therapy A. Curtis , C. Calvi , J. Tinsley et al ., Micro-scale fusion in dense relativistic nanowire array plasmas . Nature Communications 9 , 1077 ( 2018 ). doi: 10.1038/s41467-018-03445-z http://doi.org/10.1038/s41467-018-03445-z Micro-scale fusion in dense relativistic nanowire array plasmas G. Zhang , M. Huang , A. Bonasera et al ., Nuclear probes of an out-of-equilibrium plasma at the highest compression . Physics Letters A 383 , 2285 ( 2019 ). doi: 10.1016/j.physleta.2019.04.048 http://doi.org/10.1016/j.physleta.2019.04.048 Nuclear probes of an out-of-equilibrium plasma at the highest compression J.Y. Li , R.J. Wang , Y. Zhang et al ., Calculation on the associated particle method for 90 to measure the neutron yield correction factors of D-D neutron generator . Nuclear Techniques 40 , 010202 ( 2017 ). doi: 10.11889/j.0253-3219.2017.hjs.40.010201 http://doi.org/10.11889/j.0253-3219.2017.hjs.40.010201 (in Chinese) Calculation on the associated particle method for 90 to measure the neutron yield correction factors of D-D neutron generator Y.G. Liu , M. Liu , J.L Ke et al ., Heat-flow-solid coupled analyses of the deuterium-target for compact D-D neutron generator Nuclear Techniques 40 , 010202 ( 2017 ). doi: 10.11889/j.0253-3219.2017.hjs.40.010202 http://doi.org/10.11889/j.0253-3219.2017.hjs.40.010202 (in Chinese) Heat-flow-solid coupled analyses of the deuterium-target for compact D-D neutron generator H.J. He , X.H. Xu , Y.K. Lu et al ., Study of deposition behavior of deuterium in palladium using D-D nuclear reaction Nuclear Techniques 40 , 020201 ( 2017 ). doi: 10.11889/j.0253-3219.2017.hjs.40.020201 http://doi.org/10.11889/j.0253-3219.2017.hjs.40.020201 (in Chinese) Study of deposition behavior of deuterium in palladium using D-D nuclear reaction T. W. Bonner , J. T. Eisinger , Alfred A. Kraus et al ., Cross section and angular distributions of the (d,p) and (d,n) reactions in 12C from 1.8 to 6.1 MeV . Physical Review 101 , 209 ( 1956 ). doi: 10.1103/PhysRev.101.209 http://doi.org/10.1103/PhysRev.101.209 Cross section and angular distributions of the (d,p) and (d,n) reactions in 12C from 1.8 to 6.1 MeV R.W. Michelnn , D. Meyer , K. Bethce , et al ., Excitation functions for the reactions 10B(d,n)11C and 12C(d,n)13N for charged particle aciwation analysis . Nuclear Instruments and Methods in Physics Research, Section B, (Beam Interactions with Materials and Atoms) 51 , 1 ( 1990 ). doi: 10.1016/0168-583x(90)90530-8 http://doi.org/10.1016/0168-583x(90)90530-8 Excitation functions for the reactions 10B(d,n)11C and 12C(d,n)13N for charged particle aciwation analysis C.R. Brune and R.W. Kavanagh . Total cross sections and thermonuclear reaction rates for 13C (d,n) and 14C(d,n) . Physical Review C 45 , 1382 ( 1992 ). doi: 10.1103/PhysRevC.45.1382 http://doi.org/10.1103/PhysRevC.45.1382 Total cross sections and thermonuclear reaction rates for 13C (d,n) and 14C(d,n) https://tendl.web.psi.ch/tendl_2017/deuteron_html/C/DeuteronC12xs.html https://tendl.web.psi.ch/tendl_2017/deuteron_html/C/DeuteronC12xs.html M. L. Firouzbakht , D. J. Schlyer , A. P. Wolf . et al ., Cross-section measurements for the 13C(p,n)13N and 12C(d,n)13N nuclear reactions . Radiochimica Acta 55 , 1 ( 1991 ). doi: 10.1524/ract.1991.55.1.1 http://doi.org/10.1524/ract.1991.55.1.1 Cross-section measurements for the 13C(p,n)13N and 12C(d,n)13N nuclear reactions L. Csedreki , I. Uzonyi , G. Szki , et al ., Measurements and assessment of 12C(d,pγ)13C reaction cross sections in the deuteron energy range 740-2000 keV for analytical applications . Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 328 , 59 ( 2014 ). doi: 10.1016/j.nimb.2014.02.123 http://doi.org/10.1016/j.nimb.2014.02.123 Measurements and assessment of 12C(d,pγ)13C reaction cross sections in the deuteron energy range 740-2000 keV for analytical applications M. Mayer . Simnra user's guide, report ipp 9/113 . Max-Planck-Institut fur Plasmaphysik, Garching, Germany , 1 ( 9 ): 9 , 1997 . doi: 10.1063/1.59188 http://doi.org/10.1063/1.59188 Simnra user's guide, report ipp 9/113 We thank Dr. Zhou-Tong He from Shanghai Institute of Applied Physics, Chinese Academy of Sciences, for the target preparation in this work. 10.1007/s41365-019-0705-0 This work was partially supported by the Strategic Priority Re-search Program of the Chinese Academy of Sciences (Nos. XDB16 and XDPB09), the National Natural Science Foundation of China (Nos. 11890714 and 11421505), and the Key Research Program of Frontier Sciences of the CAS (No. QYZDJ-SSW-SLH002). 2019-09-07 2019-10-11 2019-10-28 2019-11-20 2019-12-01 2021-04-11 © 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. 2019 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-019-0705-0 http://dx.doi.org/10.1007/s41365-019-0705-0 . 2019 Nuclear Science and Techniques 12 30