Introduction

For hundreds of years, efforts have been made to discover laboratory reactions to induce chemical changes, since the first alchemy trial, which aim to transmute base metal to gold. The new age of “modern alchemy” started with the rise of nuclear physics at the beginning of the 20th century1. Thanks to progress in nuclear reaction studies, making elements in a laboratory has become possible and great success has been achieved since then in creating elements in a laboratory by synthesizing atomic nuclei via various nuclear reactions, such as fission, fragmentation, and low-energy fusion2. However, these conventional reactions3 have difficulties in accessing proton-rich and neutron-rich isotopes because their production is very suppressed. The cross section for fission and fragmentation drops drastically for both proton-rich and neutron-rich nuclei. The target thickness is very limited in fusion due to its low reaction energy, although fusion has a large cross section to produce proton-rich nuclei. The question of how to effectively make nuclei becomes one of the main motivations of recent nuclear reaction studies, as well as of nuclear transmutation studies of radioactive waste4. Multi-nucleon transfer reaction has been recently proposed5 for the production of neutron-rich nuclei due to its large cross sections. In contrast, no good candidate with optimized energy has been found to overcome the limitation in the production of proton-rich nuclei.

An alternative reaction for producing proton-rich nuclei via incomplete fusion (ICF) induced by loosely bound nuclei has been proposed, and this possibility has been recently demonstrated for deuteron by theoretical studies6,7,8. To produce isotopes with increased charge (ΔZ = +1), the breakup of deuteron allows ICF; namely, absorption of the proton released during the breakup in addition to complete fusion (CF) of deuterons, as displayed in Fig. 1a–c. The compound process of ICF will lead to a large enhancement in the cross section at an energy higher than that of CF and thus may increase the target thickness and maintain a large cross section, resulting in high production of proton-rich nuclei. This scenario, however, has not yet been truly addressed experimentally. We note that ICF at energies near Coulomb barrier was intensively studied experimentally and theoretically9,10,11.

Fig. 1
figure 1

Schematic illustration of the different components in deuteron-induced reaction. a Complete fusion, where deuteron interacts with a colliding nucleus as a whole set. b Incomplete fusion (of breakup proton as an example), where deuteron is broken up, and proton is absorbed and neutron is emitted. c Elastic breakup, where deuteron disintegrates into proton and neutron while the breakup nucleons have no interactions or elastically interact with the colliding nucleus. In both (a) and (b), an excited compound nucleus is formed after absorption, and thus, other nuclei can be produced from the statistical decay of compound states by emission of γ rays, neutron (n), proton (p), alpha (α), and/or other light particles. In (c), the colliding nucleus remains unexcited and thus has no contribution to the production of other nuclei

The experimental challenge in the observation of such an ICF effect arises from the cross-section measurement requiring direct isotopic identification of the reaction products. To date, ICF has mainly been investigated by bombarding target materials with beams of loosely bound nuclei6,7,9,12,13,14,15. This method requires the accumulation of γ rays or α particles, and data are therefore limited to only a few isotopes. In addition, the cross sections on deuteron were measured up to only 60 MeV6,7, which is much lower than the predicted energy of the enhancement. To clarify the effect of ICF, we report direct measurements of the absolute production cross sections of the reaction products were performed for 107Pd on proton and deuteron at 50 MeV per nucleon (MeV/u) at the Radioactive Isotope Beam Factory (RIBF) operated by RIKEN Nishina Center and the Center for Nuclear Study, University of Tokyo.

The measurements were performed using the inverse kinematics technique, in which a radioactive beam of 107Pd, produced at the in-flight separator BigRIPS16, bombarded a proton/deuteron target, and the reaction products were identified directly and unambiguously via the ZeroDegree spectrometer16. To overcome the technical difficulties in the production of 50-MeV/u beams and identification of the reaction products involving a charge-state distribution and large energy loss, several modifications were used: an energy-control method in BigRIPS17, gas hydrogen and deuterium targets at high pressure, dispersive mode of ZeroDegree spectrometer16 and a large-area plastic18 and segmented ion chamber (IC)19. Details are provided in the section “Methods”.

Results

Particle identification for reaction products

An example of a particle identification plot indicating the reaction products obtained with the ZeroDegree spectrometer is presented in Fig. 2. The high resolving power for particle identification enables clear identification of ions in different charge states, namely, fully stripped (Q = Z), hydrogen-like (Q = Z − 1), and helium-like (Q = Z − 2).

Fig. 2
figure 2

Particle identification plots for the reaction residues produced by the 107Pd beam on the deuterium target and detected by the ZeroDegree spectrometer. a Two-dimensional plot of atomic number Z versus mass-to-charge ratio A/Q. The horizontal dashed lines show the Z gate that was adopted for selecting the Pd isotopes. b The A/Q spectrum for the Pd isotopes by selecting 45.5 ≤ Z ≤ 46.5, as indicated by the horizontal dashed lines in a. The fully stripped ions could be separated from other ones in different charge states. The red solid line represents the total fit resulting from two Gaussian shapes

Isotopic distribution cross section

The isotopic distribution of cross sections on proton (σp) and deuteron (σd) measured in the present work for 107Pd at 50 MeV/u are displayed together in Fig. 3. Each panel displays the results for one element, namely Fig. 3a for Ag (Z= 47), Fig. 3b for Pd (Z= 46), Fig. 3c for Rh (Z= 45), Fig. 3d for Ru (Z= 44), Fig. 3e for Tc (Z= 43), and Fig. 3f for Mo (Z= 42). The measurement for cross section was limited down to 1 mb due to statistics. The errors shown in Fig. 3 are from statistical contributions. The systematic error includes the target thickness (2%) and the charge-state distribution (3%). For the ΔZ = +1 reaction, both newly measured σp and σd are found to be much larger than their values at 100 MeV/u, in particular for 103,104Ag, and such cross-section increases are more significant than those observed from 200 to 100 MeV/u for both proton and deuteron20, as presented in Fig. 3a. Moreover, σd at 50 MeV/u becomes much larger than σp at 100 MeV/u for 103,104Ag, indicating a different behavior from heavier Ag isotopes, where these two values are similar to each other. The large σp value is probably because the total incident energy (E) of 50 MeV is close to that of the compound process of proton fusion21. For deuteron, the cross-section increase appears at E of 100 MeV, which is much higher than the typical value of the compound process for CF of deuteron, indicating a different mechanism.

Fig. 3
figure 3

Isotopic distribution of the measured cross sections for proton (circles) and deuteron (squares) at 50 MeV/u. (a) Ag, (b) Pd, (c) Rh, (d) Ru, (e) Tc, and (f) Mo. The errors shown in the figure are from statistical contributions. The cross sections on proton (crosses) at 100 MeV/u20 are displayed for comparison. The deuteron-induced cross sections (diagonal crosses) at 100 MeV/u20 for the Ag isotopes are also displayed

Effect of ICF on cross section

To investigate the nature of the cross-section increase, both σd and σp are plotted as a function of E in Fig. 4 for 103,104,105Ag. Typically, in the ΔZ = + 1 reaction, the cross section presents a maximum in a compound process after the absorption of nucleon(s) and a smooth decrease with E due to the emission of light particles during the pre-equilibrium decay22,23. Indeed, in the case of protons, the large σp(E) values at 50 MeV for 103,104Ag reflect a compound process after fusion21. For deuterons, a significant enhancement in σd(E) at 100 MeV is highlighted by our 104Ag data compared to the linear extrapolation at high energy, revealing the characteristic of absorption of nucleon(s). As discussed before, the absorption of breakup protons can produce ΔZ = +1 isotopes and lead to a second compound process in addition to CF6,7,8. Thus, the enhancement of σd(E) in the 104Ag results of the present work provides direct experimental evidence for ICF of breakup proton.

Fig. 4
figure 4

Proton-induced and deuteron-induced cross sections for the Ag isotopes produced by 107Pd as a function of total incident energy. (a) 105Ag, (b) 104Ag, and (c) 103Ag. Cross sections on proton and deuteron are displayed by circles and squares, respectively. Filled symbols indicate the results of the present work and open symbols indicate the high-energy data at 100 and 185 MeV/u20. It is noted that the error bars are smaller than the symbol sizes. Dashed, dot-dashed, and solid lines are DEURACS calculations of incomplete fusion (ICF), complete fusion (CF), and their sum (ICF+CF), respectively

Theoretical calculations of cross sections for 103,104,105Ag using a newly developed deuteron-induced reaction analysis code system (DEURACS)8,24 are presented in Fig. 4. In DEURACS, deuteron-induced cross section is defined as the incoherent summation of contributions from the CF of deuteron and ICF of proton and/or neutron released during the breakup. The formation of compound nuclei in CF and ICF are calculated independently by the Glauber model25. The pre-equilibrium and equilibrium decays of compound nuclei made in both CF and ICF are calculated by the exciton model and Hauser–Feshbach statistical model26, respectively. Nucleon optical potentials (OPs) of Koning and Delaroche27 are applied for breakup proton and neutron in ICF; OPs of An and Cai28 are used for deuteron in CF. As described by Nakayama et al.8, in order to obtain a reliable value with a wide energy range up to 200 MeV/u, the calculated isotopic-distribution cross sections are renormalized using the deuteron-induced total reaction cross sections given by the Minomo–Washiyama–Ogata formula29, which is based on a microscopic three-body reaction model. Details on the definition of CF and ICF in DEURACS are provided in the section “Methods”.

The calculations clearly show that the σd(E) enhancement in 103,104Ag corresponds to the compound ICF process. In addition, the calculations reveal that the energy of ICF increases with the number of evaporation neutrons in a similar manner to CF. This explains that the newly measured σd(E) for 105Ag corresponds to the pre-equilibrium process, resulting in a smooth increase from high energy to 100 MeV. For all Ag isotopes, the calculations indicate that ICF plays a primary role in σd(E) towards high energy because of its moderate decrease. The overall agreement between the DEURACS results and data provides further evidence that the σd(E) enhancement is a consequence of ICF of breakup proton.

It is very instructive to compare σd(E) with σp(E) for the different processes of ICF. In the pre-equilibrium process, σd(E) is found to follow a similar tendency as σp(E). Their common tendency is consistent with the two experimental findings for the ΔZ = +1 reaction at high energy: (1) σd is approximately half of σp at the same beam energy per nucleon, as displayed in Fig. 3a and reported in the measurements for 107Pd20, 93Zr30, 90Sr31, 137Cs31, and 136Xe32. (2) σp(E) is known to be inversely proportional to the total incident energy21,33. Considering that a deuteron has two nucleons, these two findings result in the relation σd(E) ≈ 1/2 × σp(E/2) ≈ σp(E) by using the total incident energy, supporting the observation of ~200 MeV in Fig. 4. Around the compound process, σd(E) is enhanced and becomes 4–5 times larger than σp(E), which remains in the pre-equilibrium process. This observation demonstrates the potential of loosely bound nuclei for nuclei production when the reaction energy is chosen for the compound ICF process.

Discussion

The isotopic measurement on σd provides further insight on the role of CF and ICF in the Ag production. As mentioned above, CF and ICF components are predominant at respective low and high energy regions for each Ag isotope, and both excitation functions move to the higher energy side when more neutrons are evaporated. Considering the incident energy is fixed to be 100 MeV in the present study, the σd data therefore correspond to ICF for the heavy-mass isotopes and access to CF for the light-mass ones as demonstrated in Fig. 4a–c. In another word, the σd data for different Ag isotopes indicate the heavy Ag ions are mainly produced by ICF and the role of CF increases towards the light ones.

Moreover, contributions from CF and ICF components change in different elements. For the products close to the projectile in mass, σd and σp are close to each other as presented in Fig. 3b, c for heavy Pd and Rh isotopes, respectively. For these isotopes, DEURACS indicates that contributions from CF and ICF are comparable to each other as discussed in ref. 8. When more nucleons are removed from the projectile, it is found that σd becomes larger than σp and the difference increases towards the neutron-deficient side within one isotopic chain as displayed in Fig. 3b–d. These light-mass ions are mainly produced via central collision, where the production cross section depends on the energy deposited. As discussed by Wang et al. 20,31 and Kawase et al. 30, the larger σd value is due to a higher energy deposit relative to proton because both two nucleons in deuteron contribute to the reaction indicating that CF is dominant to produce these light-mass ions. Indeed, the CF contributions are found to be more than one order magnitude larger than the ICF ones for Ru, Tc, and Mo as discussed by the DEURACS calculations by Nakayama et al. 8. Such picture is further supported by the fact that σd at 50 MeV/u and σp at 100 MeV/u are similar to each other for the light-mass ions as displayed in Fig. 3c–f. In the CF, deuteron interacts with nucleus as a whole set, and thus deuteron-induced reaction at 50 MeV/u and proton-induced one at 100 MeV/u dissipate similar energy, resulting in a similar production.

Our experimental observation of ICF provides an alternative path to produce ΔZ = +1 isotopes. Taking 104Ag as an example, the total incident energy of the ICF maximum is 100 MeV, while that for CF is predicted to be 50 MeV8. This energy difference enables an increase in the range of deuteron in the Pd material by more than a factor of 334. Considering a 50% reduction in cross section, the ICF component will lead to an increase of at least 50% in production compared to the CF component. Moreover, the ICF of deuterons may result in similar production as that obtained for a proton at 100 MeV even though the range of a deuteron is shortened by 40% compared to that of a proton because σd is almost four times larger than σp. The relative mass difference between a proton and deuteron is 2, which is much larger than the difference between heavy isotopes with Z > 1. Thus, a higher production yield may be expected for ICF of a loosely bound nucleus than from CF of other isotopes. This is crucial for fusion reactions involving unstable isotopes, which are normally loosely bound, and thus, is an important step in the search for a powerful tool in “modern alchemy”.

ICF has a strong impact on the optimization of deuteron beam energy for nuclear transmutation of long-lived fission products (LLFPs) in high-level radioactive waste. It has been discussed that the deuteron-induced spallation reaction could be a candidate for the transmutation of the LLFP nucleus 107Pd20. The transmutation of 107Pd may be beneficial for the energy of ICF because it maintains both a large cross section and a certain depth in the LLFP material. Moreover, the free neutrons produced during breakup may continuously interact with the LLFP material, leading to further transmutation, together with those emitted from the evaporation of compound nuclei. A detailed simulation including the effect of ICF is desired for nuclear engineering to evaluate LLFP transmutation by deuteron.

In conclusion, our results provide direct experimental evidence for ICF, which causes a σd(E) enhancement at the energy much higher than that for CF, and further emphasizes its potential in both the production of proton-rich nuclei and nuclear transmutation of LLFP nuclei. In addition, we have shown that the cross-section measurement via inverse kinematics is essential for verifying the ICF. Extending our knowledge of reaction mechanisms induced by loosely bound nuclei plays a key role in understanding element synthesis and the application of nuclear reactions. We expect that ICF induced by heavy loosely bound nuclei, such as the production of ΔZ = +2 isotopes via ICF induced by 9Be15 and 6He12,13, can be measured in the future.

Methods

Details on the experimental setup

The radioactive beams around 107Pd were produced from in-flight fission of U at 345 MeV/u. The average primary beam intensity was ~12 particle nanoamperes (pnA). The momentum acceptance of BigRIPS was set to 0.1% by using the slits at the dispersive focal plane F1. A special configuration of degraders and detectors in BigRIPS was used based on an energy-control method, and a detailed description can be found in ref. 17. The method enabled the beam energy to be 64 MeV/u in front of the secondary targets and to select only the N = 61 isotones around 107Pd. Particles in the secondary beams were identified event-by-event by time-of-flight (TOF) and magnetic rigidities (), as described in ref. 35. The intensity and purity of the 107Pd beam was 6.0 × 102 particles per second (pps) and 16%, respectively.

The 107Pd beam impinged on the gaseous hydrogen and deuterium targets, which were in a cylindrical cell with a diameter of 36 mm36. The temperatures of the hydrogen and deuterium gases were 32.7 and 36.0 K, respectively, and the pressure was 4 atmosphere (atm) for both. The average thicknesses of the proton and deuteron targets were 23.2 and 34.5 mg/cm2, respectively. The beam energy was 50 MeV/u in the middle of the secondary targets. Data collection was also performed using a cell without any target material (empty cell) to deduce the contributions from the beam-line materials.

The reaction residues were analyzed by the ZeroDegree spectrometer in dispersive mode, where the final focus F11 was momentum dispersive. In such mode, the momentum acceptance was 4%, and the horizontal and vertical angular acceptances were 40 and 60 mrad, respectively16. To cover a broad range of fragments and study the charge-state distribution, several different magnetic rigidity settings were applied: −6%, −4%, −2%, 0%, +2%, +4%, and +6% relative to the value of the beam. At F11, a new plastic scintillation detector18 with a large effective area of 240 × 150 mm2 was developed to cover a wide momentum acceptance, and a new segmented IC19, which had a length of 750 mm along the beam direction and consisted of 15 units for the readout pad, was developed. The reaction products were stopped in the IC.

The reaction products were identified by measuring TOF, , ΔE, and total energy E. The TOF information was obtained using two plastic scintillators at the entrance of ZeroDegree spectrometer (F8) and F11, and a trajectory reconstruction was applied to determine the value using position and angle measurements at F11. Both ΔE and E were deduced from the IC19. Atomic number Z was determined using the TOF, ΔE, and E information, and the mass-to-charge ratio A/Q was obtained from the TOF– correlation. The detection efficiencies of the detectors for particle identification were higher than 95%. The resolution in Z and A/Q for the Pd isotopes was 0.60 [given here as full width at half maximum (FWHM)] and 3.7 × 10−3 (FWHM), respectively. Due to the high A/Q resolution of the ZeroDegree spectrometer, ions in the fully stripped, hydrogen-like, and helium-like charge states were clearly identified (Fig. 2b), and their respective fractions were 12%, 41%, and 43%, respectively, for Z = 46 isotopes in the deuterium-target case.

The σp and σd values were obtained by using the data for the hydrogen and deuterium targets after subtraction of the contributions of the beam-line materials using the empty-cell run.

Details on the DEURACS calculations

The CF cross section is defined as

$$\sigma _{{\mathrm{CF}}}\left( {E_{\mathrm{{d}}}} \right) = \sigma _{{\mathrm{d - ABS}}}\left( {E_{\mathrm{{d}}}} \right) \times \frac{{\sigma _{{\mathrm{A(d,x)B}}}\left( {E_{\mathrm{{d}}}} \right)}}{{\sigma _{{\mathrm{A(d,rea)}}}\left( {E_{\mathrm{{d}}}} \right)}}$$

where \(\sigma _{{\mathrm{d - ABS}}}\) represents the cross section for the absorption of deuteron. \(\sigma _{{\mathrm{A(d,x)B}}}\) is the production cross sections of nucleus B from deuteron-induced reaction on nucleus A, and \(\sigma _{{\mathrm{A(d,rea)}}}\) accounts for the deuteron total reaction cross section on A.

The ICF cross section for the breakup proton is defined as

$$\sigma _{{\mathrm{ICF - p}}}\left( {E_{\mathrm{{d}}}} \right) = \sigma _{{\mathrm{p - ABS}}}\left( {E_{\mathrm{{d}}}} \right) \times {\int} {{\mathrm{{d}}}E_{\mathrm{{p}}}f\left( {E_{\mathrm{{p}}}} \right)\frac{{\sigma _{{\mathrm{A(p,x)B}}}\left( {E_{\mathrm{{p}}}} \right)}}{{\sigma _{{\mathrm{A(p,rea)}}}\left( {E_{\mathrm{{p}}}} \right)}}}$$

where \(\sigma _{{\mathrm{{p}}} - {\mathrm{ABS}}}\left( {E_{\mathrm{{d}}}} \right)\) represents the cross section for the breakup process, where proton is absorbed, and neutron is emitted. \(f\left( {E_{\mathrm{{p}}}} \right)\) is the normalized energy distribution of the absorbed proton and it is obtained kinematically from the angular-integrate spectra of the neutron emitted during the breakup, which is calculated with the Glauber model. \(\sigma _{{\mathrm{A(p,x)B}}}\) is the production cross sections of nucleus B from proton-induced reaction on nucleus A, and \(\sigma _{{\mathrm{A(p,rea)}}}\) is the proton total reaction cross section on A. The complete description on DEURACS is provided in ref. 8.