# Laser-initiated primary and secondary nuclear reactions in Boron-Nitride

Christine Labaune, C. Baccou, V. Yahia, C. Neuville, J. Rafelski

► To cite this version:

Christine Labaune, C. Baccou, V. Yahia, C. Neuville, J. Rafelski. Laser-initiated primary and secondary nuclear reactions in Boron-Nitride. Scientific Reports, 2016, 6, pp.21202. 10.1038/srep21202. hal-01284331

HAL Id: hal-01284331

https://hal.sorbonne-universite.fr/hal-01284331

Submitted on 7 Mar 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# Laser-initiated primary and secondary nuclear reactions in Boron-Nitride

Received: 15 June 2015

Accepted: 10 November 2015

Published: 17 February 2016

C. Labaune<fnref n="1" />, C. Baccou<fnref n="1" />, V. Yahia<fnref n="1" />, C. Neuville<fnref n="1" /> & J. Rafelski²

Nuclear reactions initiated by laser-accelerated particle beams are a promising new approach to many applications, from medical radioisotopes to aneutronic energy production. We present results demonstrating the occurrence of secondary nuclear reactions, initiated by the primary nuclear reaction products, using multicomponent targets composed of either natural boron (B) or natural boron nitride (BN). The primary proton-boron reaction ($p + ^{11}\text{B} \rightarrow 3\alpha + 8.7 \text{ MeV}$), is one of the most attractive aneutronic fusion reaction. We report radioactive decay signatures in targets irradiated at the Elfie laser facility by laser-accelerated particle beams which we interpret as due to secondary reactions induced by alpha ($\alpha$) particles produced in the primary reactions. Use of a second nanosecond laser beam, adequately synchronized with the short laser pulse to produce a plasma target, further enhanced the reaction rates. High rates and chains of reactions are essential for most applications.

The proton-boron reaction is one of the most interesting candidates for fusion energy production as it produces less than 1% of the energy in neutrons, so it is called aneutronic fusion. However, it seems unreachable with conventional fusion approaches because it requires a very high temperature<fnref n="1" />, and consequently a very high laser energy, and it seems difficult to achieve a positive energy balance taking into account the nuclear energy release and the radiation losses². Radically different approaches and schemes must be sought³,⁴. Taking advantage of the fast evolution of the short, high-intensity lasers⁵, we proposed a new scheme⁶ based on fusion reactions initiated by a laser-accelerated proton beam⁷,⁸ and we demonstrated its high efficiency compared to previous approaches⁹,¹⁰. Nonetheless, it is clear that it will not be possible to reach a positive energy balance without cycles of reactions to increase the overall fusion yield. The present work demonstrates the possibility of producing reactions chains initiated by charged particles using multi-component targets irradiated by laser-accelerated particles.

## Results and Discussion

**Nuclear pB, pN and secondary α-induced reactions.** The proton-boron fusion reaction has been extensively studied and produces three alpha particles through three main channels¹¹:

$$p + ^{11}_{5}\text{B} \rightarrow \alpha_{0} + ^{8}_{4}\text{Be} + 8.59 \text{ MeV} \rightarrow \alpha_{0} + \alpha_{01} + \alpha_{02} \quad (1)$$

$$p + ^{11}_{5}\text{B} \rightarrow \alpha_{1} + ^{8}_{4}\text{Be}^{*} + 5.65 \text{ MeV} \rightarrow \alpha_{1} + \alpha_{11} + \alpha_{12} \quad (2)$$

$$p + ^{11}_{5}\text{B} \rightarrow 3\alpha + 8.68 \text{ MeV} \quad (3)$$

The predominant channels go through the formation of a compound nucleus ¹²C* followed by sequential α decays via an unbound ⁸Be in its fundamental (1) or first-excited (2) state. Channel (3) corresponds to the direct 3-body reaction that, for low-energy protons, contributes to less than 5% to the total cross-section¹². The alpha particle spectrum that emerges derives from these reactions and is between a few hundred keV to ~10 MeV given the fusion yield and available proton energy¹²,¹³, with a strong wide peak near to 4 MeV, with location somewhat dependent on proton energy. A fourth γ capture fusion path exists where the ¹²C* compound nucleus decays by γ emission, releasing ~15.9 MeV, but in comparison this is a highly improbable reaction. The cross-section of the

<footnote n="1">¹LULI, Ecole Polytechnique, CNRS, CEA, UPMC, 91128 Palaiseau, France. ²The University of Arizona, Tucson, AZ 85721-0081, USA. Correspondence and requests for materials should be addressed to C.L. (email: christine.labaune@polytechnique.fr)</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Scheme of the set-up to study interaction between a proton beam accelerated by a short, high-intensity laser pulse and a solid or plasma boron or boron-nitride target. The two laser pulses have a well-defined timing with a time-delay as shown in the insert. The particles diagnostics are displayed around the irradiated target.</figcaption></figure>

$^{11}$B(p, α)$^{2}$α reaction as a function of the centre of mass energy is characterized by two low-energy resonances at E$_{c.m.}$ = 148 keV and E$_{c.m.}$ = 614 keV$^{11,14}$, that are very interesting for producing high numbers of reactions.

Other nuclear reactions between the accelerated protons and the boron target can occur. In the experiments described in Methods, a schematic view of the set-up is shown in Fig. 1, we used natural boron that is comprised of two isotopes $^{11}$B (80.1 ± 0.2%) and $^{10}$B (19.9 ± 0.2%). The two other primary nuclear reactions of interest in this context are$^{15,16}$:

$$p + ^{11}\text{B} \rightarrow ^{12}\text{C}^{*} \rightarrow ^{11}\text{C} + n - 2.765 \text{ MeV} \quad (4)$$

$$p + ^{10}\text{B} \rightarrow ^{11}\text{C}^{*} \rightarrow ^{11}\text{C} + \gamma + 8.69 \text{ MeV} \quad (5)$$

The cross sections for these two reactions reach significant values for proton energy above 3 and 1.5 MeV respectively. However, reaction (5) is ~one hundred times less probable than reaction (4). The radioactive $^{11}$C has a half-life of (20.364 ± 0.014) min by positron decay:

$$^{11}\text{C}_{\tau=20.36 \text{ min}} \rightarrow ^{11}\text{B} + e^{+} + 0.96 \text{ MeV} \quad (6)$$

Finally, the low-energy protons can also interact with $^{10}$B to produce $^{7}$Be by the reaction:

$$p + ^{10}\text{B} \rightarrow ^{11}\text{C}^{*} \rightarrow \alpha + ^{7}\text{Be} + 1.147 \text{ MeV}, \quad (7)$$

which produces alpha particles. $^{7}$Be decays via electron capture with a half-life of 53.24 days.

Secondary nuclear reactions can occur between the *in situ* produced α-particles, generated by the (1-3) reactions and $^{10}$B$^{17,18}$:

$$\alpha + ^{10}\text{B} \rightarrow ^{14}\text{N}^{*} \rightarrow ^{13}\text{C} + p + 3.85 \text{ MeV} \quad (8)$$

which regenerates high-energy protons, and

$$\alpha + ^{10}\text{B} \rightarrow ^{14}\text{N}^{*} \rightarrow ^{13}\text{N} + n + (1/3.51 \text{ MeV, depending on the channel}) \quad (9)$$

which requires alpha energy larger than 2.5 MeV to reach significant cross section values. The radioactive $^{13}$N has a half-life of 9.97 minutes by positron decay:

$$^{13}\text{N} \rightarrow ^{13}\text{C} + e^{+} + v_{e} + 2.22 \text{ MeV} \quad (10)$$

A synthetic schematic of the primary and secondary reactions of interest in our study, initiated by the interaction of protons irradiating a natural boron target, is shown in Fig. 2a.

The novel idea of these experiments was to use compound targets to take advantage of both the wide spectrum of energy and high intensity of the laser-accelerated proton beam to initiate chains of reactions based on complementary elements and of the secondary reactions driven by the alpha particles. Some of these are producing protons that promise the continuation of chain reactions and in principle the development of energy producing devices. We thus study a reaction cycle on a multicomponent, e.g. BA, target in which the primary reactions between the protons and the component B produces an α-particle which then reacts with the component A: p + B → α + X, α + A → p + Y.
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Scheme of the primary and secondary nuclear reactions produced by the interaction between a laser-accelerated proton beam and (a) a natural boron target, (b) a boron-nitride target. In the case of the BN targets the reactions with $^{10}$B can also occur but are not shown for clarity.</figcaption></figure>

We looked for an isotope that in reactions with α-particles produces protons and does not produce neutrons and in irradiation with either α-particles or protons does not go into undesirable side chains of reactions. A weakly (1–2 MeV) endothermic reaction is fine if this keeps up the chain of reactions. Allowing for endothermic reactions produced, we found one strongly aneutronic candidate, nitrogen, which further is easy to handle combined with boron. So we chose to use boron-nitride (BN) targets. Natural nitrogen is comprised of two isotopes $^{14}$N (99.636%) and $^{15}$N (0.364%). The primary and secondary nuclear reactions between the proton beam and the boron-nitride target are shown in Fig. 2b. The key α-induced secondary reaction of interest for our purpose is$^{19,20}$:

$$ \alpha + ^{14}\text{N} \rightarrow ^{18}\text{F}^{*} \rightarrow ^{17}\text{O} + \text{p} - 1.192 \text{ MeV} \qquad (11) $$

which produces energetic protons since the energy of the α-particles transfers to a large extent to the produced proton; low-energy α cannot react considering the relatively strong Coulomb repulsion. The cross section is indeed resonant and peaked at a few MeV, just where the fusion product α-particles appear, and therefore the protons produced are energetic.

The other α-induced secondary reaction$^{21}$:

$$ \alpha + ^{14}\text{N} \rightarrow ^{18}\text{F}^{*} \rightarrow ^{18}\text{F} + \gamma + 4.415 \text{ MeV} \qquad (12) $$

is of interest in medical applications given the radioactive $^{18}$F production. $^{18}$F has a half-life of 109.8 minutes by positron emission (97%) and by electron capture (3%) and will be used to diagnose this reaction.

Furthermore, $^{11}$C can also be produced by the following reaction:

$$ \text{p} + ^{14}\text{N} \rightarrow ^{15}\text{O}^{*} \rightarrow ^{11}\text{C} + \alpha - 2.92 \text{ MeV} \qquad (13) $$

However to achieve an effective reaction rate of a few per 10,000 in solid target a proton energy larger than 6 MeV is needed$^{22}$.

More nuclear reactions are possible in this context, but we do not consider them in a first approach because either their cross section is too small at the involved energies or the proportion of present elements is too small to contribute significantly to the observations.

**α-particles characterization.** The α-particles were observed and characterized with different diagnostics. An example of the high-energy spectrum measured in the magnetic spectrometer for a solid boron target is shown in Fig. 3a. The low-energy cutoff of the spectrum is due to the Aluminum filter with 6 µm thickness placed in front of the spectrometer. ~$10^8$ sr$^{-1}$ α-particles with energy between 2 and 7 MeV are observed. Low-energy α-particles, between 500 and 900 keV, were observed in the Thomson parabola, but these particles are of less interest for the studies presented here.

α-particles were also observed with the boron-nitride targets. The spectrum of the high-energy particles is similar to that of the boron target but we observed a larger number of α-particles with BN both in the magnetic spectrometer and in the CR39 track detectors placed at different angles. These detectors were covered with 6 and
<figure><img src="image_4.png" /><figcaption>**Figure 3.** (a) Spectrum of the high-energy component of the α-particles produced by the interaction of a laser-accelerated proton beam and a boron target. The red line in the spectrum corresponds to the cutoff of the 6 µm Aluminum filter in the magnetic spectrometer. The error bars take into account the variation of the number of impacts across the direction perpendicular to the magnetic field. (b) CR39 images for B and BN targets showing the large impact diameter associated with the α-particles and the small impact diameters associated with protons.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** (a) Activation curves as a function of time for B and BN targets. The inset of part (a) is an enlargement of the beginning of the curve for the B target. (b) Fit of the BN activation curve using the sum of the contributions of ¹¹C and ¹⁸F. A Geiger-Mueller tube detector with a precision of ±5% was used, performing a series of quick measurements comprising 4–16 readouts resulting in activity measurement error that was limited by the statistical count and timing width of about half to one minute. The final measurement error for the radio activity was thus in the range of 2–10% depending on yield and activity lifespan.</figcaption></figure>

10 µm of Aluminum and so the observed particles have an energy higher than 2 and 3 MeV respectively. Typical images of CR39 covered with 6 µm of Aluminum for B and BN targets are shown in Fig. 3b.

**Activation measurements.** Examples of the time evolution of the normalized activation for B and BN targets are shown in Fig. 4. In both cases, the average slope in Logarithm scale, for t < 150 minutes, corresponds to the half-life of the ¹¹C, (20.4 ± 0.2) minutes. A very interesting feature appears at late times in the BN curve. For t > 90 min, in the case of BN, we observe a clear reduction of the slope of the curve. This is not observed in
<figure><img src="image_6.png" /><figcaption>Figure 5. Characterization of the proton beam produced by the short, high-intensity (1–2 × 10¹⁹ W/cm²) laser pulse with various targets using Thomson parabola. (a) Image of the Thomson parabola traces.</figcaption></figure>

(b) Integrated normalized proton flux per unit of solid angle, comprising all particles above energy E, as a function of the energy as deduced from the Thomson parabola analysis. The high-energy cutoff is modified by the conditions of interaction: nature of the target, intensity of the short pulse.

The B target case. By analyzing the curve, we found that the slope at the latest times was close to 100 minutes. We attributed this change to ¹⁸F whose contribution begins to appear when ¹¹C has already decreased in a significant way because the half-life of ¹⁸F is ~5 times longer than the half-life of ¹¹C. We calculated a fit of this curve using the sum of the two exponentials,

$$A(t) = A_{(^{11}\text{C})} \times \exp[-(t/29,336)] + A_{(^{18}\text{F})} \times \exp[-(t/158,36)], \quad (14)$$

where A(¹¹C) and A(¹⁸F) are the ¹¹C and ¹⁸F activations respectively. In this particular case, we found a good fit for A(¹¹C) = 500 Bq and A(¹⁸F) = 3.33 Bq, and interpreted the measured activation as the superposition of the ¹¹C and ¹⁸F decays. This corresponds to ¹¹C and ¹⁸F yields of 10⁶ and 6 × 10³ respectively. When analyzing more shots, we found that the ¹⁸F yields varied between <10³ (limit of detectability of our system) and 6 × 10³, which correspond to the error bars of our measurements.

At early times, t < 30 min, we observe that the activations measured in the experiment decrease faster, as a function of time, than the theoretical decay of ¹¹C. This is shown in the insert of Fig. 4a. This must be due to the contribution of a short-lived component that, for our conditions, could be Nitrogen-13. We observed that its proportion was larger in the B than in the BN targets, but it was not possible with the actual precision of our instrumentation to deduce an accurate yield of ¹³N. This is due to the fact that the half-lives of ¹³N and ¹¹C are too close and the access to the target takes ~14 minutes after the shot, so early essential measurements were not possible. Considering the possible nuclear reactions, the ¹³N is most likely produced by reaction (9), ¹⁰B(α, n)¹³N. As there are twice as many ¹⁰B in the B compared to the BN targets, this explains the observed higher ¹³N yields in B compared to BN targets.

We made a few tests, either by changing the proton energy spectrum or the irradiated target, to unequivocally identify the involved nuclear reactions. The proton spectrum was modified by changing the thin foil (thickness or nature) irradiated by the picosecond beam or by varying the intensity of the pulse. We first consider reactions producing ¹¹C in B targets. As a first test, we irradiated natural B targets with a proton beam accelerated out of an Aluminum foil of 10 µm thickness for which the maximum proton energy is ~3–3.4 MeV (see Fig. 5). The activation level was reduced by ~5 compared to the one obtained with a 20 µm Aluminum foil. A similar effect was observed by using Al-20 µm targets with a reduced intensity of the picosecond laser, whose main consequence was to remove the high-energy part of the proton spectrum but also to reduce the low-energy part. The activation level was then reduced by a factor ten. In both cases, as expected from its cross section¹⁵, the efficiency of the reaction ¹¹B(p, n)¹¹C was reduced, explaining the observed reduction of the activity. To check the efficiency of the ¹⁰B(p, γ)¹¹C reaction, we irradiated ¹⁰B targets with protons produced by an Al-20 µm foil at maximum intensity. The total activation was reduced by ~60 compared to the one obtained with natural boron, demonstrating the very low efficiency of the reaction (5) compared to (4) in agreement with their respective cross sections.

**Modification of the ¹¹C production in a plasma.** An interesting part of the experiment was comparing the results from B and BN targets where both solid and plasma targets were explored. Typical maximum activation levels with solid B were ~ (1000 ± 50) Bq, which correspond to ~2 × 10⁶ atoms of ¹¹C at t = 0 (time of the shot). The contributions of ¹³N and ¹⁸F are considered as negligible in a first approach. In BN targets, for all shots in solid, the activation levels were below those of B targets by a factor ~2. Most likely, this result simply reflects the fact that the number of active atoms is reduced by a factor 2 in BN targets compared to B targets. This also
<figure><img src="image_7.png" /><figcaption>**Figure 6.** (a) Ratio of the activation extrapolated to $t=0$ for a BN target, A(BN), compared to a B target, A(B), as a function of the time delay between the two laser pulses. A delay of $\Delta t = -100$ ps means that the short pulse was sent 100 ps before the long pulse. (b) Time of arrival of the protons on the target as a function of their energy for a distance of 6 mm between the two targets. The rectangles show the plasma forming pulses for the three time-delays.</figcaption></figure>

confirms that the reaction (13) is ineffective in agreement with its threshold of 6 MeV and the energy spectrum of the laser-accelerated proton beam. So, most of the observed ¹¹C was produced by reaction (4).

The surprising result was the increase of activation in the case where the BN targets were conditioned in the plasma state by the nanosecond beam. From the temporal analysis, the activation corresponds again mainly to ¹¹C production. In previous studies⁶, we already observed that the number of nuclear reactions was enhanced by plasma formation of the irradiated target. The interesting fact in this experiment is that the increase of activation is significantly greater in BN than in B plasmas. In Fig. 6a, the ratio of activation in BN versus B targets, A(BN)/A(B), is reported as a function of the time delay between the short and the long laser pulses. A delay of $\Delta t = -100$ ps means that the short pulse was sent 100 ps before the long pulse. Taking into account the distance of 6 mm between the two targets, this means that the most energetic protons arrive in the B/BN target when it is just beginning to be transformed into plasma. We observe that A(BN)/A(B) decreases from ~1.8 to 0.9 when $\Delta t$ varies from -100 ps to +500 ps and is always larger in plasma than in solid where it is ~0.5. Using a gold-10µm accelerated proton beam further increased, by ~50%, the ratio of ¹¹C production in BN in comparison to B plasmas. This latter result is likely due to the reaction (13) becoming active because of higher energy protons available with gold target.

## Discussion

Looking at the possible reactions in our experiment (Fig. 2), we interpret the observation of Fluorine-18 and Nitrogen-13 in the irradiated targets as the existence of secondary reactions driven by the α-particles produced by the primary reactions. In our experiments the α-particles are mainly produced by the ¹¹B(p,α)²α, (1-3), which in turn react with ¹⁴N through the reaction (12), ¹⁴N(α,γ)¹⁸F, to produce the observed Fluorine-18 and with ¹⁰B through the reaction (9), ¹⁰B(α,n)¹³N to produce the observed Nitrogen-13. So the first step of our experimental proof of concept was successful, showing that it is possible to produce secondary reactions involving α-particle interaction with a second component of the target. In the same way, secondary protons are likely produced by reactions (8) and (11) with an energy large enough to induce new reactions with ¹¹B.

In our BN targets, the observation of ¹⁸F radioactivity can be used as a diagnostic tool to estimate the yield of secondary protons produced in reaction (11): typically for each reaction of the type (12), where an observable radioactive isotope ¹⁸F is produced, there are up to 500 reactions (11) that produce an energetic proton. The exact ratio between (11) and (12) for in situ injected α-particle fusion spectrum in plasma environment where collective effects can be present is not known and could be further studied in future laser-plasma experiments. A similar diagnostic role is played by ¹³N activity. It is generated in reaction (9) which occurs also as a fraction of reaction (8)¹⁸. Additional experiments are needed to establish this ratio reliably, but one can presume that it is close to unity for the range of α-particle spectra available in situ.

In order to estimate the relative yield of ¹¹C produced in primary proton-induced reactions as compared to secondary ¹⁸F and ¹³N we need to remember that effective production of ¹¹C requires p at energy above 3 MeV while secondary α are products of pB fusion which occurs at much lower proton energy. The proton beam spectrum shown in Fig. 5 suggests that fusion is about a factor 100-1000 favored over ¹¹C production. However, non-laser, non-plasma beam stopping power wisdom suggests that primary reactions dominate secondary by a
factor 10,000. In addition we use ¹⁸F as diagnostic, which is suppressed by its electromagnetic production process, thus one can easily argue that one ¹⁸F is suppressed by a factor million, and this suppression is reduced to 1000–10,000 by the structure of the primary proton beam. Should ¹⁸F at a higher yield than ¹¹C/1000 be observed, this signals that secondary reactions are enhanced by plasma environment. We conceived this experiment in order to establish the strength of this effect.

The second and equally important point is the observation of the modification of the ¹¹C production in BN plasma as compared to B plasma. This result means that, beyond the plasma action, new effects occur that we attributed to the reaction cycles induced by the secondary reactions described above. It appears that in the plasma state the reaction rates are increased and this also applies to the secondary reactions. The cycle of reactions ¹¹B(p,α)2α followed by ¹⁴N(α,p)¹⁷O reiterates the reactions with ¹¹B, including the reaction ¹¹B(p,n)¹¹C producing the additional ¹¹C. This hypothesis also explains the increased number of α-particles in the case of BN compared to B targets. The effect of the time delay between the two laser pulses further demonstrates that the temperature and the ionization state of the target affects the nuclear reaction rates. In Fig. 6b, the time of arrival of the protons in the plasma as a function of their energy is shown for the three time-delays. For the case, –100 ps, all the protons of interesting energy (> 500 keV) interact with a hot plasma. For the case, +200 ps, the protons with energy below 700 keV arrive after the end of the plasma formation pulse and in a plasma with a lower temperature than in the previous case. For the case +500 ps, all the protons arrive after the end of the plasma formation. This means that the plasma entered the recombination phase with a temperature which decreases rapidly. The results are nevertheless different from those of the solid state. All this is consistent with and explains the observed decrease of A(BN)/A(B) as a function of increasing time delay.

## Conclusion

In summary, our results demonstrate that it is possible to produce a charged particle chain of reactions using boron-nitride targets irradiated by a laser-accelerated proton beam. We believe that this result is an important stepping stone towards energy production by laser fusion as well as to the production of required yields of medical radio-isotopes by laser. Our experimental conditions were not optimized for the ¹⁸F and ¹¹C productions, which explains the relatively low yields that are observed. Our results demonstrate the proof of principle that i) secondary reactions with the α-particle are possible in the target and ii) that the rates of the desired reactions can be significantly increased by the chain of reactions. In the future, improvements can be expected by optimization of isotope mixture, of the proton spectra (as was observed with gold), of the plasma formation and of the target characteristics.

## Methods

**Experimental set-up and laser-accelerated proton beam.** Experiments were conducted on the Elfie laser facility at Ecole Polytechnique. The scheme of the experimental set up is shown in Fig. 1. The short laser pulse, with 10 J in 0.35 ps at wavelength 1.06 µm, was focused with an off-axis parabolic mirror at normal incidence on a thin foil at an intensity of $1-2 \times 10^{19} \text{W.cm}^{-2}$, to produce a proton beam by the Target Normal Sheath Acceleration (TNSA) mechanism$^{7,8}$. We used three types of target: Aluminum of 10 and 20 µm thickness and Gold with 10 µm thickness. The energy spectrum of the proton beam was characterized using Thomson parabola and CR39 track-detector covered with Aluminum foils of different thickness to select the proton energy ranges. The integrated proton spectra for the three types of target are shown in Fig. 5 for two laser intensities. We observe that the low-energy part of the spectrum is similar for Al-20 µm and Au-10 µm. It is reduced for Al-10 µm or for lower short pulse intensity. The high-energy cutoff is larger for gold (~6.8 ± 0.5 MeV) than for Aluminum of 20 µm (~5.8 ± 0.5 MeV) and Aluminum of 10 µm (3.3 ± 0.5 MeV). This is in agreement with the CR39 measurements as well as with previous studies of TNSA proton acceleration. The integrated number of protons rises from $10^9 \text{sr}^{-1}$ at high energy to ~$10^{13} \text{sr}^{-1}$ at low energy. With reduced laser intensity, the energy cut-off is lowered and the integrated number of protons significantly reduced.

The proton beam was sent on a boron (B) or boron-nitride (BN) thick foil (500 µm-1 mm) placed at 20° from the proton axis at a distance of 6 mm from the proton target. In both cases, natural Boron was used, including the two isotopes ¹⁰B and ¹¹B. For part of the shots, the proton beam interacted directly with the target initially in the solid state and, for the other part, the target was previously turned into plasma by the nanosecond pulse. This pulse delivered 80 J in 0.6 ns at 1.06 µm leading to an intensity of $10^{13} \text{W/cm}^2$ on target. The short and the long pulses come from the same oscillator so a perfect synchronization was easily provided and different time-delays between the two laser pulses were used. The electron temperature of this plasma was ~0.3–0.5 keV.

**Diagnostics.** Particle diagnostics included a magnetic spectrometer, two Thomson parabolas and CR39 track-detectors. In purpose, the magnetic spectrometer and the Thomson parabola were set to collect α-particles in different ranges of energy, high- and low-energy ranges respectively. The diagnostics were installed in the horizontal plane at different locations as shown in Fig. 1. The CR39 were previously calibrated in the laboratory$^{23,24}$ and were used to distinguish the protons from the alpha particles by the size of the impact diameter, ~($1 \pm 0.5$) µm for protons and ~($2.8 \pm 0.5$)µm for α-particles. However, all these diagnostic led to an underestimation of the number of particles because the only ones that can be detected are those that come out of the target$^{25}$. So, to complete these data, we implemented activation measurements of the target after the shot as a function of time with a Geiger-Mueller tube detector with a precision of ±5%. The activation measurements allow the determination of the amount of isotope formed in the target and so provide the number of the studied reactions without depending on the output capacity of the particles off the target.
**Objective and principle of the method.** The objective of these experiments is to launch secondary nuclear reactions by a clever choice of complementary components of the target irradiated by the proton beam. The advantage of secondary reactions is the production of the elements that are needed to repeat the primary reactions and so to increase their yield. In the case of our experiments, the α-particles produced in the primary reactions react either with $^{14}$N or $^{10}$B to generate 'new' protons that have the right energy to react again with the $^{11}$B. The other secondary reactions are used as diagnostic tool. The observation of $^{18}$F unequivocally demonstrates that α-particles have reacted with $^{14}$N. This is a unique signature. In our experiments, the unique feature of the secondary reactions is that they use particles that are produced *in situ* in the target whereas the primary reactions are produced by the accelerated proton beam coming from the Aluminum foil. The observation of such reactions is of great interest for many applications where it can be cheaper to produce the reactant by interaction with a target rather than accelerating the ions.

**References**

1. Nevins, W. M. & Swain, R. The thermonuclear fusion rate coefficient for p-$^{11}$B reactions. *Nuclear Fusion* **40**, 865-872 (2000).

2. Atzeni, S. & Meyer-Ter-Vehn, J. The *Physics of Inertial Fusion*. Oxford Science, Publication (2004).

3. Hora, H. Laser-optical path to nuclear energy without radioactivity: Fusion of hydrogen-boron by nonlinear force driven plasma blocks. *Op. Comm.* **282**, 4124-4126 (2009).

4. Picciotto, A. et al. Boron-proton nuclear fusion enhancement induced in boron-doped silicon targets by low-contrast pulsed laser. *Physical Review X* **4**, 031030 (2014).

5. Strickland, D. & Mourou, G. Compression of amplified chirped optical pulses. *Opt. Commun.* **56**, 447-449 (1985).

6. Labaune, C. et al. Fusion reactions initiated by laser-accelerated particle beams in a laser-produced plasma. *Nat. Comm.* **4**, 2506-2512 (2013).

7. Wilks, S. C. et al. Energetic proton generation in ultra-intense laser-solid interactions. *Physics of Plasmas* **8**, 542-549 (2001).

8. Fews, A. P. et al. Plasma Ion Emission from High Intensity Picosecond Laser Pulse Interactions with Solid Targets. *Phys. Rev. Lett.* **73**, 1801-1804 (1994).

9. Moreau, D. C. Potentiality of the proton-boron fuel for controlled thermonuclear fusion. *Nuclear Fusion* **17**, 13-20 (1977).

10. Belyaev, V. S. et al. Observation of neutronless fusion reactions in picosecond laser plasmas. *Phys. Rev. E* **72**, 026406 (2005).

11. Ajzenberg-Selove, F. Energy levels of light nuclei A = 11-12. *Nuclear Physics A506*, 1-158 (1990).

12. Becker, H. W., Rolfs, C. & Trautvetter, H. P. Low-Energy Cross Sections for $^{11}$B(p,3α)*. Z. Phys. **327**, 341-355 (1987).

13. Dmitriev, V. F. α particle Spectrum in the Reaction p+$^{11}$B → α + 8Be* → 3α*. *Phys. Atom. Nucl.* **72**, 1165-1167 (2009).

14. Stave, S. et al. Understanding the $^{11}$B(p,α)αα reaction at the 0.675 MeV resonance. *Phys. Lett. B* **696**, 26-29 (2011).

15. Kafkarkou, A. et al. Study of the $^{10}$B(p,α) reaction between 2.1 and 6.0 MeV. *Nuclear Instruments and Methods in Physics Research B* **316**, 48-55 (2013).

16. Tonchev, A. P. et al. The $^{10}$B(p,γ)$^{11}$C reaction at astrophysically relevant energies. *Phys. Rev. C* **68**, 045803-1-12 (2003).

17. Van der Zwan, L. & Geiger, K. W. The $^{10}$B(α,n)$^{13}$N, $^{13}$N* cross section for a-energies from 1.0 to 5 MeV. *Nuclear Physics A* **216** 188-198 (1973).

18. Wilson, S. J., $^{10}$B(α,n)$^{13}$N and $^{10}$B(α,p)$^{13}$C for 2<E<8 MeV. *Phys. Rev. C* **11**, 1071-1082 (1975).

19. Bailey, S. et al. "Energy levels of $^{18}$F from 14N+α resonant reactions", *Phys. Rev. C* **90**, 024302-11 (2014).

20. Terwagne, G., Genard, G., Yedji, M. & Ross, G. G. Cross-section measurements of the $^{14}$N(α,p)$^{17}$O and $^{14}$N(α,α)$^{14}$N reactions between 3.5 and 6 MeV. *J. of Appl. Phys.* **104**, 084909-1-8 (2008).

21. Jacobs, W. W., Bodansky, D., Chamberlin, D. & Oberg, D. L. Production of Li and B in proton and alpha-particle reactions. *Phys. Rev. C* **9**, 2134-2143, (1974).

22. Kitwanga, S., Leleux, P., Lipnik, P. & Vanhorenbeeck, J. Production of $^{14,15}$O, $^{18}$F and $^{19}$Ne radioactive nuclei from (p,n) reaction up to 30 MeV. *Phys. Rev. C* **42**, 748-752 (1990).

23. Baccou, C. et al. CR-39 track detector calibration for H, He and C ions from 0.1-0.5 keV up to 5 MeV for laser-induced nuclear fusion product identification. *Review of Scientific Instruments*, **86**, 083307 (2015).

24. Baccou, C. et al. New scheme to produce aneutronic fusion reactions by laser-accelerated ions. *Laser and Particle Beams*, Cambridge University Press 0263-0346/15 (2015).

25. Kimura, S., Anzalone, A. & Bonasera, A. Comment on Observation of neutronless fusion reactions in picosecond laser plasmas. *Phys. Rev. E* **79**, 038401 (2009).

**Acknowledgements**

We acknowledge the support of the Elfie teams and discussions with G. Boutoux, J.-E. Ducret, S. Depierreux, F. Consoli and R. de Angelis. J.R. wishes to thank Ecole Polytechnique for the support of his several two-month sabbatical research visit and Christine Labaune for hospitality at the LULI laboratory where this work was carried out.

**Author Contributions**

C.B., C.N. and V.Y. carried out the experiments, installed the diagnostics and recorded the data. C.B., C.L. and J.R. analyzed the results. C.L. and J.R. conceived the idea of the scheme. C.L., C.B. and J.R. wrote this paper.

**Additional Information**

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Labaune, C. et al. Laser-initiated primary and secondary nuclear reactions in Boron-Nitride. *Sci. Rep.* **6**, 21202; doi: 10.1038/srep21202 (2016).

<figure><img src="image_8.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/