## Maximizing specific energy by breeding deuterium

## Justin Ball

Ecole Polytechnique F´d´rale de Lausanne (EPFL), Swiss Plasma Center (SPC), e e CH-1015 Lausanne, Switzerland

E-mail: Justin.Ball@epfl.ch

## Abstract.

Specific energy (i.e. energy per unit mass) is one of the most fundamental and consequential properties of a fuel source. In this work, a systematic study of measured fusion cross-sections is performed to determine which reactions are potentially feasible and identify the fuel cycle that maximizes specific energy. This reveals that, by using normal hydrogen to breed deuterium via neutron capture, the conventional catalyzed D-D fusion fuel cycle can attain a specific energy greater than anything else. Simply surrounding a catalyzed D-D reactor with water enables deuterium fuel, the dominant stockpile of energy on Earth, to produce as much as 65% more energy. Lastly, the impact on space propulsion is considered, revealing that an effective exhaust velocity exceeding that of deuterium-helium-3 is theoretically possible.

## 1. Introduction

The aim of this work is to investigate the basic question 'What fuel source can provide the most energy with the least mass?' For this question, we will look into the near future and assume that the current incremental and steady development of technology will continue. However, we will not postulate any breakthroughs or new physical processes because they cannot be predicted with any certainty. It is already well-known that the answer to our question must involve fusion fuels [1, 2], but nothing more precise has been rigorously established ‡ . A more precise answer may be useful in the future when the technology required to implement it becomes available. Specific energy (i.e. energy per unit mass) is such a basic physical property that it seems almost certain to have important applications. However, perhaps more importantly, the answer to such a fundamental question is also of intrinsic value, simply to satisfy our curiosity - to know.

- ‡ Matter-antimatter annihilation does have a higher specific energy than fusion, but it appears impractical for two reasons. First, antimatter does not occur naturally in our Solar System [3] and is inefficient to produce. It currently requires over 10 million units of input energy to create enough antimatter to generate 1 unit of energy. Second, it is difficult to store because it annihilates with normal matter. The most that has been confined by a device corresponds to just 300 Joules [2].

In this work, we will first introduce the conventional fusion fuel cycles, focusing on the catalyzed D-D cycle. Then, in sections 3 and 4 we will propose a simple extension to it using the neutron capture reaction by normal hydrogen and give an example of how it could be implemented. Next, in section 5 we will argue that this fuel cycle possesses the maximum specific energy of any cycle that appears even remotely possible. Lastly, section 6 will consider the fuel cycle's utility for spacecraft propulsion and section 7 will provide some concluding remarks.

## 2. Fusion fuel cycles

When fusion energy research started in the 1950s, a great deal of attention was given to fusion fuel cycles in order to determine which was optimal for use in future power plants [4, 5, 6, 7]. After a short period of hesitation [5, 8, 9], the deuterium-tritium (D-T) fusion reaction,

$$^ { 2 } H + ^ { 3 } H \to ^ { 4 } H e + ^ { 1 } n + 1 7. 6 \text{ MeV},$$

was identified as the most promising reaction [10, 11]. It is optimal because it has the largest fusion cross-section ( ∼ 5 barns) and this occurs at a relatively low reactant energy ( ∼ 60 keV). However, since tritium has a half-life of just 12.3 years, it is not abundant in nature. Because of this, it is generally envisioned to be produced within the power plant, rather than being an input. This can be done by surrounding the reactor with a 'blanket' that includes lithium, in order to capture the neutron produced by D-T fusion to 'breed' tritium, e.g.

$$^ { 6 } \text{Li} + ^ { 1 } n \rightarrow ^ { 3 } \text{H} + ^ { 4 } \text{He} + 4. 8 \text{ MeV}.$$

Given a small amount of tritium for start-up, this 'breeding' reaction closes the overall D-T fuel cycle,

$$^ { 2 } \mathrm H + ^ { 6 } \mathrm L i \rightarrow 2 \ ^ { 4 } \mathrm H e + 2 2. 4 \ \mathrm M e V,$$

because both deuterium and lithium-6 are abundant in nature.

The deuterium-helium-3 reaction (D-He-3) is generally considered to be the next most promising reaction. It has a maximum cross-section of ∼ 0.8 barns at ∼ 200 keV and is given by

$$^ { 2 } H + ^ { 3 } H e \to ^ { 4 } H e + ^ { 1 } H + 1 8. 4 \text{ MeV}.$$

It is attractive, but helium-3 is fairly rare on Earth [12] and the reaction does not produce a neutron with which to breed it. Thus, this reaction requires helium-3 to be found elsewhere (e.g. the moon [12, 13]) or to be produced using other nuclear reactions [12, 14].

The next best fusion reaction occurs between two deuterium nuclei (D-D). While the maximum cross-section is just ∼ 0.2 barns, which occurs at ∼ 1000 keV, this peak is

very broad. This means that significant fusion power can be generated at temperatures much lower than 1000 keV. Interestingly, the D-D reaction has two branches that occur with almost equal probability,

$$2 ^ { \ 2 } H \stackrel { \succ } { \succ } ^ { 3 } H + ^ { 1 } H + 4. 0 \text{ MeV} \\ \succ ^ { 3 } \text{He} + ^ { 1 } n + 3. 3 \text{ MeV}.$$

Note that, while the previous two reactions produce the exceptionally stable helium-4 nucleus, this reaction produces tritium and helium-3. This allows for the possibility of the 'catalyzed' D-D fuel cycle [5, 15, 16], where the tritium and helium-3 produced by D-D fusion are fused with deuterium in subsequent reactions. This can be done in the same device or the products can be extracted and fused in a separate devices. Regardless, the overall catalyzed D-D fuel cycle is given by the sum of equations (1), (4), and (5):

$$3 \, ^ { 2 } H \to ^ { 4 } H e + ^ { 1 } H + ^ { 1 } n + 2 1. 6 \, \text{MeV}.$$

Such a fuel cycle performs significantly better than D-D reactions on their own. In fact, the temperature required to achieve a self-sustaining fusion burn with catalyzed D-D is lower than for D-He-3 and the fusion power density of the plasma is higher (though the performance is still significantly worse than D-T) [15, 16].

These are the only reactions that currently appear to have practical application [17, 18]. All other fusion reactions are more difficult. They have lower fusion crosssections that occur at higher temperatures and/or involve reactants with higher electric charge. Because of this, the power lost by electron bremsstrahlung looks to exceed the fusion power for nearly all other fusion reactions [17]. This is a fairly fundamental constraint because bremsstrahlung radiation is emitted whenever charged particles collide, which is intrinsic to the process of fusion. Moreover, bremsstrahlung is generally in the ultraviolet/X-ray range of frequencies, so it looks difficult to reflect, reabsorb in the plasma, or convert to electricity at high efficiency [17]. Unless this problem can be addressed (e.g. maintaining cold electrons [19], removing electrons [20], non-equilibrium plasmas [21, 22]), it does not appear possible to produce net electricity using nearly all other fusion reactions. At the present time, even the D-T fusion reaction proves very challenging, with current fusion devices still somewhat short of net energy production [23].

## 3. Breeding deuterium

This work proposes a simple extension to the catalyzed D-D fuel cycle: to use normal hydrogen (i.e. protium) to breed deuterium with the neutron capture reaction,

$$^ { 1 } H + ^ { 1 } n \rightarrow ^ { 2 } H + 2. 2 \text{ MeV}.$$

This reaction has a cross-section of ∼ 0 3 . barns for thermal neutrons [24]. Though the cross-section may seem small compared to D-T fusion, this reaction is actually

Table 1: The specific energy of some of the easiest fusion fuel cycles.

| Fuel cycle        | Specific energy   | Fuel cycle   | Specific energy   |
|-------------------|-------------------|--------------|-------------------|
| Catalyzed D-D+D   | 6.0 MeV/AMU       | D-Li-6       | 2.8 MeV/AMU       |
| D-He-3            | 3.7 MeV/AMU       | He-3-He-3    | 2.2 MeV/AMU       |
| Catalyzed D-D     | 3.6 MeV/AMU       | T-T          | 1.8 MeV/AMU       |
| D-T               | 3.5 MeV/AMU       | D-D          | 0.9 MeV/AMU       |
| D-T w/ T breeding | 2.8 MeV/AMU       | p-B-11       | 0.7 MeV/AMU       |

much easier to accomplish. Because it occurs at room temperature and neutrons don't interact much with electrons, you can place a liquid (or even a solid) containing normal hydrogen in the path of neutrons. Liquids have a much greater number density than the plasmas typical of fusion, so the fact that the individual nuclei have somewhat smaller cross-sections can be overcome. The neutrons will enter the material, thermalize due to scattering collisions, and be captured.

Using the neutron capture reaction enables the fuel cycle

$$2 ^ { \ 2 } H \to ^ { 4 } \text{He} + 2 3. 8 \text{ MeV}.$$

This set of reactions will be referred to as the 'catalyzed D-D+D' fuel cycle. While this idea seems fairly obvious from equation (6), the practicality and advantage of deuterium breeding does not appear to have been noted. For example, references [1, 2, 25] all conclude that D-He-3 has the highest specific energy of any fuel. However, the specific energy of the catalyzed D-D+D cycle is more than 60% larger (see table 1) § . Moreover, the plasma performance would be identical to the normal catalyzed D-D (which is seen as one of the most viable fusion fuel cycle after D-T) because the only changes we have introduced occur outside the plasma. Additionally, capturing the fusion neutrons with normal hydrogen prevents them from activating other surrounding material, thus minimizing the production of radioactive waste.

Assuming the ability to achieve a conventional catalyzed D-D fuel cycle, the addition of deuterium breeding looks to be attractive and easily achievable. It can increase the specific energy of deuterium fuel by as much as 65% relative to a catalyzed D-D fuel cycle without deuterium breeding. Put another way, any given deuterium fuel supply could be used to generate as much as 65% more energy, without any additional inputs. This could be useful for a number of future applications. First, deuterium is the third most abundant isotope in the universe [26] and one in every 6420 hydrogen atoms on Earth is deuterium. Moreover, the two most abundant isotopes, hydrogen-1 and helium-4, appear impossible to fuse without using other less abundant isotopes.

§ Since D-He-3 releases all of its energy in the form of charged particles that have the potential to be converted to electricity at a high efficiency, one might think that D-He-3 could have a higher specific energy in practice. However, catalyzed D-D+D releases 3.4 MeV/AMU in charged particles alone, already almost equal to that of D-He-3. If the energy carried by neutrons can be made useful at even a low efficiency, catalyzed D-D+D will be more energy dense.

Helium-3 can be fused and is almost as abundant as deuterium in the universe, but not in our neighborhood. It is nearly non-existent on Earth and, even when totaled over all the planets in our Solar System, contains just 1/50 of the energy stored by deuterium in the Earth's oceans [2]. Thus, deuterium appears to be our dominant usable stockpile of energy. By partially breeding it, our supply can be used most efficiently.

Second, an ultra-high specific energy fuel can be advantageous in any situation that requires fuel to be transported. A clear example is industrial activities at stationary, remote outposts in outer space. For example, asteroid mining, transmission stations, or large construction projects could be best fueled by catalyzed D-D+D. Such activities would require a lot of energy, but be insensitive to the mass and volume of the power plant (especially if the raw material for the plant could be found on-site). Another particularly important application could be space propulsion, which we will consider in a later section. Of course, all these applications are sufficiently far in the future that we can only speculate. Regardless, specific energy is such a fundamental property that the fuel that maximizes it seems likely to have important applications.

## 4. An illustrative example

A simple method to implement the catalyzed D-D+D fuel cycle is to surround the D-D and D-T reactions with a blanket of normal light water. This could be pressurized, similar to a pressurized water fission reactor (PWR), or allowed to boil, similar to a boiling water fission reactor (BWR). Figure 1 illustrates such a configuration using an ITER-like tokamak as an example confinement scheme [27], but it could just as well be an inertial confinement scheme or something more exotic. The only constraint is that there be space for the deuterium breeding blanket, which figure 2 shows would have a similar thickness to the tritium breeding blankets in current D-T power plant designs [28, 29]. Using light water makes for an attractive design because it is the most common coolant in power plants and has already been used extensively in the neutron environment of light water fission reactors. Hence, its neutronic properties are well studied. Nuclear fission reactors have been observed to generate kilograms of deuterium per year. The only significant difference for a fusion reactor is the energy of the neutrons -fission neutrons are born with ∼ 1 MeV of energy while D-T neutrons have 14.1 MeV and D-D neutrons have 2.5 MeV. Nevertheless, the cross-sections of the neutron-induced reactions of hydrogen-1 and oxygen-16 are well characterized at these energies.

To demonstrate the feasibility of such a design, neutronics simulations were performed using MCNP [30]. An axial source of D-D and D-T neutrons (weighted equally) was entirely surrounded by thick cylindrical blanket of liquid light water at atmospheric pressure and temperature. The deuterium breeding ratio (i.e. the number of deuterium nuclei produced per fusion neutron) is shown in figure 2 for different blanket thicknesses. Note that oxygen-16 reacts with a very small fraction of the neutrons because it a double magic isotope, so it is exceptionally stable. Lastly, in reality a design would likely have solid metal plasma-facing components in front of the

Figure 1: A cartoon of a power plant that implements the catalyzed D-D+D fuel cycle.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00834v1-figure-1e729b5d-83db-489a-8ee4-c2b87088e004.jpg)

blanket. This structure would absorb some of the neutrons before they can reach the hydrogen. While this isn't ideal, unlike tritium breeding for the D-T fuel cycle, it is perfectly acceptable for the DBR to be less than one. If some of the neutrons are lost then the fuel cycle will require more deuterium as input, but the overall fuel cycle will still function. Note that using material to multiply neutrons (e.g. beryllium or lead) could be used, but these are substantially less abundant than deuterium and would only decrease the specific energy of the fuel cycle as a whole because of their large atomic mass.

In the design shown in figure 1, the liquid light water serves as both a coolant and a breeding blanket. The fusion reaction chamber is submerged in water, which slowly circulates past. This is similar to breeding blankets that have already been proposed for D-T power plants [29]. As it passes by the tank, the water intercepts the neutrons, causing it to heat up and produce deuterium. The heated water can then be used to power a steam cycle and generate electricity. Afterwards, the deuterium could be

Figure 2: The deuterium breeding ratio as a function of the blanket thickness, where the horizontal dashed line gives the value in the limit of a very thick blanket.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00834v1-figure-366880e2-ff14-435b-a4e1-227c4775db08.jpg)

separated from the normal water by the Girdler Sulfide process, which is currently used to extract deuterium from seawater. Periodically, normal hydrogen collected from the plasma could be combined with the oxygen leftover from the deuterium extraction. This would be added to the coolant. Unlike in the D-T fuel cycle, the speed of such a deuterium breeding cycle is not particularly important. In D-T power plants, the time required to breed, extract, and fuse the tritium determines the total power plant tritium inventory, which must be minimized for safety reasons. Deuterium, on the other hand, can be allow to accumulate for a long time and does not disappear due to radioactive decay.

Finally, it should be mentioned that, for any particular application, there are many competing design considerations and the specific energy of the fuel is just one. That is to say that the high specific energy of the catalyzed D-D+D fuel cycle is only useful to the extent that the specific energy of the fuel actually matters. For many systems, specific energy is only a secondary consideration compared to factors like neutron activation, plasma power density, and the required fusion technology. Indeed for many applications, any fusion fuel already has such high specific energy that it can be negligible compared to other necessary components. For example, in present-day fusion power plant designs [29], a water blanket as depicted above could easily have a mass 1000 times greater than the mass of deuterium that it would breed each year. Conceptual designs of DHe-3 spacecraft thrusters for travel within the Solar System [18, 25], which are designed to minimize mass, are still over 100 times more massive than the fuel they use. Of course, fuel mass would be more significant in devices that are larger, have higher power densities, and longer lifespans. Such devices would hopefully be developed in the future and would likely be required for interstellar space travel. Moreover, there may be situations in which normal hydrogen and/or structural materials are abundant, but

fusion fuel is difficult to access.

## 5. Maximizing specific energy

To argue that the catalyzed D-D+D fuel cycle (i.e. equation (8)) maximizes specific energy, we first note that nuclear forces are more than a million times stronger than the other two fundamental forces. The gravitational potential energy of two nuclei, even if they are just a femtometer apart, is less than 10 -25 eV. Electromagnetism is responsible for chemical bonds between atoms, the strongest of which are roughly 10 eV. Nuclear reactions, on the other hand, typically release millions of eV. Hence, the reactions with the highest specific energy seem almost certain to be nuclear.

The energetics of nuclear reactions are well summarized by the mass per nucleon curve shown in figure 3. Due to E = mc 2 , converting reactants with high mass nucleons to products with low mass nucleons will release energy. More precisely, the specific energy of a nuclear reaction is well approximated by the average mass of a nucleon in the reactants minus the average mass of a nucleon in the products (times the speed of light squared). Strictly speaking, this actually gives the energy released per nucleon , rather than the specific energy (i.e. the energy released per mass ). However, figure 3 shows that these two quantities are identical to within 1%. Regardless, we would like to use reactants near the top of this plot to make products near the bottom.

From figure 3, we see that converting two deuterons into one helium-4 releases almost 6 MeV/AMU, which spans over 70% of what is theoretically attainable (i.e. the reaction 56 1 H → 56 Fe+8.4 MeV/AMU). Intuitively, this seems hard to beat, especially considering how difficult fusion reactions are. Even if the only product of a fuel cycle is the most energetically favorable isotope (i.e. iron-56), the reactants still must have an average mass per nucleons greater than 936.4 MeV/AMU to beat catalyzed D-D+D. By this logic, we can rule out all fission processes. The input nucleons are too light (see figure 3).

Thus, we can restrict our considerations to fusion fuel cycles (supplemented by other nuclear reactions like neutron capture and radioactive decay). At this point, things become a bit more subjective because we must decide which fusion reactions are practical to achieve. Converting 56 protons into iron-56 has the highest specific energy, but it has only been accomplished over millions of years using a star and looks absolutely impossible to do on any sort of human scale. However, there could be other fairly viable fusion reactions that do better than catalyzed D-D+D. This is a messy empirical question, so we will perform a systematic study of 'the most comprehensive compilation of experimental nuclear reaction data' - the EXFOR database [31].

To estimate the difficulty of various fusion reactions we will use a variant of a common fusion metric: the triple product required for ignition [32]. The triple product nT τ i E is the product of the electron density, ion temperature, and energy confinement time. In order for a fuel to ignite, it must achieve a minimum triple product given by

Figure 3: The rest mass per nucleon for all isotopes with half-lives longer than a year (diamonds), a second (crosses), and a microsecond (dots). The value of iron-56 plus the specific energy of the catalyzed D-D+D cycle is indicated by the horizontal dashed line.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00834v1-figure-b7841253-7daf-4f6f-899d-d5fa7623db57.jpg)

[33]:

$$n T _ { i } \tau _ { E } = 3 \left ( 2 Z _ { 1 } Z _ { 2 } + Z _ { 1 } + Z _ { 2 } \right ) \frac { T _ { i } ^ { 2 } } { \langle \sigma v _ { i } \rangle } \frac { 1 } { f _ { c } E _ { f u s } }.$$

The lower this value is, the easier the fusion reaction is. Here Z 1 and Z 2 are the charge numbers of the two reactants, σ is the fusion reaction cross-section, v i is the relative speed between the two reactants, 〈 . . . 〉 is an average over the ion velocity distributions, E fus is the energy released by the fusion reaction, and f c is the fraction of it that is carried by charged particles. To calculate this we have assumed that both of the ion densities are n i = n/ (2 Z i ) because this choice maximizes the fusion power density. Note that if the two reactants are identical (e.g. D-D fusion), the prefactor becomes 3 Z 1 ( Z 1 +1) and the ion density is required to be n 1 = n/Z 1 by the quasineutrality condition.

For our purposes, we will make a few adjustments to equation (9). First, we will multiply by f c because we are primarily concerned about the total energy released, rather than details about what energy is deposited in the plasma and how efficiently power can be recirculated. Second, because the EXFOR database contains the results of many independent experimental measurements, the cross-section values are not evenly distributed in energy. Many unusual reactions have just a few points. Fortunately,

we only require a rough estimate of the difficulty of a given reaction. Thus, we will calculate the 'point' triple product by assuming that the ion distribution function is a delta function at a given center-of-mass energy E i . With these assumptions, equation (9) becomes

$$\ n E _ { i } \tau _ { E } f _ { c } = 3 \, ( 2 Z _ { 1 } Z _ { 2 } + Z _ { 1 } + Z _ { 2 } ) \, \frac { E _ { i } ^ { 2 } } { \sigma \, ( E _ { i } ) \, \sqrt { 2 E _ { i } / m _ { i } } } \frac { 1 } { E _ { f u s } },$$

where m i is reduced ion mass of the reactants. This metric applies equally well to a dataset with one point or a thousand. While this is a simple and useful metric to give a rough estimate, it is not perfect. It will tend to overrate fusion reactions that have sharp, high resonances in their cross-sections and underrate reactions with broad, low cross-sections (e.g. D-D). Additionally, it doesn't make much sense for endothermic reactions because a larger magnitude of E fus is a bad thing. Thus, we will compute nE τ i E f E c 2 fus for endothermic reactions.

Figure 4 shows these values for the complete contents of the EXFOR database. The column and row labels indicate the two reactants of the fusion reaction. The table boundaries were chosen such that it includes all fusion reaction with nE τ i E f c &lt; 10 26 keVs/m . 3 This number was chosen as the cutoff for potentially feasible reactions because it is more than a factor of 10 above p-B-11, the hardest reaction that some present-day researchers are pursuing [19, 21, 22]. Isotopes with no data were omitted. The first number of the first row within each cell (as well as the shading of the cell) indicates the minimum value of the point triple product across all energy values and datasets in the database. Specifically, the number is log 10 ( nE τ i E f c ), where nE τ i E f c is expressed in units of keV-s/m . 3 The number that follows the '@' sign is the ion energy (in MeV) that this minimum value occurs at. Easier fusion reactions are indicated by both of these numbers being smaller. The second line gives the products that are generated by the reaction. Any product that has a half-life of less than 1 second has been replaced by its daughter products, unless it appears as a reactant on the chart. The first number in the third line indicates the fusion energy (in MeV) released by the reaction. The number after the comma is the ratio of the fusion power density to the bremsstrahlung power density radiated by electrons. This is calculated with the formulas [17, 34]

$$p _ { \text{fusion} } = n _ { 1 } n _ { 2 } \sigma \left ( E _ { i } \right ) \sqrt { 2 E _ { i } / m _ { i } } E _ { f u s },$$

$$p _ { \text{fusion} } & = n _ { 1 } n _ { 2 } \vartheta \left ( \mathcal { E } _ { i } \right ) \vee \mathcal { E } _ { i } / n _ { i } \mathcal { E } _ { \text{fus} }, & \left \{ 1 1 \right ) \\ p _ { \text{brem} } & = \frac { 1 6 0 } { 3 } \left ( \frac { e ^ { 2 } } { 4 \pi \epsilon _ { 0 } } \right ) ^ { 3 } \frac { n ^ { 2 } \sqrt { E _ { i } } } { \left ( m _ { e } c ^ { 2 } \right ) ^ { 3 / 2 } \hbar { \, } } \\ & \times \left ( \frac { \sum _ { i } Z _ { i } ^ { 2 } n _ { i } } { n } \left ( 1 + 0. 7 9 \frac { E _ { i } } { m _ { e } c ^ { 2 } } + 1. 9 \left ( \frac { E _ { i } } { m _ { e } c ^ { 2 } } \right ) ^ { 2 } \right ) + \frac { 3 } { \sqrt { 2 } } \frac { E _ { i } } { m _ { e } c ^ { 2 } } \right ).$$

$$& ) \sqrt { 2 E _ { i } / m _ { i } } E _ { f u s }, & & ( 1 1 ) \\ & _ { \infty } ^ { - } \Big ) ^ { 3 } \frac { n ^ { 2 } \sqrt { E _ { i } } } { ( m. e ^ { 2 } ) ^ { 3 / 2 } \, \hbar { \, } } & & ( 1 2 )$$

Here e is the electric charge of a proton, glyph[epsilon1] 0 is the vacuum permittivity, m c e 2 is the electron rest energy, and glyph[planckover2pi1] is the reduced Planck constant. Note that for simplicity, the electron temperature has been assumed to be equal to the ion energy . The ion densities n 1 and n 2 are again determined in order to maximize the fusion power density, while

1

Figure 4: A summary of all measured fusion reactions in the EXFOR database [31] between light nuclei. The stable (black) and unstable (gray) reactants are given by the column and row labels, while the legend in the upper-right corner explains the contents of each reaction cell. Negative values indicate endothermic reactions (see text).

|                                                                                  | H                                                  |                                                                             |                                              | Minimum log 10                                                                       | ( nE i τ E c keV-s/m 3                           | f (                                                  | ))                                     |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
|----------------------------------------------------------------------------------|----------------------------------------------------|-----------------------------------------------------------------------------|----------------------------------------------|--------------------------------------------------------------------------------------|--------------------------------------------------|------------------------------------------------------|----------------------------------------|----------------------------------------------------|----|----|---------|----|----|----|----|------|----|----|------|------|------------|----|----|----|------|----|----|----|
| 1 H No Data                                                                      | 2 H                                                |                                                                             |                                              | (also indicated by shading) (MeV) E i                                                |                                                  |                                                      |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 2 H 28.6 @ 0.036 3 He                                                            |                                                    | 23.7 @ 0.035 3 H + 1 H 3                                                    | 23.9 @ 0.059 He + 1 n                        | 3 H 12 yr                                                                            | 21.1 @ 0.048 4 He + 1 n Products                 |                                                      |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 3 5.5,0.000 27.5 @ 6.3 4 He                                                      | 4.0,0.33 21.1                                      | @ 0.048                                                                     | 3.3,0.39                                     | 23.3 @ 0.03                                                                          |                                                  | 17.6,170                                             |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| H 19.8,0.000 12 yr                                                               | 4                                                  | He + 1 n 17.6,170                                                           |                                              | 4 He + 2 1 n 11.3,0.59                                                               | 3 He 25.6 @ 0.75                                 | (MeV) P f us                                         |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| He No Data                                                                       |                                                    | 23.1 @ 0.21 4 He + 1 H 18.4,6.2                                             |                                              | 24.7 @ 0.22 4 He + 1 H + 1 n 12.1,0.14 30.1 @ 0.21 25.1 @ 0.23 4 He + 2 H 14.3,0.060 | 4 He + 2 1 H 12.9,0.023                          | 4 He - 35.7 @ 22                                     | 6 Li                                   | P f us / P brem                                    |    |    |         |    |    |    |    |      |    |    | 30.8 |      |            |    | @  |    | 0.57 |    |    |    |
| He 6 Li - 36.3 @ 44 3 He + 2 H - 18.4, - 0.000 25.1 @ 0.34 4 He + 3 He 4.0,0.061 | 24.3 @ 0.43 2 4 He 22.4,0.46 25.1 7 Li             | 31.7 @ 0.63 6 Li 1.5,0.000 @ 0.43 + 1 H 5.0,0.068 25.2 @ 7 Be + 1 3.4,0.056 | 0.6 n 25.4 @ 0.36                            | 7 Li 2.5,0.000 25.0 @ 0.51                                                           | 1.6,0.000 27.2 @ 1.1 7 Be + 2 H 0.1,0.000        | 7 Be + 1 n - 19.0, - 0.000 31.4 @ 2.4 10 B           | 4.5,0.000 26.7 @ 1.8 7 Be + 4 He + 1 n | 7 Li                                               |    |    |         |    |    | +  | 1  |      | 7  | 4  |      |      |            | He |    |    |      | 2  | n  | Be |
| 7 Li 25.6 @ 2.6 2 4 He 17.3,0.022 31.6 @ 0.64 8 B 0.1,0.000 50 day               | No                                                 | 24.2 @ 0.62 2 4 He + 1 n 15.1,0.60                                          | 4 He + 3 H + 1 H 2.6,0.031                   | 16.1,0.097 24.2 @ 0.39 9 Be + 1 n 10.4,0.52                                          | - 33.1 @ 2.8 7 Be + 3 H - 0.9, - 0.000           | - 33.8 @ 4.5 10 B + 1 n - 2.8, - 0.001               | 1.9,0.001 No Data                      | No Data                                            |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| Be 8 1 sec                                                                       | 24.8 @                                             | Data 0.3                                                                    |                                              | No Data                                                                              | No Data                                          | No Data                                              | No Data                                | No Data                                            |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| Li 8 B No Data No Data 1 sec                                                     |                                                    | 7 Li + 3 H 4.2,0.12 No Data                                                 |                                              | No Data                                                                              | No Data No Data                                  | 25.2 @ 0.66 11 B + 1 n 6.6,0.045 27.2 @ 1.8 11 1     | No Data No Data                        | No Data No Data                                    |    |    |         | H  | C  |    |    |      |    | No |      | Data |            |    |    | +  |      |    |    |    |
| Be 10 24.7 @ 0.3 6 Li + 4 He 2.1,0.12 - 34.4 @ 1 25.1 @ 2 4 He 0.7,0.042         | 25.4 @ 0.41 7 Li + 4 He 7.2,0.030 0.27 + 2 H       | 25.5 @ 1.5 10 B + 1 n 4.4,0.022 25.8 @ 2 4 He 4.7,0.014                     | 0.87 + 3 H 25.9 @ 0.86 10 Be + 1 H 4.6,0.011 | 25.0 @ 1.1 11 B + 1 n 9.6,0.084                                                      | 26.3 @ 1.5 11 C + 1 n 7.6,0.003                  | 7.4,0.000 25.9 @ 1.3 12 C + 1 n 5.7,0.008            | No Data No Data                        | 12 C + 4 He 23.8,0.001 No Data                     |    |    |         |    |    |    |    |      |    |    |      |      | 26.3 @ 8.9 |    |    |    |      |    |    |    |
| Be 10 B Be + He 10 6 yr                                                          | 10 B + 1 n - 0.2, - 0.000 26.2 @ 0.71 7 4          | No Data 25.6 @ 1.3                                                          |                                              | No Data 26.2 @ 2.3 12 C + 1 H                                                        | No Data 27.0 @ 3.6 11 C + 2 H                    | No Data 27.7 @ 3.4 13 N + 1 n                        | No Data                                | No Data                                            |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 11 B 12 C 1.1,0.004 24.8 @ 0.15 3 4 He 8.7,0.048 28.6 @ 0.39                     | 6.5,0.016 25.0 @ 1.1 12 C + 1 H 14.5,0.061         | 11 C + 1 n 29.1 @ 54                                                        |                                              | 19.7,0.004 26.9 @ 2 13 C + 1 H 13.2,0.001                                            | 3.2,0.000 - 34.9 @ 7.1 11 C + 3 H - 2.0, - 0.000 | 1.1,0.000 27.5 @ 3.2 14 C + 1 H 0.8,0.000 33.0 @ 2.4 | No Data                                | 26.8 @ 3.5 13 C + 4 He + 1 n 10.0,0.001 27.1 @ 4.8 |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 13 C 14 C 13 N                                                                   | 25.5 @ 1.1 14 N + 1 n 5.3,0.020 25.3 @ 1.1         | 13 C + 1 H 2.7,0.000                                                        |                                              | 26.9 @ 0.89 14 C + 1 H 4.6,0.001                                                     | 27.0 @ 2.2 11 C + 4 He 1.9,0.000                 | 16 O 7.2,0.000                                       | No Data                                | 15 N + 4 He 12.4,0.000 26.7 @ 4.5                  |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 31.6 @ 0.12 14 N 7.6,0.000 - 31.9 @ 1.1                                          | 1.9,0.000                                          |                                                                             |                                              | No Data                                                                              |                                                  | 26.8 @ 2.1 16 O + 1 n 2.2,0.001                      | No Data                                | 15 N + 4 He + 1 n 7.4,0.000                        |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 14 N + 1 n                                                                       |                                                    | 15 1                                                                        |                                              |                                                                                      | No Data                                          |                                                      |                                        | No                                                 |    |    |         |    |    |    |    | Data |    | No |      | Data |            |    |    |    |      |    |    |    |
| 14 - 0.6, - 28.2 @ 0.26 15 O 7.3,0.000 10 4 yr                                   | 0.003 0.94                                         | N + n 8.0,0.032 26.2 @ 2.4 15 O + 1 n 5.1,0.003                             |                                              | No Data                                                                              | No Data 27.8 @ 5.3 13 N + 4 He                   | No Data - 35.0 @ 7.2 17 F + 1 n                      | No Data                                | No Data                                            |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| N 25.3 @                                                                         |                                                    | 26.8 @ 1.2                                                                  |                                              |                                                                                      | 10.0,0.000                                       | - 4.7, - 0.000 30.9 @ 1.3                            |                                        | No Data                                            |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 15 N                                                                             | 12 C + 4 He 5.0,0.026 29.9 @ 0.00016               | 16 N + 1 H 0.3,0.001                                                        |                                              | No Data                                                                              | No Data                                          | 19 F 4.0,0.000                                       | No Data 27.5 @ 4.7                     | No Data 27.8 @ 6                                   |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 16 O 17 O                                                                        | 17 F @ 0.64 18 F                                   | 26.4 @ 14 N + 4 He 3.1,0.002                                                |                                              | 27.1 @ 1.9 18 F + 1 n 1.3,0.000 27.4 2.3                                             | 27.3 @ 5.3 18 F + 1 H 2.0,0.000                  | 33.5 @ 1.9 20 Ne 4.7,0.000 28.0 @ 3.3                | 18 F + 4 He 6.1,0.000 No Data          | 19 F + 4 He 9.2,0.000 No Data                      |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 28.3 5.6,0.000                                                                   | 0.6,0.000                                          | No Data                                                                     | 0.93                                         | @                                                                                    |                                                  |                                                      |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
|                                                                                  |                                                    |                                                                             |                                              | 19 O + 1 H 3.5,0.000                                                                 | No Data                                          | 20 Ne + 1 n 0.6,0.000                                |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| 18 O 18 2 hr                                                                     | 25.4 @ 0.59 15 N + 4 He 4.0,0.018                  | 26.5 @ 1.2 16 N + 4 He 4.2,0.002                                            |                                              | No Data                                                                              | No Data                                          | 29.7 @ 2.4 22 Ne 9.7,0.000                           | No Data                                | No Data                                            |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| F 19 25.0 @ 16 4                                                                 | 25.3 @ 0.66 15 O + 4 He 2.9,0.020 0.83 25.1 @ 0.82 | 26.9 @ 1.6 20 F + 1 H                                                       |                                              | No Data                                                                              | 27.6 @ 2.8 18 F + 4 He                           | No Data                                              | 29.1 @ 4.2 24 Na + 1 H                 | No Data                                            |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| O + He 8.1,0.041 20 Ne                                                           | 12.8,0.031                                         |                                                                             |                                              | No Data                                                                              | 10.1,0.000                                       | 32.8 @ 0.73 22 Ne + 1 H                              |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
|                                                                                  |                                                    | 4.4,0.000                                                                   |                                              |                                                                                      |                                                  | 1.7,0.000                                            | 13.7,0.000                             |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
| F                                                                                |                                                    |                                                                             |                                              |                                                                                      |                                                  |                                                      |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
|                                                                                  |                                                    |                                                                             |                                              |                                                                                      |                                                  |                                                      |                                        |                                                    |    | No | Data    |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
|                                                                                  |                                                    |                                                                             |                                              |                                                                                      |                                                  |                                                      | Data                                   |                                                    | No |    | No Data |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |
|                                                                                  |                                                    | No Data                                                                     |                                              |                                                                                      |                                                  |                                                      |                                        |                                                    |    |    |         |    |    |    |    |      |    |    |      |      |            |    |    |    |      |    |    |    |

Figure 5: The minimum value of the point triple product for the data of figure 4.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00834v1-figure-521e4c39-b159-40c8-966e-bd5cfb24411f.jpg)

respecting quasineutrality. If there is only one ion species present the factor of n n 1 2 in equation (11) becomes n / 2 1 2 and the bremsstrahlung formula remains unchanged. Note that easier fusion reactions have larger numbers in the third line.

Lastly, some cells have more than one entry. This indicates that the reaction has multiple branches that are exothermic and have nE τ i E f c &lt; 10 26 keV-s/m . 3 If the reaction has several equally probable branches (e.g. D-D), then combining them increases the cross-section by a factor of N branch , the number of branches. To account for this in the logarithm of the point triple product, you should subtract log 10 ( N branch ). Additionally, it is appropriate to sum the different fusion-to-bremsstrahlung power ratios as the fusion power will increase due to the multiple reactions, but the power radiated by the electrons will not change.

The most eye-catching feature of figure 4 is the number of reactions that have no data. In some cases, this is because the cross-section is too small to be experimentally measured (e.g. p-p, p-He-3 [35]) or one of the reactants has a short half-life (e.g. Li-8, B-8). However, all these unknowns should be kept in mind as they make the conclusions of this work somewhat less certain. Other causes for uncertainty include the threshold that has been set at nE τ i E f c &lt; 10 26 keV-s/m . If we get better than this at fusion (or 3 can efficiently accomplish endothermic reactions), there could be useful reactions that don't appear in the figure. Additionally, some of the point triple product numbers can be changed somewhat. This is because we have made certain assumptions about what unstable products we retain or let decay (as described above). These choices affect the energetics of the reaction and can alter the triple product. The effect of this has been checked by repeating the analysis while retaining all products with a half-life greater than 1 microsecond. The conclusions that will be presented were not affected.

Digging deeper into figure 4, we see some expected results. D-T is the easiest

reaction to achieve ignition by a substantial margin and D-He-3 is next. D-D appears to be harder than T-T, but when you include both branches of D-D in the same point triple product by subtracting 0 3, you see that they are roughly even. . Moreover, T-T fusion releases a much higher fraction of its energy in the form of neutrons, which is harder to utilize. All other reactions are significantly more difficult, as is shown more clearly in figure 5. Though other reactions that have been seriously considered (e.g. D-Li-6, He-3He-3, p-Li-6, p-Be-9, p-B-11) have reasonably low point triple products compared to the rest. This gives confidence that it is a reasonable metric for comparison. Additionally, we see that the ion energies at the minimum point triple product have the right relative values to one another. However, they are systematically higher than the temperatures that detailed analyses indicate is optimal. This is understandable because a Maxwellian velocity distribution, unlike the delta function we assumed, has a long tail at high energy that has an outsized contribution to fusion energy production. Moreover, a Maxwellian distribution with an average particle energy of E i has a temperature of T i = 2 E / i 3. As a consequence, the fusion to bremsstrahlung power density ratios are lower than more realistic treatments [17]. However, their values relative to one another are still illuminating.

With the help of figures 3 and 4, we can now assess if it is possible to beat the specific energy of the catalyzed D-D+D fuel cycle. To do so, we have two options ‖ :

- (1) use inputs with an average mass per nucleon that is greater than hydrogen-2 and/or
- (2) make products with an average mass per nucleon that is less than helium-4.

Additionally, we can divide the methods to accomplish either of these into two strategies:

- (A) directly fuse input nuclei in a chain to produce larger nuclei or
- (B) use nuclei as catalysts to enable a catalytic fusion cycle.

We will consider each of the four combinations in turn.

First, we will look at option 1A: directly fusing inputs with an average mass per nucleon that is greater than hydrogen-2. Consulting figure 3, we see that the only isotope that can enable this is hydrogen-1. The free neutron is not naturally occurring as it has a half-life of just 10 minutes, so it cannot be used as an input. Stars directly fuse two hydrogen-1 nuclei into hydrogen-2 as part of the proton-proton chain, but this reaction has a cross-section so small that it has never been experimentally measured. More broadly, figure 4 shows that it not feasible to fuse hydrogen-1 with anything smaller than lithium-6. This appears to rule out option 1A because the combination of hydrogen-1 and lithium-6 has an average mass per nucleon of 934.5 MeV/c , considerably smaller 2 than that of hydrogen-2. Strictly speaking, it could be possible to use a single large nucleus in the first reaction and then repeatedly fuse the products with enough small nuclei for the average mass per nucleon of the reactants become large enough. However,

‖ We will not consider muon-catalyzed fusion because creating muons requires a lot of energy and they live for just 2.2 microseconds. We will also ignore the possibility of three-body fusion (e.g. the triple-alpha process) due to its improbability.

this does not look possible. Even if you could do the reaction 10 H + 1 6 Li → 16 O, you still wouldn't exceed the specific energy of the catalyzed D-D+D cycle. Such a reaction is not possible because p-F-19 is the only feasible reaction that produces an nucleus with a mass per nucleon below oxygen-16 and it is not possible to produce fluorine-19 with any feasible reaction.

Next, we will consider option 2A: directly fusing inputs to make products with an average mass per nucleon that is less than helium-4. Figure 3 shows that, because of the exceptional stability of helium-4, it is not enough to just create a larger nucleus we must create a much larger nucleus. Specifically, carbon-12 is the smallest nucleus that has a lower mass per nucleon than helium-4. Producing the nuclei in-between (instead of helium-4) would actually decrease specific energy. The difficulty posed by this is exacerbated by the lack of a nucleus with a mass number of 5. The longest-lived nucleus with A = 5 is helium-5, which has a half-life shorter than 10 -21 seconds [36]. Because of this gap, it is impossible to go past helium-4 by directly fusing inputs. In stars, this gap is bridged by the triple-alpha process, which fuses three helium-4 nuclei into a single carbon-12 nucleus. However, three-body fusion is far too rare to accomplish on human scales. Looking at figure 4, we see no alternatives. No reaction in the first four rows of figure 4 produces a nucleus that is larger than helium-4. Moreover, the easiest reaction in the fifth row is T-He-4, which is a billion times more difficult than D-T. There are reactions that use nuclei with A &gt; 5 that do produce carbon-12 or larger nuclei (e.g. reacting D-B-11). However, using these heavier nuclei as input is unworkable because it decreases the average mass per nucleon of the reactants by too much (as discussed in option 1A). Thus, we cannot use option 2A.

Now we will consider the possibility of a catalytic fusion cycle. In a catalytic cycle, large nuclei are used as inputs to fusion reactions, but are later regenerated by subsequent reactions to eliminate them as inputs to the fuel cycle. Such catalytic cycles play an important role in stellar fusion and some have been proposed for terrestrial fusion [37]. It is important to note that using a catalyst with A &lt; 5 does not help us accomplish our goal due to the discussions of options 1A and 2A. Additionally, no reaction in the entirety of figure 4 produces more than one nuclei with A &gt; 5. This fact directly means that option 2B is impossible. You cannot use catalysts to build products larger than helium-4, because you'll always need the largest product to regenerate the catalyst. This also means that any reaction with an A &gt; 5 reactant, but without an A &gt; 5 product cannot be used in a catalytic cycle. It becomes impossible to regenerate the heavy reactant. This is because (due to the discussion of option 2A), you can only bridge the A = 5 gap with the help of a larger nucleus and you've just lost it. Similarly, you cannot use any reaction where both reactants have A &gt; 5. It is not possible to regenerate both and you cannot tolerate having one as an input because of their low mass per nucleon.

With this context, we consider our last option, 1B: using a catalyst to construct helium-4 with hydrogen-1 as an input. This is what is accomplished by the CNO cycle that occurs in stars. However, it relies on fusing protons with carbon-12, carbon-

13, nitrogen-14, and carbon-15. All but one of these reactions are far beyond our capabilities. However, we should still search for an easier way. We know that we must start by fusing hydrogen-1 with something that has a product with A &gt; 5 (in order to have hope of closing the cycle). We see six options that do this and are at least remotely feasible: beryllium-9, boron-10, nitrogen-15, oxygen-18, fluorine-18, and fluorine-19. Let's examine them in turn.

The beryllium-9 reaction does produce a nucleus with A &gt; 5, but only in one of two branches. Moreover, reconstructing the beryllium-9 requires that only the neutronproducing branch of the D-Li-6 reaction be used (and you must wait 50 days for beryllium-7 to decay into lithium-7). However, if this can be done, it enables the overall reaction 1 H+ H+ H 2 3 → 4 He + 2 n. The hardest constituent reaction of this cycle has 1 a point triple product of 1 6 . × 10 25 keV-s/m . 3 Of course, this is almost a hundred times harder than D-D, but we're considering it to be in the realm of possibility. Another difficulty is obtaining the tritium without ruining the specific energy. The only possibility is to add the branch of D-D fusion that produces tritium, enabling tritium to be eliminated as an input. Then, deuterium breeding using normal hydrogen (i.e. equation (7)) can be used to eliminate one of the input deuterons and we can achieve the reaction 2 H + 1 2 H → 4 He. This fuel cycle has a specific energy of 6.3 MeV/AMU, somewhat higher than catalyzed D-D+D! Unfortunately, apart from all the other difficulties, we have assumed that we can exclude specific branches of three different reactions at will. Simply put, this does not appear possible to the degree that we require. The cross-sections for the different branches have similar magnitudes and depend on temperature in the same way.

The boron-10 reaction can also enable the 1 H + 2 H + 3 H → 4 He + 2 n cycle. 1 Accomplishing this requires the T-Li-7 and D-Be-9 reactions (as well as waiting for beryllium-7 to decay into lithium-7). While this involves only one reaction with branches that must be excluded, it starts from the very hard p-B-10 reaction. Regardless, attaining the 2 H + 1 2 H → 4 He fuel cycle still requires tritium, so the neutron-producing branch of the D-D reaction must be blocked. Thus, starting with boron-10 does not appear any more feasible than starting with beryllium-9.

The last four options are fairly similar. Fusing hydrogen-1 with nitrogen-15, which is part of the solar CNO cycle, produces carbon-12. Unfortunately, carbon-12 is so stable that it doesn't react with anything. Oxygen-18 produces nitrogen-15, which can only be split further, so the cycle cannot be closed. Fluorine-18 produces oxygen-15, which decays into nitrogen-15, but neither can be used to produce a larger nuclei. Finally, fluorine-19 produces oxygen-16 or neon-20, both of which are so stable that they don't react with anything.

In conclusion, we've taken quite a granular look at the possible fusion reactions and been unable to find a viable fuel cycle with a specific energy larger than that of catalyzed D-D+D. Moreover, we've discussed some arguments that suggest it is impossible, given the known fusion cross-sections. The primary features can be summarized as:

- · There are no nuclei with A = 5.

- · Hydrogen-1 cannot be fused with any other nucleus with A &lt; 5.
- · There is no viable reaction between two nuclei with A &lt; 5 that creates a nuclei with A &gt; 5,
- · Using a nucleus with A &gt; 5 as an input decreases the mass per nucleons of the reactants by too much to allow a high specific energy.
- · A nucleus with A &gt; 5 cannot be used as a catalyst because there are few viable fusion reactions involving large nuclei, so it cannot be regenerated without the ability to block certain branches of several reactions.

## 6. Space propulsion

Lastly, a particularly important application of a catalyzed D-D+D fuel cycle could be space propulsion. In space travel, two of the most important performances metrics are specific energy and specific momentum . Here specific momentum is defined to be the total directed momentum of the exhausted particles divided by the mass of the input particles, so it could also be termed an effective exhaust velocity. Specific energy is more important for propulsion schemes with dedicated propellant (e.g. [1, 18, 38]), while specific momentum is key for schemes that directly exhaust the fusion products as propellant (e.g. [18, 25]). Having dedicated propellant increases the total amount of thrust, but also increases the total combined mass of fuel/propellant that is required to achieve a given spacecraft velocity. For example, the rocket equation shows that a spacecraft must be nearly 90% fuel/propellant to reach a speed of double its specific momentum. On the other hand, high thrust is desirable for short journeys and to navigate the gravitational fields within a star system.

Many studies [1, 18, 25] have concluded that a pure D-He-3 fuel cycle is preferable to both catalyzed D-D and D-T. This is sensible because, without deuterium breeding, the D-He-3 reaction is optimal for both specific energy (see table 1) and specific momentum. However, we have just learned that deuterium breeding enables the specific energy of catalyzed D-D to surpass D-He-3. Thus, it has the potential to be useful. In this section, we will first calculate the impact of deuterium breeding on specific momentum and then discuss the overall viability of catalyzed D-D+D propulsion.

Without analysis, it is not clear if a catalyzed D-D+D fuel cycle could surpass the specific momentum of D-He-3. A D-He-3 fusion device would release virtually all of its energy in the form of charged particles, which can be efficiently redirected into a beam to maximize thrust. This enables D-He-3 fuel to achieve a specific momentum approaching 7.1% of the speed of light. On the other hand, the catalyzed D-D+D fuel cycle is theoretically capable of a specific momentum of 11.3% of the speed of light. However, this assumes that the 23.8 MeV released by the chain of reactions underlying equation (8) can be entirely transferred to the helium-4 product. In practice, this does not appear possible because much of this energy is given to neutrons, which is difficult to use efficiently. Nevertheless, how close we can get involves a complex design optimization that depends strongly on the practicality and efficiency of various

Figure 6: A cartoon fusion-powered spacecraft with neutron pusher/blanket extending over φ ∈ [0 , π/ 2] and the full [0 , 2 π ] in azimuthal angle. As an example, a magnetic mirror is used as the confinement scheme.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00834v1-figure-92c14c6b-d287-44b7-aba4-c96578c08d30.jpg)

technological components. Given that net energy production with any fusion fuel has yet to be achieved, a detailed investigation is premature. Instead, we will present a few idealized examples to convey a sense of the possibilities and limitations.

In conventional catalyzed D-D fusion the fuel cycle can be written as

$$6 \, ^ { 2 } H \to & ^ { 4 } \text{He} \, ( 3. 5 \, \text{MeV} ) + ^ { 1 } H \, ( 3. 0 2 \, \text{MeV} ) + ^ { 1 } n \, ( 1 4. 1 \, \text{MeV} ) & ( 1 3 ) \\ & + ^ { 4 } \text{He} \, ( 3. 6 \, \text{MeV} ) + ^ { 1 } H \, ( 1 4. 7 \, \text{MeV} ) + ^ { 1 } n \, ( 2. 4 5 \, \text{MeV} ) + 1. 8 3 \, \text{MeV},$$

where the kinetic energy carried by each particle is indicated in the parentheses and the energy that is carried by intermediate products (which remains on-board) is indicated by the last term. This shows that 38% of the total energy is carried off by neutrons and 4% is carried by the 3 H and 3 He fusion products, which must remain in the spacecraft to fuse. While we assume that charged products can be exhausted optimally, we will assume that neutrons are born isotropically and cannot be redirected. Consequently, only 1 / 4 of the neutron's total momentum can be utilized for thrust. This is accomplished by extending a barrier from 0 to π/ 2 in the zenith angle φ (see figure 6) in order to block all of the neutrons traveling forward and allow all of the neutrons traveling backwards to freely stream away. Even with such geometry, the normal catalyzed D-D fuel cycle achieves a specific momentum of only 5.6% of the speed of light.

Introducing a blanket of normal hydrogen that extends over the full φ ∈ [0 , π ] enables deuterium breeding and achieves the fuel cycle

$$4 \, ^ { 2 } H \to ^ { 4 } \text{He } ( 3. 5 \, \text{MeV} ) + ^ { 4 } \text{He } ( 3. 6 \, \text{MeV} ) + 4 0. 5 \, \text{MeV}, \quad \quad ( 1 4 )$$

which is equivalent to equation (8). Even assuming that the helium-4 is exhausted optimally, this fuel cycle actually has a lower specific momentum than normal catalyzed

D-D fusion - only 4.4% of the speed of light. This is because the normal hydrogen produced by fusion, which previously carried significant amounts of thrust, is now being kept on-board to breed deuterium. However, we see that there is a tremendous amount of energy now available on the spacecraft. This is bad if waste heat removal is difficult [25, 39] or good if additional power is useful to run the confinement scheme/other onboard systems. For example, one could imagine using this power to directly augment the thrust by adding an ion acceleration stage at the end of the fusion thruster. If the on-board energy could be transferred to the outgoing helium-4 with an overall efficiency of just 30%, a specific momentum of 7.2% of the speed of light could be achieved. This would require accelerating particles from 3.5 MeV to 9.5 MeV, but would enable the catalyzed D-D+D fuel cycle to surpass the specific momentum of D-He-3. An overall combined fusion product-to-acceleration efficiency of 30% does not appear impossible. Much of the on-board energy is carried by charged particles, which can theoretically be directly converted to electricity at high efficiency, while the rest could go through a steam cycle. The feasibility and efficiency of high-energy ion acceleration is very uncertain, but low-energy ion thrusters currently have efficiencies as high as 80% [40]. Note that, in the limit of a 100% overall fusion product-to-acceleration efficiency, the theoretical maximum of 11.3% of the speed of light is achieved.

Without ion acceleration, it is better to exhaust the normal hydrogen produced by fusion and breed deuterium with hydrogen that has been brought along. This gives the fuel cycle

$$4 ^ { \, 2 } H + 2 ^ { 1 } H \to & ^ { 4 } \text{He} \left ( 3. 5 \text{ MeV} \right ) + ^ { 1 } H \left ( 3. 0 2 \text{ MeV} \right ) \\ & + ^ { 4 } \text{ He} \left ( 3. 6 \text{ MeV} \right ) + ^ { 1 } H \left ( 1 4. 7 \text{ MeV} \right ) + 2 2. 8 \text{ MeV},$$

which has a specific momentum of 6.1% of the speed of light - better than the standard catalyzed D-D, but less than D-He-3. Alternatively, for missions involving very long travel times, one could modify the catalyzed D-D cycle somewhat by waiting for the tritium to decay into helium-3 before fusing it [41]. When combined with deuterium breeding it could enable the fuel cycle

$$5 ^ { \ 2 } H \to 2 ^ { \ 4 } \text{He } ( 3. 6 \text{ MeV} ) + 2 ^ { \ 1 } \text{H } ( 1 4. 7 \text{ MeV} ) + 9. 5 2 \text{ MeV}.$$

This is very similar to the standard D-He-3 reaction and has the same specific momentum of 7.1% of the speed of light. However, it has the advantages of only requiring deuterium as input and generating extra on-board energy. The disadvantages include a massive breeding blanket/radiator, waiting ∼ 12 years for tritium to decay, and an overall increased level of technical complexity.

Though even more speculative, what would significantly improve the specific momentum would be to collect normal hydrogen along the flight path. This would allow a fuel cycle to directly exhaust the hydrogen-1 produced by fusion, without having to bring it along as an input. For example, eliminating the input hydrogen-1 in equation (15) boosts the specific momentum from 6.1% to 7.6% - somewhat larger than D-He-3.

Moreover, there is still 22.8 MeV of on-board energy that could be used to power the hydrogen collection scheme and/or for ion acceleration. Using ion acceleration with a fusion product-to acceleration efficiency of 30% could boost the specific momentum to 8.9% of the speed of light.

The idea of collecting hydrogen is similar to the Bussard ramjet [42, 43], which uses large electromagnetic fields to gather particles from kilometers around ¶ . While the design detailed in this work requires a similar number of particles to be collected, it is much more practical than the Bussard ramjet because it does not depend on p-p fusion (which has an extremely small cross-section [35]), nor the collection of deuterium (which is tens of thousands of times less abundant than normal hydrogen [26]). The catch is that getting the collected hydrogen into the breeding blanket would require it to be accelerated to the speed of the spacecraft. During the spacecraft's acceleration, any savings that were gained by not bringing the fuel would be entirely canceled by drag against the particles that must be gathered. However, during deceleration it could be doubly helpful - the drag would further decelerate the spacecraft, while the collected hydrogen could be used to breed deuterium.

In summary, we see that adding deuterium breeding can substantially improve the specific momentum of a catalyzed D-D propulsion system, potentially enabling it to surpass D-He-3. Importantly, the primary advantages of a catalyzed D-D fuel cycle are retained. Catalyzed D-D uses fuel that is extremely abundant and readily available on Earth and produces plasmas that have a higher power density and ignite at a lower temperature than D-He-3 [15, 16]. However, we have seen that a D-D propulsion system and the techniques needed to improve its performance carry the drawbacks of increased technical complexity and the need for massive components (e.g. neutron shielding, waste heat radiators, ion acceleration grids, hydrogen collection systems). Thus, the ultimate practicality of a fuel cycle will depend strongly on factors like how light-weight these components are and how well they perform.

## 7. Conclusions

This work has tried to answer the question 'What fuel source can provide the most energy with the least mass?' In doing so, we have considered an addition to the catalyzed D-D fusion fuel cycle. By surrounding the reactor with a blanket containing normal hydrogen, the neutrons from fusion can be used to breed deuterium and eliminate a substantial fraction of the fuel needed to run the power plant. Given the ability to achieve a catalyzed D-D cycle, this addition appears straightforward to accomplish - one could simply immerse the reactor in a tank of liquid water. By studying experimentally measured fusion cross-sections, we have argued that such a scheme enables the maximum specific energy of any known fuel source. Apart from satisfying

¶ Catalyzed D-D+D fusion could also improve the performance of a related concept: the RAIR [44, 45, 46]. On-board fuel is used to power the spacecraft, but the propellant is gathered from interstellar space, so the specific energy (rather than momentum) is what matters.

our own curiosity, such an answer reveals how to most efficiently use the dominant energy stockpile we have here on Earth and in the universe. Moreover, such a scheme could be useful for any application that requires fuel to be transported. We have analyzed one such application, space propulsion, and seen that a catalyzed D-D thruster using deuterium breeding is theoretically capable of surpassing the performance of D-He-3.

## Acknowledgments

The author is indebted to J. Parisi and H. de Oliveira for valuable discussions and comments.

## References

- [1] G.W. Englert. Towards thermonuclear rocket propulsion. New Scientist , 16:16, 1962.
- [2] R. Zubrin. Entering Space , pages 90, 160, 161, 197, 198, 200. Jeremy P. Tharcher/Putnam, 1999.
- [3] F.W. Stecker, D.L. Morgan Jr, and J. Bredekamp. Possible evidence for the existence of antimatter on a cosmological scale in the universe. Physical Review Letters , 27(21):1469, 1971.
- [4] R.F. Post. Sixteen lectures on controlled thermonuclear reactions: Part VII. Technical Report UCRL-4231, Univ. of California Radiation Lab., 1954.
- [5] A. Simon. Project Sherwood: Orientation lectures, chapter 1. Technical Report ORNL-CF-56-8140, Oak Ridge National Lab, 1955.
- [6] R.F. Post. Controlled fusion research - an application of the physics of high temperature plasmas. Reviews of Modern Physics , 28(3):338, 1956.
- [7] G.H. Miley. Overview of alternate-fuel fusion. Technical report, Illinois Univ., Urbana (USA). Fusion Studies Lab, 1980.
- [8] F. Bello. Fusion power: The trail gets hotter. Fortune , page 135, 1957.
- [9] M.T. Robinson. On the chemistry of thermonuclear reactor breeder blankets. Technical report, Oak Ridge National Lab., Tenn., 1958.
- [10] L. Spitzer Jr, D.J. Grove, W.E. Johnson, L. Tonks, and W.F. Westendorp. Problems of the stellarator as a useful power source. Technical Report NYO-6047, Project Matterhorn, Princeton University, 1954.
- [11] W. Davis Jr. Some chemical problems associated with thermonuclear reactors. Technical Report CF-58-11-95, Oak Ridge National Lab., Tenn., 1958.
- [12] L.J. Wittenberg. Non-lunar He-3 resources. In 2nd Wisconsin Symp. Helium-3 and Fusion Power , page 19, July 1993.
- [13] L.J. Wittenberg, E.N. Cameron, G.L. Kulcinski, S.H. Ott, J.F. Santarius, G.I. Sviatoslavsky, I.N. Sviatoslavsky, and H.E. Thompson. A review of 3 He resources and acquisition for use as fusion fuel. Fusion Technology , 21(4):2230-2253, 1992.
- [14] G.H. Miley. 3 He sources for D- He fusion power. 3 Nuclear Instruments and Methods in Physics Research Section A , 271(1):197-202, 1988.
- [15] R.G. Mills. Catalyzed deuterium fusion reactors. Technical Report TM-259, 1971.
- [16] G.H. Miley, F.H. Southworth, G. Gerdin, and C. Choi. Studies of deuterium-fueled tokamak reactors. Nuclear Instruments and Methods , 144(1):9-16, 1977.
- [17] T.H. Rider. A general critique of inertial-electrostatic confinement fusion systems. Physics of Plasmas , 2(6):1853, 1995.
- [18] J.L Hilton, J.S. Luce, and A.S. Thompson. Hypothetical fusion propulsion rocket vehicle. Journal of Spacecraft and Rockets , 1(3):276, 1964.
- [19] E.J. Lerner. Prospects for P 11 B fusion with the dense plasma focus: New results. Current Trends in International Fusion Research , page 25, 2008.

- [20] X.-P. Huang, F. Anderegg, E.M. Hollmann, C.F. Driscoll, and T.M. O'Neil. Steady-state confinement of non-neutral plasmas by rotating electric fields. Physical Review Letters , 78(5):875, 1997.
- [21] N. Rostoker, M.W. Binderbauer, and H.J. Monkhorst. Colliding beam fusion reactor. Science , 278(5342):1419-1422, 1997.
- [22] H. Hora, S. Eliezer, J. Wang, G. Korn, N. Nissim, Y.-X. Xu, P. Lalousis, G.J. Kirchhoff, and G.H. Miley. Laser boron fusion reactor with picosecond petawatt block ignition. IEEE Transactions on Plasma Science , 46(5):1191-1197, 2018.
- [23] M. Keilhacker, A. Gibson, C. Gormezano, P.J. Lomas, P.R. Thomas, M.L. Watkins, P. Andrew, B. Balet, D. Borba, C.D. Challis, and others. High fusion performance from deuterium-tritium plasmas in JET. Nuclear Fusion , 39(2):209, 1999.
- [24] B. Hamermesh, G.R. Ringo, and S. Wexler. The thermal neutron capture cross section of hydrogen. Physical Review , 90(4):603, 1953.
- [25] J.F. Santarius. Lunar He-3, fusion propulsion, and space development. In 2nd Conference on Lunar Bases and Space Activities , page 75, 1992.
- [26] J.M. O'Meara, D. Tytler, D. Kirkman, N. Suzuki, J.X. Prochaska, D. Lubin, and A.M. Wolfe. The deuterium to hydrogen abundance ratio toward a fourth QSO: HS 0105+ 1619. The Astrophysical Journal , 552(2):718, 2001.
- [27] M. Shimada, D.J. Campbell, V. Mukhovatov, M. Fujiwara, N. Kirneva, K. Lackner, M. Nagami, V.D. Pustovitov, N. Uckan, J. Wesley, and others. Chapter 1: Overview and summary. Nuclear Fusion , 47(6):S1, 2007.
- [28] A.R. Raffray, L. El-Guebaly, S. Malang, I. Sviatoslavsky, M.S. Tillack, X. Wang, ARIES Team, et al. Advanced power core system for the ARIES-AT power plant. Fusion Engineering and Design , 80(1-4):79, 2006.
- [29] B.N. Sorbom, J. Ball, T.R. Palmer, F.J. Mangiarotti, J.M. Sierchio, P. Bonoli, C. Kasten, D.A. Sutherland, H.S. Barnard, C.B. Haakonsen, J. Goh, C. Sung, and D.G. Whyte. ARC: A compact, high-field, fusion nuclear science facility and demonstration power plant with demountable magnets. Fusion Engineering and Design , 100(378), 2015.
- [30] X-5 Monte Carlo Team. MCNP - a general n-particle transport code, version 5 - volume I: Overview and theory. LA-UR-03-1987, Los Alamos National Laboratory , 2003.
- [31] N. Otuka, E. Dupont, V. Semkova, B. Pritychenko, A.I. Blokhin, M. Aikawa, S. Babykina, M. Bossant, G. Chen, S. Dunaeva, et al. Towards a more complete and accurate experimental nuclear reaction data library (EXFOR): international collaboration between nuclear reaction data centres (NRDC). Nuclear Data Sheets , 120:272-276, 2014.
- [32] J. Wesson. Tokamaks. page 10. Oxford University Press, fourth edition, 2004.
- [33] J.D. Lawson. Some Criteria for a Power Producing Thermonuclear Reactor. Proceedings of the Physical Society. Section B , 70(1):6-10, 1957.
- [34] J.D. Huba. NRL Plasma Formulary . Naval Research Laboratory, 2013.
- [35] E.G. Adelberger, S.M. Austin, J.N. Bahcall, A.B. Balantekin, G. Bogaert, L.S. Brown, L. Buchmann, F.E. Cecil, A.E. Champagne, L. De Braeckeleer, et al. Solar fusion cross sections. Reviews of Modern Physics , 70(4):1265, 1998.
- [36] E.M. Baum, H.D. Knox, and T.R. Miller. Nuclides and isotopes: chart of the nuclides . KAPL, 2002.
- [37] J.R. McNally Jr. Fusion chain reactions - I. Nuclear Fusion , 11(2):187, 1971.
- [38] G. Miley, L. Wu, H. Momota, and J. Santarius. On D-He-3 fusion reactor development for fusion propulsion. In 5th International Energy Conversion Engineering Conference (IECEC) , page 4749, 2007.
- [39] S.H. Maslen. Fusion for space propulsion. IRE Transactions on Military Electronics , 1051(2):5257, 1959.
- [40] E.Y. Choueiri. New dawn for electric rockets. Scientific American , 300(2):58-65, 2009.
- [41] J. Kesner, D.T. Garnier, A. Hansen, M. Mauel, and L. Bromberg. Helium catalysed D-D fusion

in a levitated dipole. Nuclear Fusion , 44(1):193, 2003.

- [42] R.W. Bussard. Galactic matter and interstellar flight. Astronautica Acta , 6(4):179-194, 1960.
- [43] D.P. Whitmire. Relativistic spaceflight and the catalytic nuclear ramjet. Acta Astronautica , 2(5-6):497-509, 1975.
- [44] A. Bond. An analysis of the potential performance of the ram augmented interstellar rocket. Journal of the British Interplanetary Society , 27:674, 1974.
- [45] C. Powell. Flight dynamics of the ram-augmented interstellar rocket. Journal of the British Interplanetary Society , 28:553, 1975.
- [46] A. Jackson. Some considerations on the antimatter and fusion ram augmented interstellar rocket. Journal of the British Interplanetary Society , 33:117-120, 1980.
Justin Ball
Maximizing specific energy by breeding deuterium
Specific energy (i.e. energy per unit mass) is one of the most fundamental
and consequential properties of a fuel source. In this work, a systematic study
of measured fusion cross-sections is performed to determine which reactions are
potentially feasible and identify the fuel cycle that maximizes specific
energy. This reveals that, by using normal hydrogen to breed deuterium via
neutron capture, the conventional catalyzed D-D fusion fuel cycle can attain a
specific energy greater than anything else. Simply surrounding a catalyzed D-D
reactor with water enables deuterium fuel, the dominant stockpile of energy on
Earth, to produce as much as 65% more energy. Lastly, the impact on space
propulsion is considered, revealing that an effective exhaust velocity
exceeding that of deuterium-helium-3 is theoretically possible.