The Effect of Nuclear Elastic Scattering on Temperature Equilibration Rate of Ions in Fusion Plasma

A plasma with two different particle types and at different temperatures has been considered, so that each type of ion withMaxwellBoltzmann distribution function is in temperature equilibrium with itself. Using the extracted nuclear elastic scattering differential cross-section from experimental data, solving the Boltzmann equation, and also taking into account themobility of the background particles, temperature equilibration rate between two different ions in a fusion plasma is calculated.The results show that, at higher temperature differences, effect of nuclear elastic scattering is more important in calculating the temperature equilibration rate. The obtained expressions have general form so that they are applicable to each type of particle for background (b) and each type for projectile (p). In this paper, for example, an equimolar Deuterium-Hydrogen plasma with density n = 5 × 1025 cm is chosen in which the deuteron is the background particle with temperature (also electron temperature) T b = 1 keV (usual conditions for a fusion plasma at the ignition instant) and the proton is the projectile with temperature T p > T b . These calculations, particularly, are very important for ion fast ignition in inertial confinement fusion concept.


Introduction
In a fusion plasma, because of interactions, collisions, and gain and loss powers, we encounter an unstable plasma [1,2].The thermal equilibrium between charged particles is crucial for understanding the overall energy balance in a fusion plasma, where the ignition and burn of the plasma are strongly temperature dependent [3].Since a plasma is composed of charged particles, they are affected by Coulomb force in the plasma.Coulomb force is a long-range force and, thus, may also be involved in the short-range collisions and the long-range collective effects.The problem of temperature equilibration was studied first by Landau [4,5] and Spitzer (LS) [6].The LS model is applicable to dilute, hot, fully ionized plasmas where the collisions are weak and binary.Recently, there have been several theoretical and computational works which aim to study temperature equilibration [7][8][9][10][11][12].Among the theoretical studies, Brown, Preston, and Singleton (BPS) by solving the Boltzmann equation for short-range collisions and the Lenard-Balescu equation for long-range collective effects and also considering the quantum corrections have achieved more accurate results [9].In a hot and dense fusion plasma, in addition to the Coulomb interactions, the nuclear elastic scattering (NES) is able to role via short-range interactions (hard collisions).The effect of the NES on the stopping power of a charged particle has been investigated in previous studies [13][14][15].It has been indicated that the stopping power due to NES is greater at the higher projectile energies.This fact motivated us to study the effect of NES on thermal equilibration rate for a fusion plasma.A fusion plasma is formed from charged particles at high temperatures in which they interact with each other to achieve balance.In this paper, for the simplest case, a plasma with two different ions at different temperatures has been considered, so that each type of ion with Maxwell-Boltzmann distribution function is in temperature equilibrium with itself.Here, for example, proton at higher temperature is considered as projectile and deuteron at lower temperature is assumed as background.By solving the Boltzmann equation and taking into account the thermal motion of plasma ions, the thermal 2 Advances in High Energy Physics equilibration rate due to NES is calculated.The results are compared with Coulomb calculations and the conditions in which each contribution is dominated are discussed.

Elastic Collision
Consider a projectile with mass   and velocity k  interacting with a background particle with mass   and velocity k  in laboratory frame.(Vectors are written in bold type and nonbold letters indicate quantity magnitude.)In an elastic collision, V  = (  k  +   k  )/  is the velocity of the center of mass (C.M.) and relative velocity of the system is k  = k  − k  .After the collision, the projectile and target velocities consider k   and k   , respectively.Therefore, k   = k   − k   is the relative velocity of two particles after collision.These quantities are related to each other by where   =   +   is the total mass.The conservation laws for linear momentum and kinetic energy in the C.M. system indicate that the relative speeds of the particles do not change before and after collision as a result of an elastic collision |k  | = |k   | = V  .As a result, the magnitude of the relative velocities remains unchanged and only the angle Θ (scattering angle in C.M. system) between k  and k   is changed.Since the magnitude and direction of C.M. velocities (before and after a collision) will not change, V  can be selected as the reference.Using the above relations between V  and k  velocities, the kinetic energy of the projectile is changed as follows: where   =     /  is reduced mass of the projectile and target system.If  is the angle between V  and k  , according to scattering angle in C.M. system Θ, then angle between V  and k   is ( − Θ).The change in kinetic energy of the projectile is equal to Since the scattering in the C.M. frame is axially symmetric about the relative speed k  , transverse components average is zero in the scattering process, and so the second term on the right of this equation will be removed.

Temperature Equilibration Rate due to Coulomb Interaction
The problem of temperature equilibration was first addressed by Landau and Spitzer, who used the Fokker-Planck equation to derive an equilibration rate for one ion species temperature   , given another ion species temperature   [4][5][6]: where  LS  is defined as where  , are the ion species masses,  , are the ion charges, and  , are the number densities.The LS Coulomb logarithm ln Λ LS  is defined as follows: Recently, there have been several theoretical and computational works which aim to study temperature equilibration [7][8][9].The computational studies have all been done with classical molecular dynamics and have focused on temperature equilibration in a multi-eV hydrogen plasma [10][11][12].The results of these studies showed that deviations from LS approach are to be expected, even for moderate values of ln Λ LS  .Among the theoretical works, Brown, Preston, and Singleton (BPS) produced an analytic calculation for Coulomb energy exchange processes for a fusion plasma [9,16,17].These precise calculations are accurate to leading and next-to-leading order in the plasma coupling parameter and to all orders for two-body quantum scattering within the plasma.In general, the energy density exchange rate between the two charged particles  and  in a plasma can be written as follows: where  Coul  is called the rate coefficient and it can be written as a sum of three terms: where the first two arise from classical short and long distance physics and the latter term arises from short distance twobody quantum diffraction.The first term in (8), the shortdistance classical scattering contribution, reads where  ≃ 0.57721 . . . is Euler constant.The second term in (8), the long-distance, dielectric term that accounts for collective effects in the plasma, is given by Here,  is an arbitrary wave number so that the total result does not depend upon .However, sometimes choosing  to be a suitable multiple of the Debye wave number of the plasma simplifies the formula.We write the inverse temperature of the projectile as   =  −1  and the plasma species  as   =  −1   , which we measure in energy units.Debye wave number   of this species is defined by where   is the number density of species .The total Debye wave number   is defined by The function () is related to the leading-order plasma dielectric susceptibility in which it may be expressed in the dispersion form: where the limit  → 0 + is understood.The spectral weight,  total (V) is defined by where According to the BPS calculations, the general case is obtained by adding a quantum correction, which is related to the short distance two-body quantum diffraction and to the classical result.Therefore, the third term in (8),    is a coefficient which is calculated from the quantum correction and is defined as follows: the strength of the quantum effects associated with the scattering of two plasma species  and  is characterized by the dimensionless parameter: where the square of the thermal velocity in this expression is defined by and () is the logarithmic derivative of the gamma function.

Temperature Equilibrium due to Nuclear Elastic Scattering
Since nuclear force is a short-range force, it can only be involved in short-range collisions and, therefore, does not contribute to long-range cumulative effects of plasma.Boltzmann equations are described generally as short-range encounters.The equation for the phase-space density   (P  ) We suppress the common space and time coordinates r,  and write the collision term involving species  in the following form: This relation describes the scattering of the particles of masses   and   and the scattering from the initial momenta p  =   k  , p  =   k  to the final momenta p   =   k   , p   =   k   with the scattering amplitude (,  2 ) depending on the center-of-mass energy, , and the squared momentum transfer,  2 .It is convenient to employ this quantummechanical notation for several reasons.It explicitly displays the complete kinematical character of a scattering process, including the detailed balance symmetry.Furthermore, it shows that the collision term (19) vanishes when all the particles are in thermal equilibrium with the generic densities (p) ∼ exp[−V 2 /2] because of the conservation of energy enforced by the delta function.In case we have a spatial uniformity (homogeneity), include the gradient removed Advances in High Energy Physics from (19) and, thus, the rate of change of energy is obtained as follows: Placement of relation (20) in (21) and rate of change of energy density of species  in collisions with plasma species  will be written as follows: Now, in general, the cross section for the scattering of particles  and  into a restricted momentum interval Δ is given by We can write (23) in terms of the speed of the center of mass V  and relative speed k  .We use the conversion ) is the Jacobian.Using these definitions, we find that Equation ( 24) represents the rate of change of energy density of projectile particles  in collisions with plasma species .Now, consider the particles distribution functions  and ,   (p  ), and   (p  ) where is the Maxwellian distribution function.Then temperature Equilibration rate can be calculated as follows; where the coefficient  NI  is given by where the integral and   = 1/2  V 2  is the total energy of the projectile and the target is in the center of mass system.Extracting of NES differential cross-section from experimental values, ( NI  /Ω) rate of temperature equilibration is archived.It is seen that ( 26) when the two types of ion temperature are equal,  NI  / = 0 which indicates the temperature equilibrium.

Results and Discussion
Now, we will consider the thermal equilibration rate between the proton and deuteron in an equimolar plasma with density  = 5 × 10 25 cm −3 .The deuteron is considered as the background at temperature   = 1keV and the proton is the projectile at temperature   , greater than the background temperature (  >   ).It is assumed that each type of particle is in equilibrium with itself so that the equilibrium is described by Maxwell-Boltzmann distribution function.The density and temperatures are chosen so that the their values are close to initial conditions of an inertial confinement fusion fuel at ignition moment.The thermal equilibration rate due to Coulomb interactions is obtained from ( 5) and ( 7) for LS and BPS relations, respectively (Figure 1).The portions of singular term (see (9)), regular term (see (10)), and quantum term (see ( 16)) of BPS relation are shown, separately.It is seen that the values obtained from the BPS calculations would give lower values than LS calculations.However, the most important point is that both equilibration rates are further reduction at the higher temperature differences.This event is predictable of course, because the Coulomb differential cross section (Rutherford's differential cross section) is proportional with the total energy in center of mass system (  ) as  Coul /Ω ∝  −2   .The total energy in the center of mass system increases with increasing the temperature difference between the two ions, and, thus,  4) and ( 7) for LS and BPS relations, respectively.the rate of Coulomb interaction decreases due to reduction of its cross-section.
In Figure 2, the ratio of the experimental differential cross-section ( exp  /Ω) to the Coulomb differential crosssection ( Coul  /Ω) is plotted versus the scattering angle in the center of mass system (Θ) for different energies, for proton-deuteron elastic scattering.The labeled energies in the figure are the total energy in the center of mass system (  ).As can be seen, this ratio is different for different energies and also is associated with fluctuations.However, in general, this ratio is increased more by increasing the energy and this means that the experimental differential crosssection is more deviated from the Coulomb differential crosssection.Deviation increase reflects the growth of the contribution of NES differential cross-section ( NI  /Ω) at higher energies (the superscript NI denotes the sum of nuclear and interference terms that, for simplicity, are called nuclear elastic scattering (NES)).
By extracting the NES differential cross-section from experimental data ( NI  /Ω =  exp  /Ω −  Coul  /Ω), the integral  NI (  ) (see (28)) can be calculated (Figure 3).The integral  NI (  ) is submitted in (27) and the thermal equilibration rate due to NES can be obtained.Figure 4 shows the thermal equilibration rate due to NES (NI) and Coulomb interaction (LS and BPS) and also the total term that is obtained from sum of nuclear and Coulomb contributions (NI + BPS).It is seen that the thermal equilibration rate due to Coulomb interaction has the main contribution in the lower temperature differences.With increasing temperature, the Coulomb interaction contribution decreases and increases in the share of NES.Thus, at temperature differences, Δ ≥ 3.278 MeV; the NES is the dominant contribution to the thermal equilibration rate.
The rapid development of laser-accelerated ion beams is developing the ion fast ignition (IFI) scheme in inertial confinement fusion concept [19].In the IFI scheme, an intense laser pulse accelerates an ion beam to multi-MeV temperatures.The accelerated ion beam deposits its energy in the ignition region of fuel pellet that is named "hot spot" and provides the required energy for ignition [20,21].For this reason, the diver ion beam is named ignitor.The power is deposited through an ion beam energy exchange with the background plasma particles with few-keV temperatures.In the energy exchange process, some of the ion beam energy is delivered to the plasma ions, and another part of energy is delivered to the plasma electrons.The electrons are the source of power loss in a fusion plasma that it is better to be kept at low temperature.At the high temperature difference (multi-MeV) between ion beam and background plasma particles, the rate of energy exchange via Coulomb interaction is reduced and, in turn, the rate of energy exchange via NES is increased.The NES causes the ion temperature to be further increased, resulting in reduced power loss.The laser-accelerated ion beams can be produced with different distribution functions, such as exponential, Maxwellian, and quasimonoenergy distribution functions [22,23].Furthermore, some of ion beams, such as the deuteron beam can be on their way to interact with the background ions [24,25].Therefore, research in this area requires more detailed calculations that will be covered in the future.

1 )Figure 1 :
Figure 1: The thermal equilibration rate due to Coulomb interactions is obtained from (4) and (7) for LS and BPS relations, respectively.

Figure 2 :
Figure 2: The ratio of the experimental differential cross-section ( exp  /Ω) to the Coulomb differential cross-section ( Coul  /Ω) is plotted versus the scattering angle in the center of mass system (Θ) for different energies   , for proton-deuteron elastic scattering.

Figure 3 :
Figure 3: The integral  NI (  ) (see (28)) versus the total energies in the center of mass coordinate   .

1 )Figure 4 :
Figure 4: The thermal equilibration rate due to NES (NI) and Coulomb interaction (LS and BPS) and also the total term that is obtained from the sum of nuclear and Coulomb contributions (NI + BPS).