In this chapter, we investigate the stopping power of an ion in a magnetized electron plasma in a model of binary collisions (BCs) between ions and magnetized electrons, in which the two-body interaction is treated up to the second order as a perturbation to the helical motion of the electrons. This improved BC theory is uniformly valid for any strength of the magnetic field and is derived for two-body forces which are treated in Fourier space without specifying the interaction potential. The stopping power is explicitly calculated for a regularized and screened potential which is both of finite range and less singular than the Coulomb interaction at the origin. Closed expressions for the stopping power are derived for monoenergetic electrons, which are then folded with an isotropic Maxwell velocity distribution of the electrons. The accuracy and validity of the present model have been studied by comparisons with the classical trajectory Monte Carlo numerical simulations.
- ion stopping
- magnetized plasma target
- binary collisions
There is an ongoing in the theory of interaction of charged particle beams with plasmas. Although most theoretical works have reported on the energy loss of ions in a plasma without magnetic field, the strongly magnetized case has not yet received as much attention as the field-free case. The energy loss of ion beams and the related processes in magnetized plasmas are important in many areas of physics such as transport, heating, magnetic confinement of thermonuclear plasmas, and astrophysics. The range of the related topics includes ultracold plasmas [1, 2], the cooling of heavy ion beams by electrons [3, 4, 5, 6, 7, 8, 9, 10, 11, 12], as well as many very dense systems involved in magnetized target fusions , or heavy ion inertial confinement fusion (ICF).
For a theoretical description of the energy loss of ions in a plasma, there exist some standard approaches. The dielectric linear response (LR) treatment considers the ion as a perturbation of the target plasma, and the stopping is caused by the polarization of the surrounding medium. It is generally valid if the ion couples weakly to the target. Since the early 1960s, a number of calculations of the stopping power (SP) within LR treatment in a magnetized plasma have been presented (see Refs. [13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37] and references therein). Alternatively, the stopping is calculated as a result of the energy transfers in successive binary collisions (BCs) between the ion and the electrons [37, 38, 39, 40, 41, 42, 43, 44, 45]. Here, it is necessary to consider appropriate approximations for the screening of the Coulomb potential by the plasma . However, significant gaps between these approaches involve the ion stopping along magnetic field B and perpendicular to it. In particular, at high values, the BC predicts a vanishingly parallel energy loss, which remains at variance with the nonzero LR one. Also, challenging BCLR discrepancies persist in the transverse direction, especially for vanishingly small ion projectile velocity when the friction coefficient contains an anomalous term diverging logarithmically at [23, 24]. For calculation of the energy loss of an ion, two new alternative approaches have been recently suggested. One of these methods is specifically aimed at a low-velocity energy loss, which is expressed in terms of velocity-velocity correlation and, hence, to a diffusion coefficient . Next, in Ref.  using the Bhatnagar-Gross-Krook approach based on the Boltzmann-Poisson equations for a collisional and magnetized classical plasma, the energy loss of an ion is studied through a LR approach, which is constructed such that it conserves particle number locally.
An alternative approach, particularly in the absence of any relevant experimental data, is to test various theoretical methods against comprehensive numerical simulations. This can be achieved by a particle-in-cell (PIC) simulation of the underlying nonlinear Vlasov-Poisson Equation [10, 31]. While the LR requires cutoffs to exclude hard collisions of close particles, the collectivity of the excitation can be taken into account in both LR and PIC approaches. In the complementary BC treatment, the stopping force has been calculated numerically by scattering statistical ensembles of magnetized electrons from the ions in the classical trajectory Monte Carlo (CTMC) method [7, 10, 37, 38, 39, 40, 41]. For a review we refer to a recent monograph  which summarizes all theoretical and numerical methods and approaches also discussing the ranges of their validity.
The very recent upheaval of successful experiments involving hot and dense plasmas in the presence of kilotesla magnetic fields (e.g., at ILE (Osaka), CELIA (Bordeaux), LULI (Palaiseau), LLNL (Livermore)) remaining nearly steady during 10–15 ns strongly motivates the fusion as well as the warm dense matter (WDM) communities to investigate adequate diagnostics for their dynamic properties. This opens indeed a novel perspective by allowing magnetic fields to play a much larger if not a central role both in ICF and WDM plasmas. In this context proton or any nonrelativistic ion stopping is likely to provide an option of choice for investigating genuine magnetization features such as anisotropy, when the electron plasma frequency turns significantly lower than the cyclotron one . In addition, an experimental test of proton or alpha particle stopping in a magnetized plasma is currently envisioned (see, e.g., Ref.  for a preliminary discussion). The parameters at hand are a fully ionized hydrogen plasma with a density up to and temperature between 1 and 100 eV. The steady magnetic field can be up to 45 T strong. A preliminary examination based on comparing electron Debye length with corresponding Larmor radius indicates that to experience a strong influence of the magnetic field, the electron density should be comparable with a few . We expect these endeavors to lead to the very first unambiguous and genuine identification of an experimental magnetic signature for nonrelativistic ion stopping in plasmas.
Motivated by these recent developments, our purpose is to investigate the SP of an ion moving in a magnetized plasma in a wide range of the value of a steady magnetic field. The present paper is based on our earlier studies in Refs. [8, 24, 44, 45] where the second-order energy transfers for individual collisions of electron-ion [8, 24, 44] of any two identical particles, like electron-electron , and finally of two gyrating arbitrary charged particles  have been calculated with the help of an improved BC treatment. This treatment is—unlike earlier approaches of, e.g., Refs. [9, 42]—valid for any strength of the magnetic field. As the first application of the theoretical BC model developed in Refs. [8, 24, 44, 45], we have calculated in Ref.  the cooling forces on the heavy ion beam interacting with a strongly magnetized and temperature anisotropic electron beam. It has been shown that there is a quite good overall agreement with both the CTMC numerical simulations and the experiments performed at the ESR storage ring at GSI [48, 49, 50].
In Section 2 we introduce briefly a perturbative binary collision formulation in terms of the binary force acting between an ion and a magnetized electron and derive general expressions for the second-order (with respect to the interaction potential) stopping power. In contrast to the previous investigations in Refs. [8, 24, 44, 45], we here consider the (macroscopic) stopping force which is obtained by integrating the binary force of an individual electron-ion interaction with respect to the impact parameter and the velocity distribution function of electrons. That is, the stopping force for monoenergetic electrons is folded with a velocity distribution. The resulting expressions involve all cyclotron harmonics of the electrons’ helical motion and are valid for any interaction potential and any strength of the magnetic field. In Section 2.4 we present explicit analytic expressions of this second-order stopping power for the specific case of a regularized and screened interaction potential [51, 52] which is both of finite range and less singular than the Coulomb interaction at the origin and which includes as limiting cases the Debye (i.e., screened) and the Coulomb potentials. For comparison of our expressions with previous approaches, we consider in Section 3 the corresponding asymptotic expressions for large and small ion velocities and strong and vanishing magnetic fields. The analytical expressions presented in Section 2.4 are evaluated numerically in Section 4 using parameters of the envisaged experiments on ion stopping . In particular, we compare our approach with the CTMC simulations. The results are summarized and discussed in Section 5. The regularization parameter and the screening length involved in the interaction potential are briefly specified and discussed in Appendix A.
2. Theoretical model
2.1. Binary collision (BC) formulation
Let us consider two point charges with masses and charges , respectively, moving in a homogeneous magnetic field . We assume that the particles interact with the potential with , where is the permittivity of the vacuum and is the relative coordinate of the colliding particles. For two isolated charged particles, this interaction is given by the Coulomb potential, i.e., . In plasma applications is modified by many-body effects and the related screening and turns into an effective interaction. In general, this effective interaction, which is related to the wake field induced by a moving ion, is non-spherically symmetric and depends also on the ion velocity. For any BC treatment, however, this complicated ion-plasma interaction must be approximated by an effective two-particle interaction . This effective interaction may be modeled by a spherically symmetric Debye-like screened interaction with a screening length , given, e.g., by the Debye screening length (see, e.g., ), in case of low ion velocities and an effective velocity-dependent screening length for larger ion velocities (see [53, 54, 55]). Further details on the choice of the effective interaction are given in Ref. .
In the presence of an external magnetic field, the Lagrangian and the corresponding equations of particle motion cannot, in general, be separated into parts describing the relative motion and the motion of the center of mass (cm) . However, in the case of heavy ions, i.e., , the equations of motion can be simplified by treating the cm velocity as constant and equal to the ion velocity , i.e., . Then, introducing the velocity correction through relations , where is the relative electron-ion velocity and are the unperturbed electron and relative velocities, respectively, the equation of relative motion turns into
Here, is the force exerted by the ion on the electron, is the electron cyclotron frequency, and at . In Eq. (1) is the unit vector perpendicular to the magnetic field; the angle is the initial phase of the electron’s helical motion; is the relative velocity of the guiding center of the electrons, where and (with ) are the unperturbed components of the electron velocity parallel and perpendicular to , respectively; and is the cyclotron radius. In Eq. (1), the quantities and are defined by the initial conditions. In Eq. (2) is the ion-electron relative coordinate.
2.2. The perturbative treatment
We seek an approximate solution of Eq. (2) in which the interaction force between the ion and electron is considered as a perturbation. Thus, we are looking for a solution of Eq. (2) for the variables and in a perturbative manner , where are the unperturbed ion-electron relative coordinate and velocity, respectively, and are the nth-order perturbations of and , which are proportional to .
The parameter of smallness which justifies such kind of expansion can be read off from a dimensionless form of the equation of motion Eq. (2) by scaling lengths in units of the screening length , velocities in units of the initial relative velocity , and time in units of . Then, it is seen (see Ref.  for details) that the perturbative treatment is essentially applicable in cases where , that is, when the (initial) kinetic energy of relative motion , is large compared to the characteristic potential energy in a screened Coulomb potential. Or, expressed in velocities, the initial relative velocity v0 must exceed the characteristic velocity , that is, here demarcates the perturbative from the non-perturbative regime. If this condition is met not only for a single ion-electron collision but in the average over the electron distribution, e.g., by replacing v0 with the averaged initial ion-electron relative velocity , i.e., , we are in a regime of weak ion-target or, here, weak ion-electron coupling, which allows the use of perturbative treatments (besides BC also, e.g., linear response (LR)). For nonmagnetized electrons this is discussed in much detail in Refs. [53, 54]. Even though the particle trajectories are much more intricate in the presence of an external magnetic field, the given definitions and demarcations of coupling regimes are basically the same for magnetized electrons. That is, the applicability of a perturbative treatment is essentially related to the charge state of the ion and the typical range of the effective interaction, but not directly on the strength of the magnetic field. The latter may affect the critical velocity only implicitly via a possible change of the effective screening length with .
The equation for the first-order velocity correction is obtained from Eq. (2) replacing on the right-hand side of the exact relative coordinate by with the solutions and
Here, we have introduced the following abbreviations:
and have assumed that all corrections vanish at .
2.3. Second-order stopping power
We now consider the interaction process of an individual ion with a homogeneous electron plasma described by a velocity distribution function and a density . We assume that the ion experiences independent binary collisions (BCs) with the electrons. The total stopping force, , acting on the ion is then obtained by multiplying the binary force by the element of the flux relative flux , integrating with respect to time and folding with velocity distribution of the electrons. The impact parameter introduced here in the electron flux is defined by and is the component of perpendicular to the relative velocity vector with . As can be inferred from Eq. (1), represents the distance of the closest approach between the ion and the guiding center of the electron’s helical motion.
The resulting stopping power, , then reads
which is an exact relation for uncorrelated BCs of the ion with electrons. We evaluate this expression within a systematic perturbative treatment (see Ref.  for more details). First, we introduce the two-particle interaction potential , and the binary force is written using Fourier transformation in space. Furthermore, the factor in the Fourier transformed binary force is expanded in a perturbative manner as , where and are the unperturbed and the first-order corrected relative coordinates (Eqs. (1) and (3)), respectively. Next, we consider only the second-order binary force and the corresponding stopping force with respect to the binary interaction since the averaged first-order force (related to ) vanishes due to symmetry reasons [8, 24, 44, 45, 47]. Within the second-order perturbative treatment, the stopping power can be represented as
From Eq. (6) it is seen that the second-order stopping power is proportional to . Inserting now Eqs. (1) and (3) into Eq. (6), assuming an axially symmetric velocity distribution , and performing the s integration, we then obtain
where is the Bessel function of the th order; and are the components of parallel and transverse to , respectively; and and are the electron velocity components parallel and transverse to , respectively. This general expression (7) for the stopping power of an individual ion has been derived within second-order perturbation theory but without any restriction on the strength of the magnetic field .
2.4. The SP for a regularized and screened coulomb potential
For an electron plasma with an isotropic Maxwell distribution, the velocity distribution relevant for the averaging in Eq. (7) is given by
where the thermal velocity is related to the electron temperature by (here, the temperature is measured in energy units). Inserting Eq. (8) into expression (7) and assuming now a spherically symmetric potential yields after performing the velocity integrations (see Ref. ), the stopping power
Here, we have introduced the thermal cyclotron radius of the electrons , and , , where and are the ion velocity components transverse and parallel to , respectively. For the Coulomb interaction , the full two-dimensional integration over the s-space results in a logarithmic divergence of the integration in Eqs. (7) and (9). To cure this, cutoff parameters and must be introduced (see, e.g., Refs. [8, 24, 47] for details). These cutoffs are related to the screening of the interaction in a plasma target and the incorrect treatment of hard collisions in a classical perturbative approach. As an alternative implementation of this standard cutoff procedure, we here employ the regularized screened interaction with the Fourier transform
where represents a Debye-like screened interaction (see Section 2.1) which is additionally regularized at the origin [51, 52] and thus removes the problems related to the Coulomb singularity in a classical picture and prevents particles (for ) from falling into the center of the potential. The parameter related to this regularization is here considered as a given constant or as a function of the classical collision diameter .
Here, we have introduced the dimensionless quantities , . is the angle between and , , , , and
Eq. (11) for the SP is the main result of the outlined BC treatment which will now be evaluated in the next sections.
3. Comparison with previous approaches
Previous theoretical expressions for the stopping power which have been extensively discussed by the plasma physics community (see, e.g., Refs. [3, 8] for reviews) basically concern the two limiting cases of vanishing and infinitely strong magnetic fields. We therefore investigate the present approach for these two cases, first for arbitrary interactions and electron distributions as given by Eq. (7) and later for the more specific situation of the regularized interaction (10) and the velocity distribution (8) as given by Eq. (11).
3.1. General SP Eq. (7) at vanishing and infinitely strong magnetic fields
At vanishing magnetic field , and the argument of the Bessel function in Eq. (7) should be replaced by . Then, denoting the second-order SP at vanishing magnetic field as and assuming spherically symmetric potential with , one obtains
where is the generalized Coulomb logarithm:
Taking the bare Coulomb interaction with , Eq. (16) diverges logarithmically at and , and two cutoffs and must be introduced as discussed in Section 2.4. In this case the generalized Coulomb logarithm takes the standard form .
The asymptotic expression of Eq. (15) at high ion velocities can be easily derived using the normalization of the distribution function which results in
At an infinitely strong magnetic field , the term in Eq. (7) proportional to and the argument of the Bessel function vanish since the cyclotron radius . In this limit, denoting the SP as and assuming a spherically symmetric interaction potential, we arrive at
The corresponding high-velocity asymptotic expression is given by
Eqs. (15) and (19) and their asymptotic expressions for high velocities in Eqs. (18) and (20), respectively, agree with the results derived by Derbenev and Skrinsky in Ref.  in case of the Coulomb interaction potential, i.e., with . Using instead a regularized interaction potential and thus the Coulomb logarithm, allows closed analytic expressions and converging integrals and avoids any introduction of lower and upper cutoffs “by hand” in order to restrict the domains of integration. Moreover, employing the bare Coulomb interaction may, as pointed out by Parkhomchuk , result in asymptotic expressions which essentially different from Eqs. (19) and (20), which is related to the divergent nature of the bare Coulomb interaction (see Ref. ).
3.2. Some limiting cases of Eq. (11)
Next, we discuss some asymptotic regimes of the SP (Eq. (11)) where the regularized interaction (Eq. (10)) and the isotropic velocity distribution (Eq. (8)) have been assumed. In the high-velocity limit where , only small t contributes to the SP (Eq. (11)) due to the short time response of the electrons to the moving fast ion. In this limit we have . The remaining integration can be performed explicitly. This integral is given by .
Here, the function is determined by Eq. (14), and is the generalized Coulomb logarithm (Eq. (17)). The remaining expressions do not depend on the magnetic field, i.e., , as a natural consequence of the short time response of the magnetized electrons. In fact, and and the related integration (Eq. (21)) are also valid for vanishing magnetic field . Integration by parts turns Eq. (11) into
where is the error function and is again the scaled ion velocity. The SP (Eq. (22)) is isotropic with respect to the ion velocity and represents the two limiting cases of high velocities at arbitrary magnetic field and arbitrary velocities at vanishing field. Of course, expression (22) can be also obtained by performing the remaining integration in the nonmagnetized SP (Eq. (15)) using the isotropic velocity distribution (Eq. (8)) and .
A further increase of the ion velocity finally yields
which completely agrees with the asymptotic expression (18) in case of . Inspecting Eq. (23) shows that the SP does not depend explicitly on the electron temperature at sufficiently high velocities, while may still be involved in the generalized Coulomb logarithm .
At and small velocities , the SP (Eq. (22)) becomes
Now, we consider the situation when the magnetic field is very strong and the electron cyclotron radius is the smallest length scale, , and the SP is only weakly sensitive to the transverse electron velocities and, hence, is affected only by their longitudinal velocity spread. In this limit and are obtained from Eq. (11) after straightforward calculations:
Up to the definition of the Coulomb logarithm (i.e., versus ), the expressions are identical to those obtained by Pestrikov .
respectively, where (with ) is the modified Bessel function. It is also constructive to obtain the angular averaged stopping power. From Eq. (25) one finds
where is the exponential integral function.
In the high-velocity limit with , the SP (Eq. (25)) becomes
The SP given by Eq. (30) (or Eq. (20) with ) decays as the corresponding SP (Eq. (23)) like with the ion velocity. But here, the parallel SP (Eq. (27)) vanishes exponentially at which is a consequence of the presence of a strong magnetic field, where the electrons move parallel to the magnetic field. If the ion moves also parallel to the field (i.e., ), the averaged stopping force must vanish within the BC treatment for symmetry reasons.
Finally, we also investigate the case of small velocities at strong magnetic fields. Considering a small ion velocity in Eq. (25), we arrive at
where is Euler’s constant. Now, it is seen that the SP, , leads at low ion velocities and for a nonzero to a term which behaves as . Thus, the corresponding friction coefficient diverges logarithmically at small . This is a quite unexpected behavior compared to the well-known linear velocity dependence without magnetic field (see asymptotic expressions above). Finally, at the logarithmic term vanishes and the SP behaves as .
4. Features of the SP (Eq. (11)) and comparison with CTMC simulations
In this section we study some general properties of the SP of individual ions resulting from the BC approach by evaluating Eq. (11) numerically. We consider the effect of the magnetic field on the SP at various temperatures of the plasma. The density and the temperatures , 10 or 100 eV of the electron plasma, are in the expected range of the envisaged experiments on proton or alpha particles stopping in a magnetized target plasma  (see corresponding Figures 1–3). As an example we choose proton projectile for our calculations. In all examples considered below, the regularization parameter thereby meets the condition , i.e., , and does not affect noticeably the SP (Eq. (11)) at low and medium velocities as shown in Appendix A (see also Ref.  for more details).
For a BC description beyond the perturbative regime, a fully numerical treatment is required. In the present cases of interest, such a numerical evaluation of the SP is rather intricate but can be successfully implemented by classical trajectory Monte Carlo (CTMC) simulations [37, 38, 39, 40]. In the CTMC method, the trajectories for the ion-electron relative motion are calculated by a numerical integration of the equations of motion (Eq. (2)). The stopping force is then deduced by averaging over a large number (typically –) of trajectories employing a Monte Carlo sampling for the related initial conditions. For a more detailed description of the method, we refer to Refs. [8, 44, 45]. Both the analytic perturbative treatment and the non-perturbative numerical CTMC simulations are based on the same BC picture and use the same effective spherical screened interaction . The following comparison of these both approaches thus essentially intends to check the validity and range of applicability of the perturbative approach as it has been outlined in the preceding sections.
5. Stopping profiles and ranges
5.1. General trends
The parameter analysis initiated on Figures 1–3 at and is implemented for monitoring a possible experimental vindication through a fully ionized hydrogen plasma out of high-power laser beams available on facilities such as ELFIE (Ecole Polytechnique) or TITAN (Lawrence Livermore) . The given adequately magnetized targets (in the 20–45 T range) would then be exposed to TNSA laser-produced proton beams out of the same facilities, in the hundred keV-MeV energy range .
Therefore, we are looking for the most conspicuous effect of the applied magnetized intensity B on the proton stopping.
Moreover, that anisotropy evolves only moderately between and .
Another significant feature is the extension to any of the scaling . For instance, SP at and , at a given , is equivalent to that for and .
As expected, B effects impact essentially the low-velocity section (, = target electron thermal velocity) of the ion stopping profile. One can observe, increasing with B, a shift to the left of SP maxima, as shown in Figure 4 at tesla, for the profiles displayed in Figure 3, with -averaged SP remaining close to .
Switching now attention to corresponding ranges, down to projectile at rest , one witnesses on Figure 5 the counterpart of the above-noticed SP behavior.
In a low projectile velocity , one gets the largest B effects and the smallest proton ranges attributed to the highest B. The fan of B ranges then merges on a given point, located between 10 keV and 100 keV at , and then inverts itself with increasing B featuring now increasing ranges. Moreover, the aperture of the fan of ranges increases steadily with .
5.2. Specific trends
The projected experimental setup  could manage constant, static, and homogeneous B values up to 45 T. So, we are let to investigate range limits within which significant B effects can be observed.
Obviously, is expected to show quantitatively larger B impact than .
Giving attention to proton ranges of T dependence in a low-density plasma at and 100 eV, respectively, (Figure 6), one witnesses the smallest ranges for , increased by four orders of magnitude between 1 and 100 eV while remaining essentially unchanged for . Turning now to at , one can see that the given SP remains quasi-isotropic, hardly -dependent, except at extreme magnetization (). Discrepancies between B = 0 and 20 T remain visible only for . does not feature anymore the highest stopping when . Also, SP exhibits a few top wigglings. Upshifting T at 10 eV yields back SPs very similar to these displayed on Figure 2 (, ) Figure 7.
Corresponding proton ranges (, ) are shown in Figure 8.
Experimentally, accessible and very small ranges are thus documented for . Here, and 20 T data remain everywhere distinguishable.
5.3. Very-low-velocity proton slowing down
Up to now we limited our investigation to proton stopping by target electrons. In the very-low-velocity regime , the target protons can also contribute significantly as evidenced on Figure 9. This topic will be more thoroughly addressed in a separate presentation.
We developed and extensively used a kinetic approach based on a binary collision formulation and suitably regularized Coulomb interaction, to numerically document for any value of the applied magnetization B, the stopping of a proton projectile in a fully ionized hydrogen plasma target. Both ion projectile and target plasma parameters have been selected in order to fit a planned ion-plasma interaction experiment in the presence of an applied magnetic field . It should be pointed out that we restricted the target plasma to its electron component. It therefore remains to include the target ion contribution to proton stopping , thus featuring a complete low-velocity ion slowing down.
More generally, we expect that the present investigation, experimentally geared as it is, could help to bridge a long-standing and persisting gap between theoretical speculations and experimental facts in the field of nonrelativistic ion stopping in magnetized target plasmas.
The work of H.B.N. has been supported by the State Committee of Science of the Armenian Ministry of Education and Science (Project No. 13-1C200). This work was supported by the Bundesministerium für Bildung und Forschung (BMBF) under contract 06ER9064.
Our results (Eq. (11)) were derived by using the screened interaction . As already mentioned, the use and the modeling of such an effective two-body interaction are a major but indispensable approximation for a BC treatment where the full ion-target interaction is replaced by an accumulation of isolated ion-electron collisions. The replacement of the complicated real non-spherically symmetric potential, like the wake fields as shown and discussed in Ref. , with a spherically symmetric one is, however, well motivated by earlier studies on a BC treatment at vanishing magnetic field (see Refs. [53, 54, 55]). It was shown by comparison with 3D self-consistent PIC simulations that the drag force from the real nonsymmetric potential induced by the moving ion can be well approximated by an BC treatment employing a symmetric Debye-like potential with an effective velocity-dependent screening length . In these studies also a recipe was given how to derive the explicit form of , which turned out to be not too much different from a dynamic screening length of the simple form . Here, is the Debye screening length at , is the electron plasma frequency, and is a thermal velocity of electrons. Although no systematic studies about the use of such an effective interaction with a screening length have been made for ion stopping in a magnetized electron plasma, the replacement of the real interaction by a velocity-dependent spherical one should be a reasonable approximation also in this case. The introduced dynamical screening length also implies the assumption of a weak perturbation of the electrons by the ion and linear screening where the screening length is independent of the ion charge , which coincide with the regimes of perturbative BC (see, e.g., Ref. ). Therefore, we do not consider here possible nonlinear screening effects.
Next, we specify the parameter which is a measure of the softening of the interaction potential at short distances. As we discussed in the preceding sections, the regularization of the potential (Eq. (10)) guarantees the existence of the s integrations, but there remains the problem of treating accurately hard collisions. For a perturbative treatment, the change in relative velocity of the particles must be small compared to , and this condition is increasingly difficult to fulfill in the regime . This suggests to enhance the softening of the potential near the origin of the smaller . Within the present perturbative treatment, we employ a dynamical regularization parameter [44, 45], where and . Here, is the averaged distance of the closest approach of two charged particles in the absence of a magnetic field, and is some free parameter. In addition we also introduced in . In Refs. [44, 45], this parameter is deduced from the comparison of the second-order scattering cross sections with an exact asymptotic expression derived in Ref.  for the Yukawa type (i.e., with ) interaction potential. As we have shown in Refs. [44, 45] employing the dynamical parameter , the second-order cross sections for electron-electron and electron-ion collisions excellently agree with CTMC simulations at high velocities. Also, the free parameter is chosen such that , where is the distance at . From the definition of , it can be directly inferred that does not play any role at low velocities, while it somewhat affects the size of the stopping force at high velocities when . More details on the parameter and its influence on the cooling force are discussed in Ref. .