Relations between Shot Noise, Gain Bandwidth, and Saturation of Instabilities

: There are numerous instabilities present in charged particle beams that undergo exponential growth and reach saturation. In various applications, such as free-electron lasers or micro-bunching light sources, achieving saturation is desirable. Conversely, there are applications where these instabilities are utilized as linear broad-band amplifiers for signals embedded in the charged beam. In the latter scenario, the saturation of an instability induces non-linear distortions in the imprinted signal, thereby limiting the useful range of such amplifiers. Accurate evaluation of these instabilities necessitates a complete and comprehensive modeling approach that includes shot noise within the beam. Unfortunately, such modeling is not always feasible or practical. In this paper, we introduce a methodology utilizing the frequency and bandwidth of the instability as key parameters. Through this, we derive an estimation for the range of linear instability growth. Our derivation is conducted in a model-independent manner, making it applicable to a broad spectrum of instabilities. To validate our approach, we employ established and thoroughly benchmarked simulations with a free electron laser (FEL) code as well as self-consistent

Meanwhile, for a number of theoretical and practical estimations of system capabilities, it is sufficient to know the central frequency, bandwidth, and achievable gain before the saturation of an instability.The central frequency and the bandwidth of the instability can be frequently estimated from the beam parameters.For example, the low and upper bands of the instability can be determined from basic beam and system parameters.For a broad-band space charge-driven longitudinal instability [21], the low band (f min ) is typically determined by the transverse beam size, σ ⊥ , scaled by the beam's relativistic factor (β • γ, β = 1 − γ −2 ): where γ is the Lorentz factor and c is the speed of light.Similarly, energy (σ γ ) and angular (σ θ ) spreads in the beam define the upper band of the instability, where L is the length and R 56 is the longitudinal dispersion parameter of the system.The situation is even simpler with narrow bandwidth instabilities such as FEL.The FEL central frequency (f 0 ) and the wavelength (λ 0 ) are determined by the resonant FEL condition: where λ w and a w = eB 0 λ w 2 √ 2πm e c are the FEL wiggler's period and dimensionless wiggler parameter, correspondingly, with B 0 denoting the peak magnetic field in the wiggler, m e the electron mass and e the electron charge.Similarly, the relative bandwidth of FEL instability is determined by a dimensionless parameter ρ (see, e.g., [22]): At the same time, similar simple estimations for the maximum pre-saturation gain are not available for a generic instability.This is one of the reasons for our attempt to find an estimation for attainable linear gain using a minimum set of parameters, relying on rather generic assumptions about the onset of saturation.
It is known that any e-beam instabilities (further on, we use an electron beam as an example of any charged beam) can be described by a set of self-consistent Maxwell and Vlasov equations [23] (see English translation in Ref. [24]), [25].In the classical limit, Maxwell's equations are linear while the Vlasov equation is not, hence, the non-linear terms in the latter are fully responsible for the saturation of the instability.On the other hand, it is established that the Vlasov equation can be linearized when the density modulation is significantly smaller than the initial beam's density [26][27][28].In other words, the Vlasov equation becomes non-linear when the density modulation becomes comparable with the initial beam's density (|δn| ∼ n 0 ).
Our hypothesis is that |δn| ∼ n 0 corresponds to the onset of saturation of the instability and the saturation of its gain and the appearance of signal distortions in the case of an amplifier.In Section 2, for compactness, we consider a one-directional instability while dedicating Appendix A for a detailed derivation in the 3-dimensional (3D) case.
We represent the response of the Maxwell-Vlasov linearized equations in the form of the Green function [29,30].Using properties of the Green function, we derive limitations to the Green function and spectral gain of the system by saturation arising from shot (Poisson) noise in the beam.
In Section 3, we compare our estimation with the results of numerical simulation of gain saturation for FEL and plasma-cascade instabilities.In Section 4, we discuss both the usefulness as well as the limitations of the proposed method.

Limitation for the Green Function
Let us assume that instability is growing from an initial infinitesimally small density perturbation: where n 0 is assumed to be uniform at the typical correlation time of the problem.In this case, the linear response to the initial perturbation can be described by a Green function, which naturally depends on the azimuth along the beam's trajectory s [31,32]: Mathematically, the Green function is a response to a Dirac δ-function perturbation in the beam's density.In the case of the lossless propagation of the beam, the number of particles is preserved:, applying an important constraint on the integral of the Green function: To estimate the maximum possible growth of the instability, we assume that the amplified signal is small when compared with the initial noise in the beam and can be neglected.Let us consider a portion of the beam contacting N particles, extending from −T to T, with a constant density: where I is the beam's current and e is the charge of particles.Let us consider the minimal possible level of noise in the beam, which is presented by the Poisson shot noise: with random times t i of the particles and a well-established statistical relation: where brackets ⟨. ..⟩ indicate averaging on the microscopic scale ∆t, which contains large number of particles: ∆N=∆t • n 0 ≫ 1.It also indicates averaging over an ensemble of possible implementations of random distributions [33].In other words, we are using a traditional method of three scales to connect microscopic and macroscopic phenomena [34,35]: where τ G is the duration of the Green function, which diminishes at the 2T boundaries, i.e., In other words, we are considering instabilities with bandwidths significantly broader than the macroscopic frequency content in the beam structure.
Since the beam noise in Equation ( 5) has zero value, one can introduce a modified Green function with zero average value (using Equation (3)): with the following property: (8)     Physics 2024, 6 924 which can be straightforwardly proved: It is worth noticing that relations (7) indicate that In other words, there is little difference between the two Green functions G, G we introduced in this paper.
Let us introduce a properly defining duration time for the Green function: Now one can calculate the root mean square (RMS) noise at the location s: Since the positions of particles are uncorrelated and the modified Green function has a zero average value, one can prove that: For i ̸ = j, the ensemble average is equal to zero because one can expand the modified Green function into the Fourier series with all terms averaging zero: g m e i(m(t−t i )+k(t−t j )) Ωt = 0.
Therefore, the statistically averaged density perturbation can be expressed as follows: It is our hypothesis that non-linear effects (saturation) will arise when RMS density perturbations become comparable with the average density of the beam: This gives us an estimation for the maximum possible value of the Green function (which has the dimensionality of the particle's density): Since the correlation time is connected to the bandwidths of the instability (this relation is exact for Gaussian pulses: one has a connection between the beam density, the bandwidth of the instability, and the maximum values of the Green function: This estimation can be evaluated further for frequency dependent gain: where g(ω) is dimensionless.In contrast with t-dependent the Green function, the g(ω) can be called a gain of the instability.Since the Green function does not extend beyond the {−T,T} interval, one can rewrite Equation (10) as where we define the central frequency of the instability ω c and its effective bandwidth ∆ω G , with an expected value of Hence, we are obtaining an estimation for the maximum attainable linear gain: or in practical units, There are two typical Green function shapes: one is nearly a "single wiggle," and the other is a nearly periodic oscillation with a smooth envelope (as shown in Figure 1).A broadband response with σ ω ∼ ω c is typical for microbunching [2,3] or plasma-cascade [20] instabilities, while high-frequency but narrower bandwidths with σ ω ≪ ω c are typical for FEL instability [4][5][6][7][8][9].
While the above estimations are valid for any generic Green functions, it is possible to develop them further for a narrow-band instability.In this case, the Green function can be presented by a smooth complex (dimensionless) gain envelope, g τ : where represents the local portion of the particles participating in the FEL process, known as a relative bunching factor in FEL theory [22].Equation ( 17) allows us to estimate maximum values for g τ : which, in addition to a numerical factor, has the product of two key dimensionless parameters.The first parameter is the number of particles in the wavelet: where ν is the beam's velocity and λ 0 is the central wavelength of the oscillations (radiation).
The second parameter is the relative bandwidth of the instability, σ ω /ω c : For a "single wiggle" type of Green function, |g τ | max represents the number of electrons collected at the central peak in response to the initial distortion by a single particle.In this case, the numerical factor √ π ∼ 1 in Equation ( 18) should be ignored.It can be used as an estimation for a broad range of instabilities, which in practical units states and for relativistic beams moving close to the speed of light as a function of the radiation wavelength: Hence, the maximum attainable linear gain is proportional to the square root of the beam current.The gain is also proportional to the square root of the radiation wavelength, i.e., is inversely proportional to the square root of the central frequency of the instability.In all cases, the gain is proportional to the square root of the relative bandwidth of the instability.In other words, a broad-band instability can have a higher linear gain when compared with a narrow-band instability with the same central frequency (wavelength).

FEL Simulations
For the first numerical test of our model-independent estimations, we selected a wellstudied and well-understood FEL instability.We considered a statistically representative set of initial conditions at the FEL's entrance using the 3D FEL code, Genesis [36], for a number of random-shot noise sets.We simulated the evolution of the bunch and the EM wave in the FEL both with and without a δ-function-like perturbation using identical input for the shot noise.The length of the FEL was intentionally chosen to be sufficiently long so that saturation is observed at its exit.
We performed Genesis simulation in time-resolved mode to simulate the FEL amplification and saturation process.Genesis has an excellent shot noise generator and allows adding a δ-function-like perturbation in the form of an additional bunching factor in a single wavelet.Specifically, in our simulations, we considered that the FEL response (i.e., its coherence length) is much shorter than the length of the electron bunch.Thus, we sliced the electron bunch into wavelets (with steps of wavelength, λ o ), with each slice represented by its local bunching factor.
Genesis is a perfect tool to simulate the evolution of the bunching factors in a 3D FEL.With two runs (one with and one without input signal) and using identical shot noise, the response to the perturbation at the exit of the FEL was calculated by subtracting the bunching factors caused by the shot noise from that in the presence of the input signal.This subtraction was conducted for each random shot-noise seed.
Specifically, to extract the response function, we generated two electron-bunch distributions for each shot-noise set.The first distribution has a random Poisson shot noise, generated by Genesis.This is a typical and well-tested setting for self-amplified spontaneous emission (SASE) FEL [37] simulations.The second distribution is identical to the first one with the exception of a single wavelet located at the middle of the bunch wherein we added a relatively small fixed bunching factor as an input signal.Technically, we use a quiet start to generate one set slice with an initial bunching factor at the level of interest, e.g., 10 −4 , and superimpose this slice on the middle of the first distribution.We therefore propagated these two bunches (one only with the shot noise and the other with the shot noise and the perturbation) through the FEL to record the evolution of the bunching factors in each wavelet as a function of its longitudinal position in the wiggler.Being the difference between two complex numbers, such an FEL response is a complex function and can be described by both the amplitude and the phase.Finally, we normalized the response to the amplitude of the perturbation.
One critical consideration is that the amplitude of the perturbation had to be chosen so that the response is well above the numerical noise of the code but below the level when non-linear effects from the input signal in the FEL's dynamics start to appear.We verified the latter by the following steps: first, we increased the amplitude of the initial perturbation while observing the evolution in the FEL with a quiet start.When non-linearity (early saturation) was observed, we marked the amplitude as the threshold level.Then the amplitude of the initial perturbation was reduced from the threshold level by three orders-of-magnitude for the rest of the simulations.
To observe the evolution, the FEL response was calculated at various locations along the FEL.It is noteworthy that before reaching saturation, the FEL response is identical to the Green function as defined above.In other words, we also have the means to calculate its correlation length.
The bunching factor, at location s, defined as without the initial perturbation (b 1 ) and with the perturbation (b 2 ), can be expressed as: Thus, the bunching caused by the perturbation is the difference between these two results, which can be written as or specifically: The arctan2 function provides the correct values of the phase within the range of [−π, π].However, it introduces an artificial "jump" from −π to π when the phase crosses the boundary.This jump could result in superficial discontinuities of the bunching phase (modulo 2π).We monitored this and applied the correction (modulo 2π) whenever such a jump occurred.The gain of the initial perturbation then is where b 0 = b s (0) is the initial bunching factor of the signal.
Simulating the FEL process in Genesis using the entire electron bunch would require large number of slices (e.g., a full bunch could contain 10 5 or more slices).Simulating the full-bunch FEL would be time-consuming and unnecessary.
In FEL, the information propagates (slips) for one wavelength per wiggler period.Hence, the number of slices needed for extracting the response of the initial perturbation (the FEL's response and Green function) is equal to the number of wiggler periods, which is on the order of 100.Naturally, we added some margin to eliminate the influence of the boundaries.More specifically, 1000 slices with 16,384 macro particles per slice were used in the simulation.Additionally, 32 random shot-noise sets in Genesis were implemented for statistical purposes.
We studied four different FEL amplifiers with wavelengths spanning from infrared (IR) to X-rays: (1) an IR FEL amplifier, (2) a visible FEL amplifier, (3) a vacuum ultraviolet (VUV) FEL amplifier, and (4) a hard X-ray FEL amplifier.Table 1 lists the detailed parameters used in the simulations.
We studied each of the four cases shown in Table 1.For conciseness, we describe one case-the visible FEL-amplifier-in detail.Initially, we verified that a single-slice perturbation with bunching at 10 −4 does not saturate in our FEL amplifier in the absence of shot noise (see Figure 2).Secondly, we confirmed that the addition of an initial perturbation does not affect the saturation of the SASE signal.Then, we performed simulations with thirty-two random initial shot-noise sets to extract the overall evolution of bunching in the FEL (in Figure 3).For each location along the FEL's length, we calculated the average values of the FEL response as well as the RMS values of its amplitude and phase variations.Since the Green function is unique in the linear regime, the fluctuations in it are minimal (as shown in Figure 4).Hence, we evaluated its saturation by observing the increase in the fluctuation's phase and amplitude.Naturally, they coincided with FEL power and the transition of bunching from exponential growth into saturation at about 20 m (Figure 3).Either the FEL response or the Green function cannot be solely described by the amplitude of its peak value; there are also some other important features at that location, such as its peak and the phase of the bunching modulation.Figure 5 shows details of the FEL response to the perturbation as a function of location (path length) in the FEL.The entire process (slippage, amplification, and saturation of the FEL response) can be visualized by the evolution of the response's profile in the FEL.24)), expressed in units of FEL wavelength (wavelets).Responses are calculated for different locations in the FEL as indicated.
Figure 5 illustrates the evolution of the bunching amplitude induced by the initial perturbation as a function of the slippage (in units of radiation wavelength) at several locations in the FEL wiggler.The initial single-wavelet perturbation has a bunching amplitude of 10 −4 and zero phase.It is placed at the middle of the bunch with natural and random shot-noise, whose contribution was subtracted.The only non-trivial part comes from the non-linear interaction between the shot noise and the coherent signal induced by the perturbation.
Starting from a blip in the middle, the information (and the bunching) expands forwards one wavelength per wiggler period.Naturally, the FEL electromagnetic field (radiation), propagating with the speed of light, carries the information forward.Meanwhile, the peak bunching remains at the origin of the perturbation (at an initial amplitude of about 10 −4 ).Since we are interested in the evolution of the bunching beyond the initial perturbation, we plot all the data except that at the origin.
At around 12 m, the peak of the response extends above the level of the initial perturbation-this is the onset of the exponential growth of the response, as well as the crest of the bunching envelope propagating with the group velocity of the FEL, i.e., between the velocity of the electrons and the speed of light.Later, when the signal approaches saturation, the response also saturates and becomes less even.
While informative, Figure 5 does not provide a straightforward answer about what the maximum gain is.To study this, we plotted both the amplitude (from Figure 3) and the phase (from Figure 4) of the bunching factor with the maximum amplitudes as functions of the wiggler's length.In other words, we calculated a peak value of the response (excluding the origin) and the phase for each location in the wiggler.
While detailing the evolution specifically, this choice has one defect-the maximum can jump from one location to another.Specifically, as shown in Figure 6, for a while (s < 12 m), the maximum of the response is located at the first wavelet, but at about s = 12 m, it jumps forward.This explains the change in the response's functional behavior as depicted in Figure 4, as well as the phase jump in Figure 4.With this information, we can further explore the observed features of FEL evolution shown in Figures 3 and 4.
In Figure 4, the evolution of the average bunching factor (green curve) is typical of SASE FELs, with an initial formation period, followed by exponential growth, and finally, saturation with |b SASE | ∼ 1.The evolution of the FEL response to perturbation is different.It has a period of initial coherence build-up, followed by exponential growth, and then saturation.The maximum amplification of the initial perturbation is reached just before the FEL's saturation (at approximately 20 m) and reaches 27.3.
Figure 4 shows that after the initial build-up, the phase of the evolution of the response function is smooth until about 20 m into the FEL, followed by jumpy behavior after that.The RMS spread of the phase rapidly increases in that region as well.Its origin is the non-linear response of the FEL that mixes the response to the perturbation with random shot noise.In other words, saturation removes the information about the input of the FEL amplifier making use of this amplified signal for feedbacks difficult, if not impossible, to implement.
Figure 6 shows both the amplitude and the phase of the FEL's response at a few selected characteristic locations in the wiggler.
Figure 6a shows the initial single-wavelet perturbation itself.It is worth noticing that the initial density perturbation travels at the speed of the electron beam, while the response moves with a slightly higher velocity, i.e., the group velocity of the FEL.When the electrons propagate in the FEL wiggler, three wave-packets develop from the initial perturbation: the density, the energy modulation, and the electromagnetic (EM) radiation.The radiation wave-packet slips in front of one FEL wavelength per FEL wiggler period.After the initial process of correlations build-up between the wave-packets, exponential growth starts.At some point, the response overtakes the initial perturbation (Figure 6b).
The maximum of the bunching amplitude is reached at the onset of saturation.Figure 6c shows the FEL response at this location, which has a fullwidth at half-maximum (FWHM) length of about 150 FEL wavelengths.Finally, as shown in Figure 6d, the response saturates, in which both the amplitude and phase are randomized by the non-linear interaction with the shot noise of the beam.Naturally, this is a situation where amplifying a desirable signal while preserving its information is no longer achievable.

Simulations of Plasma-Cascade Instability (PCI)
To verify that our estimation ( 16) is applicable to broad-band instabilities, we explored the plasma-cascade instability [20] using a different simulation tool-SPACE [38-40]a parallel, relativistic, 3D electromagnetic (EM), particle-in-cell (PIC) code capable of simulating interactions between relativistic particle beams.We simulated the PCI, a broadband microbunching instability in the system shown in Figure 7.The exponential growth of the longitudinal density modulation in this system is caused by strong periodic (or semi-periodic) transverse focusing and the beam's density modulation.The latter results in plasma frequency modulation and corresponding parametric resonance instability [20].Specifically, the PCI system for this study is composed of eight strong-focusing solenoids, comprising seven cells, with a gain per cell of approximately 4.36 [20].The solenoids focus the electron beam transversely, with a narrow waist located between them.The parameters of the beam and the solenoids are listed in Table 2. Figure 8 shows the RMS beam size evolution of the electron beam propagating through the PCI system.Figure 9 depicts the frequency dependence of the gain that the system would have for an infinitesimally small input signal in the absence of short noise.While this level of gain cannot be achieved in practice in this system, it represents the Fourier spectrum of the system's Green function ( 14) and defines the central frequency and the bandwidth of the instability.
Hence, with a central frequency f c = 20.8THz and a bandwidth ∆f G = 12.6 THz, according to Equation ( 16) one should expect that the maximum gain in this PCI system cannot exceed 1986 at a frequency of 20.8 THz and 993 at a frequency of 30 THz, as we numerically evaluated.According to this estimate, we expected to reach saturation in six-seven cells of the PCI system.
To explore saturation in the PCI system, we performed SPACE-code simulations for density modulation at a frequency of 30 THz.In these simulations, we used a representation number equal to one-where one computational particle represents one real electron-and random positions of the particles to obtain real Poisson noise in the beam sample.Similar to the FEL simulation, we performed two simulation runs to find the effect of shot noise on the seed signal.The first run was performed without an input signal, with the electron beam having random (Poisson) noise.In the second simulation run, we added a weak 30 THz density modulation to the shot noise distribution identical to the first run.The difference between density modulations in the two runs, recorded in various locations along the system, represents the evolution of the 30 THz seed system in the system.We used three different seeds of the shot noise in such simulations.Figure 10 shows the gain evolution for a 30 THz initial signal along the PCI system.The observed maximum gain of the seed (towards the end of seven-cell PCI system) ranges from 558 to 1070 for different seeds.One thing worth noticing is that different random seeds have little influence on the gain of the signal before approaching saturation.To detail how the density modulation evolves when saturation is close, we show the longitudinal density evolution resulting from the first sample of the shot noise (Figure 11, left) and that resulting from adding the 30 THz seed in (Figure 11, right).The sinusoidal seed shape of the signal is well maintained in the early stages, where the noise background is far from saturation.However, significant distortion of the input (seed) signal appears as early as the sixth cell of the PCI system.It becomes completely deformed after seven cells.Similar to what has been discussed in the FEL case, it is important to note that we used a seed amplitude much smaller than that of the shot noise to avoid any non-linearity caused by the seed signal itself (Figure 11a,b).As depicted in Figure 11k, the seed signal is amplified without significant distortions until about s = 12 m in the PCI system.However, distortion of the amplified seed signal reaches approximately its amplitude at s = 13.3 m when the peak-to-peak density modulation caused by shot noise reaches or exceeds approximately 30% of the average density (see Figure 11m).The seed signal is completely distorted at s = 15.3 m, where the peak-to-peak density modulation reaches a value comparable to the average density (see Figure 11o).
To confirm that our observations are consistent, we observed the same behavior for all three random seeds of the shot noise.The amplified seed signal indeed preserves its shape at s = 12.12 m, where shot noise causes relative peak-to-peak density modulation from 18% to 25%.At the same time, the seed signal is completely distorted at s = 15.3 m, where noise-induced peak-to-peak density modulation is comparable to or even exceeds the average beam density.To be concise, we only include plots for distributions at s = 12.12 and s = 15.3 m for sets two and three in Figures 12 and 13, respectively.To summarize, we studied the evolution of a 30 THz sinusoidal signal in the PCI system with different random seeds for the background noise.The maximum gain we observed ranges from 558 to 1070 as compared to our theoretical prediction of 993.The background noise has little effect on the gain evolution until saturation is approached.
While near saturation, significant distortions of the signal would occur when the relative peak-to-peak density modulation, caused by the beam noise, reaches approximately 30% of the average beam density.

Results and Discussion
We derived a set of equations for estimating the gain saturation of any 1D instability.It is quite remarkable that this simplified estimate compares favorably with the results of 3D simulations of FEL and PCIs.
Our numerical studies of the visible FEL amplifier as detailed above have shown that the maximum gain in bunching is limited to about twenty-seven (averaged over all thirty-two random seeds).We conducted a similar analysis for each of the other three FELs listed in Table 1 and found their maximum non-saturated gain as well as their correlation lengths.Table 3 compares these values with the estimates given by Equations ( 19)- (21).In short, our FEL simulations are in agreement with our approximate Equation ( 21) for the maximum unsaturated gain.We repeated these simulations with thirty-two completely different sets of initial random noise to eliminate the possibility of the agreements being a result of a fortunate choice of random shot-noise sets.The average of all simulated gains with different random seeds was still consistent with our theoretical predictions.
For the case of a broad-band plasma-cascade instability at 30 THz, we developed dedicated 3D particle-in-cell simulation (using code SPACE) and found that the maximum gain ranges from 558 to 1070 for various random seeds.This result is in agreement with what our prediction of gain at 993 using Equation (16).
For a simplified 1D estimate, derived without accounting for many important 3D effects, the agreements are remarkably good.
Therefore, our derived estimations for the attainable linear gain using charged beam instabilities have quite good agreement with actual values for instabilities with both narrow and broad bandwidths.This is a reason why they can be used as a tool in investigating systems where linearity is essential, for example, coherent electron cooling or seeded FELs.
Nevertheless, it is also worth noticing that for high-fidelity amplification, the gain of such a system should be at least a factor 2 below the maximum values indicated by our estimation.In this case, the density modulation amplified noise is approximately 25% of the average density and the signal quality is preserved.
Since the noise in the beam plays a key role in the saturation and distortions of the input signal, Equations ( 19)-( 21) can be extended to estimate the effect of increased noise in the beam: by adding the term for the noise power from the random (Poisson) distribution (P 0 ) to that of the noise power in the real beam (P b ).

Conclusions
In this paper, we have derived simple formulas to estimate the maximum attainable gain in charged beam instabilities in the presence of shot noise.In addition to the peak current of the electron beam, a necessary parameter for this estimation is the correlation length of the amplifier response, the central frequency, and the bandwidth of the instability.We compared the estimation provided by the formula with direct 3D FEL and PCI simulations for various wavelengths using the Genesis code in the time-resolved mode.We observed quite good agreement between the simulations and our theoretical estimate.

Figure 1 .
Figure 1.Schemtic picture of the time dependence of the Green functions for broad-band (σ ω ∼ ω c ′ , red line) and narrow-band (σ ω ≪ ω c ′ , blue line) instabilities.ω c is the central frequency of the instability.Horizontal and vertical scales chosen for illustrative purposes.

Figure 2 .
Figure 2. Evolution of FEL response to a δ-function-like perturbation without shot noise.FEL response grows exponentially after initial stage where slippage occurs to transfer bunching information from initial seed (located in middle of bunch) forward.

Figure 3 .
Figure3.Evolution of amplitudes of bunching factor in visible FEL amplifier (averaged over bunch for thirty-two random shot-noise seeds): evolution of total bunching factor's amplitude, originating from shot noise (in green) and evolution of FEL response to a δ-function-like perturbation (in red).Peak value of FEL response (18) is plotted as a function of its longitudinal position in FEL.

Figure 4 .
Figure 4. Phase of peak of FEL response (23)-average phase (red) and RMS phase variation (green)as a function of position in FEL wiggler.Thirty-two random sets of shot noise were used to generate data.

Figure 5 .
Figure 5. Amplitudes of FEL response as a function of s (see Equation (24)), expressed in units of FEL wavelength (wavelets).Responses are calculated for different locations in the FEL as indicated.

Figure 6 .
Figure 6.Bunching as a function of wavelet (slice) number at four locations in FEL wiggler: (a) near entrance (s = 0.6 m) (where initial δ-function-like perturbation is placed at middle of bunch), (b) at s = 11.4 m (where amplitude of response is about to overtake initial perturbation), (c) at onset of saturation, at s = 19.8m, and (d) at s = 22.0 m, well into saturation phase.We used thirty-two random shot noise seeds to derive average and RMS values.Red color indicates average value of the amplitude.Green color shows RMS spread of bunching amplitudes.

Figure 7 .
Figure 7. Schematic of plasma-cascade instability (PCI) system beamline, where electron beam (in red) traverses through eight strong solenoids (in blue).

Figure 8 .
Figure 8. RMS beam size evolution in PCI system.Typical ratio between beam sizes in focusing solenoids and in waists is approximately 10.This results in an approximately 10-fold modulation of plasma frequency[20].

Figure 9 .
Figure 9. Unsaturated gain curve for seven-cell PCI system.Peak gain is 15,190 at frequency f c = 20.8THz and bandwidth ∆f G = 12.6 THz.

Figure 10 .
Figure 10.Evolution of 30 THz signal density modulation amplitude in seven-cell PCI system for three random noise seeds.

Figure 11 .
Figure 11.Evolution of longitudinal density distribution with average electron density of n 0 = 1.04 × 10 12 m −1 and sample of shot noise (left), and result from initial 30 THz density modulation (right, in presence of noise) for different locations as indicated.Top: initial density distributions.Considering the oscillating nature of density modulation amplitude in PCI[20], we selected locations where seed signal has maximum amplitudes.Red dashed-dotted curves in signal evolution are sinusoidal fits which indicate shape of initial signal.Distortion is identified when signal shape strongly deviates from sinusoidal fit (when s > 12.12 m).

Figure 12 .
Figure 12.Longitudinal density distributions at s = 12.12 m for three random samples of shot noise of amplified shot noise (left) and from initial 30 THz density modulation (right).Signal shape is preserved, although shot noise's density modulation fluctuates due to different random seeds.

Figure 13 .
Figure 13.Longitudinal density distributions at s = 15.3 m for three random samples of shot noise. of amplified shot noise (left) and from initial 30 THz density modulation (right).Signal's shape fluctuates significantly when shot noise is approaching saturation, i.e., peak-to-peak density modulation is close to average beam density.It is also worth noticing that one can find correlation between locations of peaks in amplified shot noise and distortions of the signal.

Table 2 .
Beam and system parameters for PCI simulation.