Microbunched Electron Cooling with Amplification Cascades

The Microbunched Electron Cooling (MBEC) is a promising cooling technique that can find applications in future hadron and electron-ion colliders to counteract intrabeam scattering that limits the maximum achievable luminosity of the collider. To minimize the cooling time, one would use amplification cascades consisting of a drift section followed by a magnetic chicane. In this paper, we first derive and optimize the gain factor in an amplification section for a simplified one-dimensional model of the beam. We then deduce the cooling rate of a system with one and two amplification cascades. We also analyze the noise effects that counteract the cooling process through the energy diffusion in the hadron beam. Our analytical formulas are confirmed by numerical simulations for a set of model parameters.


I. INTRODUCTION
Microbunched coherent electron cooling (MBEC) of relativistic hadron beams has been proposed by D. Ratner [1] as a way to achieve cooling rates higher than those provided by the coherent cooling using a free electron laser [2]. The mechanism of MBEC can be understood in a simple setup shown in Fig. 1. An electron beam with the same relativistic γ-factor as a) Amplification section in the hadron beam, co-propagates with the hadrons in a section of length L m called the "modulator". In this section, the hadrons imprint microscopic energy perturbations onto the electrons via the Coulomb interaction. After the modulation, the electron beam passes through a dispersive chicane section, R (e,1) 56 , where the energy modulation of the electrons is transformed into a density fluctuation referred to as "microbunching". This chicane is followed by an amplification section consisting of a drift of length L d and another chicane R (e,2) 56 , as shown in Fig. 1a. If the length of the drift is equal to one-quarter of the plasma oscillation period in the electron beam, 1 4 λ p , and the chicane strength is properly optimized, the density fluctuations in the electron beam generated by the chicane R Theoretical analysis of MBEC without amplification has been carried out in a recent study [3]. In this paper, we extend the analysis of Ref. [3] to include the amplification sections. Following the approach developed in Ref. [3] we adopt a general framework in which we look at the dynamics of the fluctuations in both beams. We assume that before the beams start to interact, their density and energy fluctuations can be described as uncorrelated shot noise. In the process of interaction, the fluctuations in the electron and hadron beams establish correlations, and when the beams are recombined in the kicker the fluctuating electric field in the electron beam acts in a way that decreases the energy spread in the hadron beam. As in Ref. [3], for the hadron-electron, as well as electron-electron interactions we adopt a model in which the particles are replaced by thin disks with a Gaussian transverse charge distribution.
The paper is organized as follows. In Sec. II we summarize the Coulomb interaction between thin slices with transverse Gaussian distribution of charge. In Sec. III we study plasma oscillations in a beam consisting of thin Gaussian slices. In Sec. IV we demonstrate that an initial sinusoidal modulation of small amplitude in a beam is amplified after the passage through a quarter of plasma wavelength drift and a subsequent chicane. The amplification factor derived in this section is then used, in Secs. V and VI, for calculation of the cooling rate in an MBEC cooling system with one and two amplification sections, respectively. In Sec. VII we discuss the wake field associated with the amplified cooling which is related to the effective energy exchange of two hadrons located at a given distance z. In Sec. VIII we present results of computer simulations of the cooling rates. In Sec. IX the noise and saturation effects in the cooling process are studied, and in Sec. X numerical estimates of the hadron cooling are presented for the eRHIC electron-ion collider design. We concluded this paper with the summary in Sec. XI.
We use the Gaussian system of units throughout this paper.

II. INTERACTION OF CHARGED GAUSSIAN SLICES
As was already mentioned in the Introduction, we treat the Coulomb interaction between particles as if a hadron were a disk of charge Ze with an axisymmetric Gaussian radial distribution with the rms transverse size equal to the rms transverse size of the beam. The electron is also modeled by a Gaussian disk of charge −e with the same transverse profile. A similar Gaussian-to-Gaussian interaction model was used in 1D simulations of a longitudinal space charge amplifier in Ref. [4].
In this model, a hadron of charge Ze at the origin of the coordinate system exerts a force f z on an electron at coordinate z, where Σ is the rms beam radius and the function Φ is defined by the following expression [5], with erfc the complementary error function. The function Φ is odd, Φ(−x) = −Φ(x); its plot can be found in Ref. [3]. Neglecting the relative longitudinal displacements of hadrons and electrons in the modulator, the force (1) causes the relative energy change ∆η of an electron located at coordinate z, where L m is the length of the modulator and r e = e 2 /m e c 2 is the classical electron radius, and we use the notation η for the energy deviation ∆E of a particle normalized by the nominal beam energy γmc 2 , η = ∆E/γmc 2 . Eq. (3) can also be considered as a Green function for the energy modulation of electrons induced by a delta-function density perturbation in the hadron beam.
In our analysis we will assume that the beam radius in the amplification sections, Σ p , may be different from the beam radius in the kicker and the modulator, Σ. For the electronelectron interaction in these sections we use Eq. (1) with Z = −1 and Σ → Σ p , In what follows, we will also need the Fourier transform of the function Φ. Because of the antisymmetry of the function Φ its Fourier transform is purely imaginary, so we define function H as The plot of function H(κ) is shown in Fig. 2. For large values of the argument, κ 1, the function H asymptotically approaches 1/2κ. One can also find an approximation for H near the origin, κ 1, H ≈ κ(2 ln κ − γ E ), where γ E = 0.577 is the Euler constant.

III. PLASMA OSCILLATIONS IN A GAUSSIAN BEAM
To calculate the increase in the amplitude of the fluctuations in the electron beam when it propagates through the amplification sections, we first need to analyze the beam plasma oscillations in the drift. A similar problem has been studied in Refs. [6], however, our treatment is simpler because we will assume a cold plasma and neglect the transverse degrees of freedom in the beam. In our analysis, we use the Vlasov equation for the distribution function in the longitudinal phase space, f (z, η, t), normalized so that ∞ −∞ f dη = n, with n the number of particles in the beam per unit length. Here z = s − vt is the longitudinal coordinate in the beam with s the distance measured along the direction of the beam propagation in the lab frame, and v the nominal beam velocity. We represent the distribution function as f = n 0 F 0 (η) + δf (z, η, t), where F 0 (η) is the equilibrium beam distribution function, n 0 is the nominal linear particle density, and δf describes small-amplitude time-dependent fluctuations in the beam, |δf | F 0 . We consider the fluctuations with the longitudinal scale much smaller than the bunch length and carry out our analysis in a small vicinity of a given location in the bunch where the variation of the distribution function F 0 with z can be neglected; for this reason the coordinate z is omitted from the arguments of the function The linearized Vlasov equation for the perturbation of the distribution function, δf , is whereη is the energy change per unit time. The rate of energy change is expressed through the longitudinal force in the electron beam, where and F z is given by Eq. (4). Making the Fourier transform of Eqs. (6) and (7) and using the notations we obtain with the effective impedance ζ(k) given by where the function H is defined by Eq. (5). At large values of k, k γ/Σ p , we have H(κ) ∼ 1/κ, so the impedance in this region can be estimated as ζ ∼ e 2 /Σ 2 p γmck.
In our analysis we will assume that the second term in Eq. (10) is much smaller than the third one. The conditions for such assumption are estimated at the end of this section.
Neglecting the second term we can integrate the Vlasov equation over time, To get an equation for δn k (t) we integrate Eq. (10) over η, where is the averaged perturbation of the energy η by δf k . Note that the large third term in Eq. (10) does not contribute to this equation, so we have to keep the contribution from the second term. We can also obtain an equation for δq k by integrating Eq. (10) with weight η, where we have neglected the contribution from the second term in Eq. (10). Combining Eqs. (13) and (15) we find which is the equation for plasma oscillations in the beam. A somewhat different derivation of plasma oscillations in a relativistic beam is given in Ref. [7]. Note that there is no Landau damping effects [6] in this equation, which means that our assumption of the smallness of the second term in Eq. (10) is equivalent to the cold plasma approximation. It follows from Eq. (16) that the plasma frequency ω p (k) is given by the following equation: with and κ p = kΣ p /γ. The plot of the function 2κ p H(κ p ) -this function is equal to the normalized plasma frequency ω p /Ω -is shown in Fig. 3. One can see that the short-period plasma oscillations with k γ/Σ p , the plasma frequency is approximately equal to Ω p , while for the long-wavelength oscillations ω p decreases with the wavelength.
We can now find the condition when the second term in Eq. (10) is much smaller than the third one. In our estimates we will assume k ∼ γ/Σ p which gives for the impedance ζ ∼ e 2 /Σ p γ 2 mc. Using the estimates F 0 ∼ 1/σ 2 e , δf k ∼ δn k /σ e and η ∼ σ e we find that the ratio of the third term to the second one is approximately equal to which has to be much greater than one. Here I e = ecn 0 is the peak beam current and I A = mc 3 /e is the Alfvén current. In the next section, we will see that the same parameter (19) appears in the expression for the gain factor of the amplification section.

IV. GAIN FACTOR IN AN AMPLIFICATION CASCADE
The solution of Eqs. (13) and (15) is where δn k (0) and δq k (0) are the initial values of the density and energy perturbations in the beam. Let us compare the first and the second terms in this equation taking into account that plasma oscillations occur after the fist chicane, R (e,1) 56 , of the cooling section. In these estimates we again assume k ∼ γ/Σ p and the optimal value for the chicane strength, R (e,1) 56 ∼ Σ p /γσ e (see order of magnitude estimates in Ref. [8]). The magnitude of δn k (0) is estimated as δn k (0) ∼ kR (e,1) 56 δq k (0) ∼ δq k (0)kΣ p /γσ e , so that the ratio of the second term to the first one in Eq. (20) is of the order of cσ e /γΩΣ p . This combination of parameters turns out to be equal to the inverse of the square root of the parameter in Eq. (19), and, by assumptions, is much smaller than one. Hence we can neglect the second term in Eq. (20), Substituting this result in Eq. (12) we obtain Estimating the relative magnitude of the two terms on the right-hand side in this equation, as it was done in the derivation of Eq. (19), we find that the ratio of the first term to the second one is given by the same parameter (21) and hence we can neglect the first term in Eq. (23). Using this expression for δf k , we can find the linear density perturbation in the beam, δn δn where we have replaced the time by the length of the drift divided by the speed of light, The last expression in this formula is calculated for a Gaussian distribution function F 0 = (2π) −1/2 σ −1 e e −η 2 /2σ 2 e with σ e the rms relative energy spread in the electron beam. It makes sense to define the gain factor G of the amplification section as G = δn (2) k /δn k (0). Using Eq. (11) for the impedance and Eq. (17) for the plasma frequency, after simple calculations, we find where 56 σ e γ/Σ p . The dependence of G versus the transverse size of the beam Σ p is mostly determined by the ratio H(kΣ p /γ)/Σ p and for a given value of k this function, and hence the gain factor, increases when Σ p becomes smaller. As a function of the chicane strength, the gain factor reaches maximum at q p = 1/κ p with the maximum value of G equal to where 'e'≈ 2.71 is the base of the natural logarithm. The plot of 2H(κ p )/eκ p is shown in For a drift length equal to a quarter of the plasma wavelength, L d = πc/2ω p , the sin function is equal to one. Note that in this case the gain factor (26) is negative (if q p > 0). It means that an amplification section also introduces a 180 degrees phase shift in harmonics of the plasma oscillations relative to their values at the beginning of the section.

V. COOLING RATE WITH ONE AMPLIFICATION SECTION
For the cooling time measured in the revolution periods in the ring, N c , without the amplification sections (that is in the absence of the drift L d and the chicane R (e,2) 56 in Fig. 1a), the following expression was derived in Ref. [3]: where Σ is the rms transverse beam size in the modulator and the kicker, σ h is the relative rms energy spread of the hadron beam, 56 σ h γ/Σ, and κ is the normalized wavenumber k, κ = kΣ/γ. The impedance Z is given by with 56 σ e γ/Σ and L m and L k the lengths of the modulator and the kicker, respectively.
With one amplification section, a density perturbation with a wavenumber k is amplified by the factor (26) and hence we need to multiply Z in Eq. (28) by G. Denoting the product Z p = ZG, we find, where Note the factor r in the normalization of the drift L d -since Ω scales as 1/Σ p this extra factor r makes l independent of the transverse size of the electron beam in the amplification section. Replacing Z in Eq. (28) by Z p we obtain where Note the negative sign in Eq. (32)-it means that in order to have cooling one of the chicanes (or all three of them) should have a negative value of R 56 (assuming that the sin function in the integral (33) is positive). This is due to the fact that an amplification section flips the phase of the density harmonics, as was indicated at the end of the previous section.
Because this expression is symmetric with respect to the interchange of the three variables |q 1 |, |q 2 | and |q h |, the maximum of I 1 is attained when they are all equal, |q h | = |q 1 | = |q 2 | = q.
We then have We numerically maximized I 1 with respect to q for given values of the parameters r and l. maximum value of I 1,max is approximately equal to 0.042 which gives for the cooling rate Comparing this result with Ref. [3], we conclude that an amplification section adds a factor of 0.75 1 σ e to the cooling rate.

VI. TWO AMPLIFICATION SECTIONS
The setup with two amplification sections is shown in Fig. 1b. The cooling rate for a two-section amplification is obtained by replacing the impedance Z in Eq. (29) by ZG 2 , which gives the following expression for the cooling rate: where (in this expression, as in Eq. (34), we have assumed that all dimensionless values of the chicane strength are equal, q h = q 1 = q 2 = q 3 = q). We numerically maximized I 2 with respect to q for several values of the parameters r and a range of the values of l, with the result shown in Fig. 7. Fig. 8 shows Taking as a reference value r = 0.2, we find from Fig. 7 that the maximum value of I 2,max ≈ 0.026. Substituting this value into Eq. (37) we obtain for the cooling rate Comparing this result with Ref. [3], we conclude that two amplification sections add a factor of 0.33A 2 , which is approximately equal to the squared amplification factor of one section (36).
Because the cooling rates (35) and (39) (which scale as I 3/2 e and I 2 e , respectively) depend on the local electron beam current that varies within the electron bunch, in application to practical problems, one has to average these equations taking into account the finite electron bunch length which we denote by σ and I 2 e that a hadron sees over many passages through the electron beam: Here, the average is meant as an integral over the longitudinal position z, using the hadron ) as a weighting function. For instance, we have For an electron beam several times shorter than the hadron one, we can neglect σ by I 2 e z in (39) (we recall that in these equations we have assumed the ratio Σ p /Σ = 0.2) we obtain for the cooling rate with one amplification section and with two amplification sections Note that in derivation of Eqs. (41) and (42) we ignored the fact that the normalized length l in Eqs. (33) and (38) also depends on the local beam current (because the plasma frequency Ω scales as √ I e ). This can be justified by the fact that if the argument of the sin function is chosen to be equal π/2 for the peak beam current (so that sin has a maximum value at this current), its variation near the maximum value is not so important as the direct dependence of the cooling time versus I e that is taken into account in Eqs. (41) and (42). Nevertheless, the derived equations should be considered as a rough approximation to the cooling rates.

VII. EFFECTIVE WAKE FIELD IN MBEC
In our preceding analysis we used the concept of effective impedance Z p (κ) which is obtained from the impedance without amplification (29) by multiplying Z by the gain factor G one (for one amplification section) or two (for two amplification sections) times. While Z p (κ) gives all that is needed to calculate the cooling rates, it is also instructive to analyze the wake field that is associated with this impedance. The wake field is defined by the following equation [3] w It has a meaning of the energy change (normalized by the square of the hadron charge, (Ze) 2 ) induced by one hadron on a second one located at distance z away due to the interaction through the electron beam (that passes through one or two amplification sections). For the case without an amplification, this wake field was calculated and plotted in Ref. [3] where it was shown that it is positive for z > 0 and negative for z < 0 with a negative derivative w (0) > 0 at the origin.
Plots of the wake field for one (w 1 ) and two (w 2 ) amplification sections for the optimized values r = 0.2, l = 0.85, and q e = 0.87 are shown in Fig. 9. The normalization factors for these wakes are different: w 01 = 2 3/2 AI e L m L k /πΣ 3 γ 2 I A σ e and w 02 = 4A 2 I e L m L k /πΣ 3 γ 2 I A σ e , where A is defined in Eqs. (31) and has a meaning of the amplification factor in one section (for the eRHIC parameters from Table I  this derivative negative) -this corresponds to the fact mentioned in Section IV that the amplification factor of one section is negative. Also note that both wakes make oscillations and change sign away from the origin (at z ≈ ±3.1Σ/γ) -a feature absent in the wake of Ref. [3]. The reason of such oscillatory behavior lies in the finite bandwidth of the gain factor G (see Fig. 4) which is localized in the region of small values of κ p . This sign change of the wake means that, for a given value of R (h) 56 , hadrons with a large energy deviation will be shifted longitudinally into the region where the cooling force changes sign and leads to further increase of the relative energy (the so called anti-cooling effect). Effects of these nature have been studied for classical stochastic cooling (see, e.g., [9]) and for the optical stochastic cooling [10]; they impose a constrain on the value of R  i ) and then through one or two amplification sections, as shown in Fig 1. Finally, in the kicker, the hadron energy is changed from η (h) , where now f z,i denotes the force acting on the hadron from ith electron. This procedure is repeated M times and the cooling rate is estimated as an average over M runs of the difference (η (h) + ∆η (h) ) 2 − σ 2 h . In the drift sections of the amplification cascades we use the following equations of motion for the electrons, where Σ p is the rms size of the beam in the drift. We scale the energy deviation,p = η Σ p /r e , the coordinatez = zγ/Σ p , and the distance,s = ct/l, withl = γ(Σ 3 p /r e ) 1/2 , so that the equations of motion become dimensionless, The plasma frequency for the wavelengths with κ 1 is given by Ω in Eq. (18) which means that the plasma period in variables is where ν = n 0 Σ p /γ.
A different normalization of the variables is used in the modulator and the kicker: the energy deviation η is normalized by the rms energy spread of the electron beam, q = η/σ e , and z is normalized by the transverse size Σ of the beam in the modulator, ζ = zγ/Σ. With this normalization the energy change of an electron on length L m due to an interaction with a hadron is We denote by A 1 the factor in front of the sum in this equation. Note the relations between the variablesz andp in the amplifier and ζ and q in the kicker and the modulator: We first simulated the amplification of initial perturbations of small amplitude in the electron beam, as discussed in Section IV. We used N e = 10 5 electron macroparticles and the length of the "electron bunch" ∆z = 20Σ p /γ in the simulations. An initial density perturbation with the dimensionless wavenumber κ p and a relative amplitude of the density and parameter A from Eq. (31) is related to the amplification factor of one cascade. Two more parameters are the dimensionless strengths of the chicanes, q e and q h , defined in Sections VI and VIII. Finally, there is a ratio r = Σ p /Σ of the transverse size of the electron beam in the amplification section and in the modulator and the kicker.
Calculating numerical values of ν, A, A 1 and A 2 for the eRHIC parameters from Table I and assuming the electron peak current of I e = 30 A, we find ν = 1.5 × 10 6 , A = 24.5, Simulations with these values are difficult due to a required large number of macroparticles and small values of the interaction strengths, so we used larger values for A 1 and A 2 and smaller values for ν and A: ν = 5 × 10 4 , A = 10, Because With the dimensionless simulation parameters given by Eq. (51) we also simulated the cooling process with one amplification section. In this simulation, we used N e = 10 5 electron macroparticle and the length of the "electron bunch" ∆z = 10Σ/γ in the simulations. The averaging was done over M = 5 × 10 4 runs. The plot of the simulated cooling times as a function of the dimensionless chicane strength q is shown in Fig. 11 by blue squares. The solid curve is calculated using Eq. (32) with q h = q 1 = p 2 = q. One can see that Eq. (32) is in good agreement with the simulations which we consider as a confirmation of the correctness on our analytical results. In these simulations we assumed the ratio Σ p /Σ = 1.
In another set of simulations we used the ratio Σ p /Σ = 0.2 and varied the length l of the amplification section. The result of these simulations is shown in Fig. 12. The theoretical curve is calculated with the same Eq. (32). Again, we find a reasonably good agreement between theory and the simulation.

IX. SATURATION AND NOISE EFFECTS
In this section, we discuss the effects of diffusion and nonlinear behavior (saturation) for an MBEC configuration that utilizes amplification stages. In Ref. [3], it was shown that diffusion effects due to the noise in the hadron beam can be quantified by means of a diffusion coefficient D h , which is given by Here, Z is the impedance without amplification cascades, given by (29). For the cooling to overcome the diffusion, we need to satisfy the following requirement: With one stage of amplification, we need to replace Z with Z p = ZG, where G is the amplification gain of Eq. (26). For the case of two amplification stages, we instead replace Z by ZG 2 . Here, we will focus on the latter case, for which the effects of diffusion and saturation are more important. After working out the algebra, we obtain The equation given above assumes that all chicane strengths, q, are equal to each other which is indeed the optimized configuration. To take the finite length of the electron beam into consideration, we follow an averaging procedure entirely analogous to the one used for the cooling. In particular, after noting that D h ∝ I h I 4 e , we obtain where we have neglected the dependence of l with respect to I e . Moreover, we have made use of the property Combining this with Eq. (42), we obtain the formula Using a similar procedure, we can derive the diffusion rate of the hadrons due to the intrinsic noise in the electron beam. For the energy perturbation of the hadrons due to the electrons in the kicker, we use a formula analogous to Eq. (48) from Ref. [3] i.e. with The electron density perturbation δn e is now due to the shot noise in the e-beam, so we have δn e (z)δn e (z ) = n 0e δ(z − z ) for the case of no amplification. A general definition of the diffusion coefficient was given in Ref. [3] as Adapted for the new coefficient, this expression yields To accommodate the amplification effect, we first substitute into Eq. (62), to obtain Next, we note that, when including amplification stages, the Fourier quantityG k ∝ for the case of two amplification stages. Substituting the expression for the gain, we obtain The averaged diffusion rate is yielding a new diffusion ratio (again for r = 0.2, l = 1.0 and q = 1.1) Finally, we would like to address the issue of possible nonlinear behavior in the amplification cascade, an effect which can be important if the gain is large enough. Recalling our earlier analysis, we observe that the linearization of the Vlasov equation is valid provided that |δn| n 0e , where δn = ∞ −∞ dk e ikz δn k /2π is the density modulation of the electron beam. Let us assume, for the moment, that no amplification stages are present and δn represents the density perturbation after the electron chicane. In Ref. [3] it was shown that is the density perturbation of the hadrons in the modulator and F(k) ≡ Zn 0e g 0 (k)ζ 0 (k). The other functions mentioned here are given by and These are basically the expressions of Eqs. (53), (55) and (57) from Ref. [3], with some minor notation changes. In the case of amplification stages, we have F(k) ≡ Zn 0e g 0 (k)ζ 0 (k)G S (k), where (again) S is the number of stages and G is the gain factor. In all these cases, we have F(−k) = F * (k) so the linearity condition can be re-written as where we have defined a saturation measure I sat and made use of the property δn (M ) k δn (M ) k = 2πn 0h δ(k + k ) regarding the initial noise in the hadron beam. Collecting all the necessary terms, the result for two stages becomes

X. ESTIMATES FOR THE ERHIC COLLIDER
As a numerical illustration of the general theory developed in the previous sections we will estimate the optimized cooling rate for the nominal parameters of the electron-hadron collider eRHIC [11]. The parameters of the proton beam in eRHIC and hypothetical parameters of the electron beam in the cooling system are given in Table I.
Substituting parameters from Table I  Modulator and kicker length, L m and L k 40 m RHIC ring of 13 µs, this corresponds to 2.7 hours and 9 minutes cooling time, respectively.
For the two-stage case, l opt ≈ 1.0 for r = 0.2 so the length of the amplification section is L = l I A /I e Σγ 3/2 ≈ 83 m.
Using the results of the previous section, we can also estimate the diffusion and saturation effects for the eRHIC parameters. For the diffusion caused by the noise in the proton beam, from Eqs. (57), we find the ratio r 1 ≈ 0.90, and for the diffusion due to the noise in the electron beam, Eq. (68) yields r 2 ≈ 8 × 10 −2 . While both conditions, r 1 , r 2 < 1 are satisfied, the margin for r 1 is not large. From Eq. (72) we also find I max sat ∼ 0.85 which means that in this regime the nonlinear effects are essential. Thus, for the 9 min cooling time, the hadron diffusion and saturation neglected in our study are considerable. The situation can be mitigated by choosing a smaller chicane strength and slower cooling rate. For instance, using q = 0.3, we obtain a cooling time of 50 minutes, with r 1 ≈ 0.11, r 2 ≈ 4 × 10 −2 and I max sat ∼ 0.15. For these set of parameters, the linear theory of this paper provides a good approximation to reality.
A relatively small value q = 0.3 also helps with the cooling of hadrons at the tail of the beam energy distribution. As was discussed in Section VII, for the particles that are shifted longitudinally more than 3.1Σ/γ, the effective wake function changes sign and their energy spread increases with time (the anti-cooling effect). With q = 0.3, the strength of the hadron chicane is R (h) 56 = 0.3Σ/γσ h , and such particles lie at the far tail of the energy distribution, η > 10σ h , where their effect can be neglected.

XI. SUMMARY
In this paper, we derived the cooling rate for the longitudinal, or momentum, cooling using a simple 1D model that treats particles as charged disks interacting through the Coulomb force. Extending analysis of Ref. [3], we studied the cooling with one and two amplification sections in the system. In contrast to Ref. [3], where the noise effects are small, adding one or two cascades of signal amplification through a quarter of plasma wavelength drifts and chicanes also amplifies the noise. We have analyzed the role of the diffusion caused by the amplified noise in the electron beam. We also derived formulas that allow estimation of nonlinear effects in the amplification. These effects limit the maximum amplification level that can be used in an MBEC cooling device.
In our analysis, we assumed a round cross section of the beams with a Gaussian radial density distribution. This assumption can be easily dropped and other transverse distributions (e.g., with unequal vertical and horizontal sizes) used for the particle interaction.
This will only change the specific form of the interaction potential (47), with the rest of the calculations of the cooling rate remaining the same.
There are several effects that are neglected in our model. Clearly, the transverse dynamics due to the beam focusing is ignored, as well as longitudinal displacement of particles due to this focusing. We also ignored plasma oscillations in the electron beam in the modulator and the kicker regions. This is justified if the length of the modulator and the kicker is smaller than a quarter of the plasma period in the electron beam.
Finally, we note that the 1D theory can also be extended to include the effects of the transverse cooling. This type of cooling is achieved through the introduction of the dispersion in the modulator and the kicker regions, as it was proposed for the optical stochastic cooling scheme [12,13]. A preliminary consideration of the horizontal emittance cooling in MBEC has been carried out in Ref. [14].