Published online by Cambridge University Press: 25 March 2004
Bremsstrahlung is one of the most important energy loss mechanisms in achieving ignition, which is only possible above a threshold in temperature for a given fusion reaction and plasma conditions. A detailed analysis of the bremsstrahlung process in degenerate plasma points out that radiation energy loss is much smaller than the value given by the classical formulation. This fact seems not useful to relax ignition requirements in self-ignited targets, because it is only relevant at extremely high densities. On the contrary, it can be very positive in the fast ignition scheme, where the target is compressed to very high densities at a minimum temperature, before the igniting beamlet is sent in.
The classical approach to produce ignition in fusion plasmas is to heat them up to a temperature high enough for the fusion-born charged particle energy to overcome radiation emission. This approach is followed both in magnetic confinement, where plasmas are underdense, and inertial confinement, where plasmas are compressed to very high densities.
Although standard schemes (Bodner 1981, Brueckner & Jorna, 1973; Kidder, 1979; Meyer-ter-Vehn, 1982; Martìnez-Val et al., 1994; Tahir & Hoffmann, 1994; Atzeni, 1995, 2002) of inertial fusion research do not take into account plasma degeneracy, inertial compression could produce them according the lines established in the shock-multiplexing implosion scheme (Kidder, 1974; Yamanaka et al., 1986, 1987; Yamanaka & Nakai, 1986), which got 600 times compression with moderate laser beam energy (Nakai et al., 1991) in a series of experiments carried out at the end of the 1980s, particularly with spherical polymer shells, where a compression factor of 1000 times was achieved (Azechi et al., 1991). However, the interest of inertial fusion research focused on different implosion schemes, yielding moderate densities (∼20 g/cm3) with high electron temperatures (∼1 keV), which correspond to a classical plasma (Yaakobi et al., 1999). Nevertheless, the same laboratory had obtained higher densities in compressed targets 10 years before (McCrory et al., 1988) depending on the tailoring of the implosion pulse and the target structure.
In the standard approaches formerly cited, targets were self-ignited either from a central hot spark or from a larger domain (volume ignition). A different scheme was proposed in the fast ignition model (Tabak et al., 1994), where the target is highly compressed at a temperature as low as possible, and ignition is triggered by sending a powerful ultrashort beam of the required energy. In this case, the interaction phase between the ignition beam and the target is of primary importance (Deutsch, 1997; Kato et al., 1997; Piriz & Sanchez, 1998; Boreham et al., 1999; Key, 2001; Kodama et al., 2001; Norreys et al., 2000; Roth et al., 2001) but is worth pointing out that the compressed target can be in a degenerate state when the igniting beam is launched. Although the igniting region will be heated to a very high temperature in a very short time, a significant fraction of the reacting plasma can still be at a degenerate state. It is therefore very important to analyze radiation emission from degenerate plasmas, in the framework of fast ignition.
Radiation losses in inertial targets are dominated by bremsstrahlung. Inverse Compton can additionally transfer energy from electrons to photons (Eliezer & Martinez-Val, 1998), but this mechanism is not significant if the electron temperature is moderate or low, as is the case in degenerate plasmas. A comprehensive review of atomic physics processes in very dense plasmas can be found elsewhere (More, 1993).
The possibility of using degenerate plasmas for inertial fusion was recently outlined in a macroscopic approach (Leon et al., 2001) but a more accurate and detailed analysis for the relevant physical mechanisms inside the plasma was needed. In the following, radiation loss by bremsstrahlung emission is studied from the starting point of electron-photon energy transfer through its integration in the Fermi–Dirac distribution. It must be taken into account that most of the low level energy states in the electron statistics are already occupied, and a very broad range of electron transitions are forbidden.
This subject was already studied by Totsuji (1985) with a general purpose, not specifically oriented to fusion applications. Besides that, in our analysis, a clear connection between the degenerate and the classical case is obtained. It is clearly seen that there is a transition temperature, for a given plasma density, below which the degenerate emission rate applies. On the other hand, bremsstrahlung emission follows the classical formula for higher temperatures.
After analyzing in detail bremsstrahlung emission in degenerate plasmas, the corresponding equation is applied to the definition of ignition temperature. In targets compressed to very high densities, as those sought in the fast ignition concept, it is found that the ignition temperature is lower than the classical one, because of bremsstrahlung emission reduction.
The spectral distribution of radiation emitted by a monoenergetic electron beam of unit intensity scattered by an ion of charge Ze, is (Miyamoto, 1980) in SI units
for hυ ≤ ½me2u2, being 0 for hυ ≥ ½me2u2. In the foregoing expressions, υ is the photon frequency, and u stands for the electron speed. No relativistic correction has been included because the electron temperature Te (in energy units) is much lower than me c2.
It must be noted that Bdυ has units of J·m2, and it can be integrated to give the radiation power density w(ν) dν emitted from a plasma with an ion density ni and an electron density ne. In the case of having multiple ion species, a summation over the ion species must be extended. We first recall that the classical integration is
where umin = (2hν/me)1/2 is the minimum electron speed to generate an hν photon, Te is given in energy units and f (u,Te) is the classical Maxwellian distribution.
In the degenerate case, the energy electron distribution follows the Fermi–Dirac expression (see, e.g., Eliezer et al., 1986):
where the Fermi energy, εF (known also as the chemical potential) plays a fundamental role, and depends on electron density and Te.
For Te << εF0,
where
It is worth noting that the probability of state with energy ε to be occupied is
In particular, for Te = 0, all energy states ε < εF are fully occupied, whereas all states ε > εF are empty. When Te increases (still in the degenerate case, and therefore Te < εF) transition from ε < εF to ε > εF becomes not so sharp, and some electrons can be found with ε > εF. Because of the forbidden electron transitions in the degenerate case, Eq. (3) is no longer valid, and ne f (u,Te) must be substituted by a term properly taking into account that most of the low energy states are fully occupied, and an electron cannot fall into them after emitting a photon. Hence, ne f (u,Te) in Eq. (2) must be replaced by
where ne(ε,Te) is given by Eq. (4).
The total radiation loss expressed in terms of power density is
where the spectral power density is now
and
The integrand of Eq. (10) is depicted in Figure 1. The degeneracy effect is clearly seen.
The spectrum is obtained by integrating Eq. (10). Defining the degeneracy parameter by
one gets for the spectrum
and F1/2(η) is the corresponding Fermi–Dirac integral function. The total bremsstrahlung power density in a degenerate plasma, characterized by ne and Te is obtained by integrating Eq. (13)
Fermi–Dirac integrals are given below, with the corresponding limits for η >> 1 (degenerate case)
If bremsstrahlung emission in the degenerate case is compared to the value of the classical case, we find the exact result
For η >> 1, which corresponds to a fully degenerate case, the last equation can be rewritten as
For instance, if Te = 0.1εF in the degenerate plasma, its bremsstrahlung loss will be about 7% of the radiation loss from a classical plasma at the same temperature. However, this comparison is not totally accurate. If the plasmas have the same temperature, the density of the degenerate case will be much larger than that of the classical case. If they have the same density, they will have different temperatures (much lower in the degenerate case). A more comprehensive analysis is presented in the following section.
In our analytical model, the fully degenerate case corresponds to η >> 1 [Eq. (12)]. Note that the electron density can be written as
If ne and Te are given, the value of the Fermi–Dirac integral F1/2(η) is derived from them, using the last equation. The value of η is then calculated from Eq. (15). In very degenerate cases, the approximate value (3/2F1/2)2/3 can be used. Once η is known, the Fermi energy εF(Te) is directly derived from Eq. (12).
When using the degenerate model, the classical limit is theoretically achieved for η → −∞. This is so because of the difference between Fermian and Maxwellian energy distributions. It can be checked that the degenerate model goes exactly to the classical one as η → −∞. Equation (17) equals one in this limit.
In the classical case, bremsstrahlung emission goes as Te1/2. In the degenerate case, it goes as Te2. Figure 2 shows the bremsstrahlung power density for various density values, as a function of Te. Below a certain Te value, for a given density, bremsstrahlung emission follows the degenerate case expression. On the contrary, for upper temperatures, it corresponds to the classical case. A short Te range is observed as the transition regime from one to another. Of course, the transition temperature increases with density.
Additional information is given in Figure 3, where the ration of the radiation emission in the degenerate case to the value of the classical case is depicted. It is seen that the ratio becomes 1 above the transition temperature, and decreases very fast as temperature decreases below it. Such a feature will be analyzed later on in order to find out whether it is possible to reduce the fusion ignition temperature by compressing the targets up to degenerate levels.
Spectral differences are also observed between classical and degenerate cases, as can be seen in Figures 4 and 5. In Figure 4, spectral distribution of photons is given for a classical case and a degenerate case, both of them with an electron density of 1.81 × 1023 cm−3. Fermi energy is in this example εFO = 11.68 eV. In the degenerate case, the temperature is 3.5 eV (η = 3). In the classical one, it is 30 eV. The spectral shapes (in semilog scale) are the same at medium and high photon energies, but they are different at low values, where the degenerate case does not have a constant slope, which is a typical feature of the bremsstrahlung spectrum from classical plasmas. Therefore, it is worth pointing out that the degenerate case has a special signature, as can be seen in a more clear way in Figure 5. The density chosen in these calculations can represent a solid surface plasma of aluminum irradiated by femtosecond lasers of the appropriate power (less than 1014 W/cm2). The aluminum layer does not have time to expand, and Z = 3 is produced by the conduction band electrons. This spectral information can be useful to properly understand radiation emitted from that surface, and can also help analyze compression conditions of inertially confined targets.
We will follow the standard definition of the ignition temperature as the value for which the fusion-born charged particle power equals the radiation loss power. A calculation has been carried out for the stoichiometric deuterium–tritium plasma, featured by the following fusion reactivity value (in m3/s):
In this formula, T is given in kiloelectron volts. Ions are not degenerate and have a Maxwellian velocity distribution.
In Figure 6, bremsstrahlung power and fusion-born alpha-particle power are plotted versus temperature for different density values. For the lower density (1025 electrons per cm3), the crossing between both lines (bremsstrahlung power and fusion power) belongs to the classical region, clearly identified by the slope of bremsstrahlung emission (as T1/2). On the contrary, for the highest electron density in the figure, 1029 cm−3, the crossing point corresponds to the degenerate case, and it is slightly below 1 keV, that is, somewhat lower than the standard ignition temperature, 4.5 keV. It is also seen from the picture that the electron density has to be over 1027 cm−3 for the degeneracy effects to be relevant to decrease the ignition temperature. For this density, the transition temperature between both bremsstrahlung regimes lies slightly below the classical ignition temperature, that is 4.5 keV.
It is worth pointing out that 1027 cm−3 is an extremely high electron density, more than 20,000 times DT solid density. In inertial fusion schemes based on self-ignited targets, such a compression factor seems to be out of reach, because it would need a very powerful and very symmetric spherical implosion with the formation of a central hot spark.
However, targets compressed for fast ignition could be degenerate because the main objective is to precompress the full target to very high densities at very low temperatures. During ignition, a complex interaction phase will develop, and a fraction of the ignited region will become hot and classical, but a significant fraction of the igniting region will actually start burning when it still is a degenerate plasma, because of the significant reduction in the radiation loss. Therefore, ignition requirements in the fast ignition concept could be significantly relaxed by the aforementioned reduction in the bremsstrahlung emission power. To produce degenerate plasmas in the compressed targets, fuel preheating must be avoided during the compression phase. Otherwise, such very high density values as those formerly cited in our calculations could not be attained.
Note that in the former analysis, ion and electron temperatures have been assumed to be equal. However, if a charged particle beam is used as an igniting trigger (Roth et al., 2001) the ion temperature in that region will be larger than the electron temperature, and this factor will also contribute to keep electron degeneracy for a longer while, thus reducing radiation losses. In fact, the interaction of charged particles with a degenerate plasma (Skupsky, 1977; Maynard & Deutsch, 1982) will become very important in this context and would have to be properly accounted for in the study of fast ignition in precompressed degenerate targets.
Degeneracy effects could also be of importance in the propagation of fusion burning waves in highly compressed targets (Chu, 1972; Eliezer & Martinez-Val, 1998) where a double front could develop, because electron heat conduction is not inhibited in degenerate plasmas. The former analysis of bremsstrahlung emission would therefore be relevant in physical cases as the one depicted in Figure 7, where a first wavefront of electron temperature moves ahead of the fusion burning front (at higher ion temperature). In highly compressed targets, the colder fuel could be degenerate and could reabsorb a fraction of the radiation emitted from the reacting zone. It is worth remembering that in the classical case, an optically thick region is given by the condition that the total radiation loss can not exceed the black body radiation limit:
where Ts is the surface temperature and Ta the average volume one, b being a term to characterize the bremsstrahlung emission (including the density factors) and σ the Stefan–Boltzmann constant. If the same temperature is assumed for Ta and Ts (which is an oversimplified hypothesis), it is found for the classical case that the spark radius for which the plasma becomes optically thick is
which is related to the fact that the radiation mean free path in a classical plasma scales as T7/2, assuming thermal equilibrium between electrons and radiation.
For the degenerate case, a similar analysis leads to
where b′ is the bremsstrahlung parameter for the degenerate formulation. Therefore, in degenerate plasmas, radiation mean free path scales as T2, and it is dominated by bound-free electron transitions in photon absorption, reciprocal to the emission process characterized by Eqs. (9) and (10). Radiation reabsorption could also help the propagation of a fusion burning wave across a plasma highly compressed up to degenerate levels (before the fusion onset), although the main effect will be a dramatic reduction in bremsstrahlung emission.
Bremsstrahlung is strongly inhibited in degenerate plasmas, because of the very broad range of electron energy transitions forbidden by Pauli's exclusion principle. The foregoing analysis has shown that actual values for radiation losses decrease very rapidly as the plasma becomes more and more degenerate. In the theoretical case of reaching ultrahigh compressed densities, the actual ignition temperature value would be lower than the classical one. It seems that density levels associated to this phenomenon are beyond the scope of inertial fusion implosions for self-triggered targets. However, bremsstrahlung emission reduction can be important in fast-ignited targets, which can have ignition requirements much lower than those predicted by classical bremsstrahlung models.