N. Rostoker, M. W. Binderbauer, and H. J. Monkhorst (1) revisit several ideas—the use of the reactions of protons with boron-11 (p-^{11}B) to produce fusion energy, (2) the use of the field-reversed configuration (FRC) for magnetic confinement of particles and energy (*E*) (3), and the use of nonthermal ion distributions to enhance the fusion reactivity (4)—in their proposal for a colliding beam fusion reactor (CBFR). While there are unresolved issues with each of these choices, I focus here on nonthermal ion distributions.

The fundamental difficulty with nonthermal ions is apparent when one compares the fusion cross section (σ_{fusion} ≈ 1 barn for p-^{11}B at *E*
_{cm} = 580 keV, where*E*
_{cm} is the energy in the center-of-mass frame) to the effective cross section for many small-angle Coulomb scattering events that combine to produce a scattering angle of 90° rms in an incident beam (σ_{eff} ≈ 60 barns for protons scattering on ^{11}B at *E*
_{cm} = 580 keV). Highly nonthermal systems, like the colliding beam reactor proposed by Rostoker *et al*. would relax to local thermal equilibrium before a significant amount of fusion power could be produced. Alternatively, the nonthermal ion distribution could be maintained by cycling sufficient power through the system.

With this possibility in mind, it is instructive to examine the operating point presented by Rostoker *et al*. (*n*
_{p} = 4 × 10^{15}/cm^{3},*n*
_{B} = 1 × 10^{13}/cm^{3},*T*
_{p} = 25 keV, and *T*
_{e} = 20 keV), where *n*
_{p} is the proton density,*n*
_{B} is the boron density,*T*
_{p} is the temperature of the protons, and *T*
_{e} is the temperature of the electrons. We consider the problem in the frame-of-reference of the protons. Maintaining *E*
_{cm} = 580 keV would require that the boron beam velocity be *u*
_{B} ≈ 1.1 × 10^{9} cm/s. Coulomb collisions with the boron beam would heat the proton distribution at a rate (5) of 1/2*m*
_{B}
*u*
_{B}
^{2}
*n*
_{B}υ_{ɛ}
^{B/p}≈ 500 W/cm^{3}, where *m*
_{B} = 11 AMU is the boron mass and υ_{ɛ}
^{i/j} is the energy transfer rate from collisions between particles of species “i” and “j.” As Rostoker *et al*. suggest, heat would be removed from the protons by collisions with the (colder) electrons at a rate of 3/2*n*
_{p}(*T*
_{p} −*T*
_{c}) υ_{ɛ}
^{p/c} ≈ 350 W/cm^{3}. The relatively small net heating of the proton distribution could be eliminated if *T*
_{e} were reduced to 18.6 keV.

The largest term in the electron power balance is direct heating by the boron beam, which would proceed at the rate of 1/2*m*
_{B}
*u*
_{B}
^{2}
*n*
_{B}υ_{ɛ}
^{B/e}≈ 2.0 kW/cm^{3}. Maintaining the electrons at*T*
_{e} ≈ 18.6 keV would require that*P*
_{electron} ≈ 2.5 kW/cm^{3} be extracted from the electron distribution. While bremsstrahlung could provide only 44 W/cm^{3} of this electron cooling, experience suggests that it would not be difficult to find other channels for electron heat loss. The 2.5 kW/cm^{3} extracted from the electrons must be balanced by ion heating.

The fusion power (*P*
_{fusion}) calculated by Rostoker *et al*. at the operating point is*P*
_{fusion} =*n*
_{p}
*n*
_{B}
*Y*
_{p-11}
_{B}〈σ_{p-11}
_{B}
*v*〉 ≈ 45 W/cm^{3}, where*Y*
_{p-11B} = 8.7 MeV is the energy yield per fusion event (and *v* is the velocity of the proton and boron nuclei). The difference between*P*
_{electron} and *P*
_{fusion}sets a lower limit on the external ion heating power. Thus, the fusion gain (which is normally defined as the ratio of the fusion power to the external heating power), *Q ≡ P*
_{fusion}/*P*
_{heat} ≤*P*
_{fusion}/(*P*
_{electron}− *P*
_{fusion}) ≈ 0.02 would be much lower than the value of 2.7 stated by Rostoker *et al*. (1).

Highly nonthermal fusion schemes generally suffer from an unattractive power balance (4), and so magnetic fusion research has focused on systems that are in local thermal equilibrium. In such systems, high fusion gain (*Q* >> 1) remains a distinct possibility.

**R**ecognizing the unfavorable power balance of a thermal proton-boron (p-^{11}B) plasma, N. Rostoker *et al*. (1) propose restricting the proton energy relative to the boron to be near the resonance in the fusion cross section at E_{0}±ΔE = (580 ± 140) keV. Although this beam-beam configuration would avoid the large power required (2) in the migma (3) to maintain the proton distribution against self-collisions, a large power input would nevertheless be required to replace the directed energy lost to frictional heating of the proton and boron beams. The classical formula (4) for this power density iswhere ln Λ ≈ 15 (or larger) is the Coulomb logarithm and Z_{B} is the atomic number of boron, 5. Because the fusion power density is *P*_{fusion} = (8.7 MeV) 〈συ〉 *n*_{p}*n*_{B}, if we set the reactivity equal to its value at the resonance, enhanced by a factor of 1.6 for spin polarization, 〈συ〉 = 1.3 × 10^{−21} m^{3} s^{−1}, then we find the ratio to be P_{fusion}/ P_{fric}≈ 0.12, which shows that the proposal for a reactor with net electrical power output is unrealistic.

The power balance would be at least another factor of three less favorable than this estimate because the coupling of the ions through the electrons would be stronger than the direct coupling if*T*
_{e}
* < E*
_{0}/15 = 40 keV. If the electron coupling were decreased by raising*T*
_{e}, *T*
_{p} must also rise. Rostoker *et al*. suggest that the*T*
_{p} has only to be less than 140 keV, but this is not sufficiently cool. In order for protons with velocity*T* must be κ_{B}
*T*
_{p} ≲ (Δ*E*)^{2}/2*E*
_{0} = 17 keV. Rostoker *et al*. published a similar unworkable reactor design earlier (5) with κ_{B}
*T*
_{p} = 200 keV.

Next to power balance, a serious problem of the CBFR is equilibrium. The plasma volume envisaged would be a long, thin cylindrical shell with thickness/radius Δ r/r ≈ 0.08. Such a configuration would not be in axial equilibrium because the tension of the field lines curving around the shell at the ends would provide a powerful compressive force. For 2-dimensional equilibria, Δ r/r ≈ 0.5 (6), so the highly localized profiles required to prevent radial particle and energy losses would not be maintained.

*Response*: The main point of the comment by Nevins is that direct heating of electrons from Coulomb scattering by the boron beam (2.0 kW/cm^{3}) is the largest term in the electron power balance, and that we neglected it in our article (1). The electrons would have to be maintained at a low T in order to cool the proton beam. The problem is with the calculation of the electron T , which is not detailed in our article or in the comment by W. M. Nevins. The calculation of heating of electrons by boron scattering does give a large result. However, electrons lose energy to the boron ions because their velocity would be higher and there would be a collisional drag effect. This energy loss would cancel the energy gain. The electron temperature*T*_{e} is defined bywhere *n*_{e}is the electron density,*V⃗*_{e} is the average electron velocity, f_{e}(r⃗,υ⃗) is the electron distribution function, and*N*_{e} is the line density of electrons, which is constant. The calculation of*dT*_{e}*/dt* involves the Fokker-Planck equation and is complicated. If it is done in a straightforward way with the above definition, the result will be as Nevins describes. However, if we note that in a steady state*dV*_{e}*/dt* = 0 and calculate only *d〈υ ^{2}〉/dt*, the cancellation takes place, and electrons would only be heated by the protons. Similar considerations show that protons are cooled by electrons and boron.

We have revisited ideas such as “the use of p-^{11}B as a fuel.” However, we have developed new and systematic calculations and the conclusions are not the same as they were with the use of that proposed fuel.

Carlson employs a classical generic formula for the power density required to overcome the friction between proton and boron beams. This formula is inadequate for the Colliding Beam Fusion Reactor. The magnetic field is important, and it is distinguished by its absence in this formula. The complete formula can be derived by taking the appropriate moment of the Vlasov/Fokker-Planck equation. First, the equilibrium conservation of momentum equation iswith i = (1, 2) for protons and boron. Here n_{i} is the particle density; e_{i} is the charge; m_{i} is the mass; T_{i} is the temperature; E_{r} and B_{z} are the electric and magnetic field, respectively; and V_{i} is the velocity in the azimuthal direction. The kinetic equation obtained from the same moment, but including the collision operator, is
with a similar equation forV⃗_{2}. The dot product ofV⃗_{1} with the above equations yields an energy (power) equation and the magnetic field seems to vanish. The power is
The first term is precisely the expression for power density employed by Carlson. The second term involves the radial velocity *V*
_{ir} and is positive definite. It can be estimated with some approximations;*V*
_{ir} <<*V*
_{iθ} and*n*
_{i}
*V*
_{ir} = −*D*
_{i}(∂*n*
_{i}/∂_{r}), where *D*
_{i} is the diffusion coefficient; D_{1} ≃ a_{1}
^{2}/ t_{12}, where a_{1} = V_{1}/Ω_{1}is the gyro-radius and Ω_{1} = e_{1}| B_{z}
*|/m*
_{1}c . Similar expressions obtain for boron. Although V_{ir} may be neglected as compared with*V*
_{iθ}, it may not be neglected in the second term. This term can be estimated from the equilibrium equation. One can see that the magnetic field that previously seemed to vanish has returned. The second term is positive definite and the magnitude n_{1}m_{1}(a_{1}/L_{1})(V_{1}
^{2}/t_{12}) is similar to that of the first term. L_{1} is the scale length of the equilibrium, that is (1/ L_{1}) = (1/n_{1})|dn_{1}/dr |. To determine this power quantitatively requires a considerable amount of work (2). It requires a new development in classical transport theory because earlier studies assume a_{i}<< L_{1}, which would not be the case in the CBFR. This calculation should also include electrons and the fusion products. The result is that*P*
_{fric} would be tolerable.

Concerning the resonance, we have made detailed calculations. If the beam temperatures are less than one-half of the half width of the resonance, the reactivity should be greater than one-half of the maximum reactivity for zero temperature beams. The result stated by Carlson that the beam temperature must be less than 17 keV seems to contradict this. However, no results for 〈συ〉 are given for comparison.

The equilibrium calculations to which Carlson refers are not appropriate for the CBFR. Axial equilibrium requires an axial T , or the FRC will contract in the axial direction. It has been observed experimentally that FRCs have long axial equilibria. We have previously considered long, thin cylindrical shell models because they simplify many calculations.