Raman backscattering saturation due to coupling between ω p and 2 ω p modes in plasma

Raman backscattering ( RBS ) in plasma is the basis of plasma-based ampli ﬁ ers and is important in laser-driven fusion experiments. We show that saturation can arise from nonlinearities due to coupling between the fundamental and harmonic plasma wave modes for suf ﬁ ciently intense pump and seed pulses. We present a time-dependent analysis that shows that plasma wave phase shifts reach a maximum close to wavebreaking. The study contributes to a new understanding of RBS saturation for counter-propagating laser pulses.

Nonlinear Raman amplification has been demonstrated experimentally e.g. seed pulses have been amplified to 2 10 13 »´W cm −2 (from an initial intensity I 6 10 0 11 »´W cm −2 ) and 1. 3  amplification and pump depletion are accompanied by seed pulse compression, which leads to high intensities and an increased waveform steepness due to the high oscillation amplitude of the plasma wave. Nonlinear plasma-wave theory predicts that steepened plasma density perturbations include harmonics with amplitudes that increase relative to the fundamental plasma wave [27][28][29]. An early experiment provided evidence of plasma wave harmonic modes [30]. Moreover, a numerical study of nonlinear plasma waves has shown that second harmonic plasma waves strongly influence the time evolution of the fundamental through a resonant wavewave interaction [31].
In this paper, we investigate nonlinear growth and saturation of RBS in a parameter regime characterized by high density, cold, homogeneous and underdense plasma with intense initial seed and pump laser pulses. Parameters have been chosen to: (i) allow the seed to quickly enter the nonlinear pump depletion regime; (ii) ensure rapid saturation of nonlinear RBS to avoid development of Raman forward, modulational and filamenation instabilities; (iii) set a high wave breaking limit for the slow plasma wave (responsible for RBS, where v c c 2 , which has a significantly lower wave breaking limit compared with the fast plasma wave (v c p » ); and (iv) avoid Landau damping, which is higher in low density, warm plasma [32]. In this regime, beating between the pump and seed pulse drives large amplitude plasma waves that include harmonic modes [30,33]. We show, for the first time, how coupling between the p w and 2 p w plasma wave modes feed back into the p w mode, leading to a plasma wave phase shift. This phase shift then competes with pump depletion and wavebreaking as mechanisms of RBS saturation. It should be noted that the present work is different from the study of the fourwave instability in [34], where coupling between slow and fast plasma waves is considered for two counterpropagating laser pulses with 2 , 0 s p w w w w D = -= which does not allow 2 p w modes to be produced. We first present a theoretical model taking only the time-dependence into account, to identify and analyze the origin of the plasma wave phase shift. The structure of this paper is as follows: In section 2, a theoretical model to investigate coupling between fundamental and harmonic plasma wave modes is presented. Section 3 uses a 1D averaged particle-in-cell (PIC) code (aPIC) and a numerical analysis to study the nonlinearity arising from coupling between p w and 2 p w plasma wave modes. We also show the interplay between wavebreaking, pump depletion and the nonlinearity due to coupling between p w and 2 p w plasma wave modes at saturation.
Conclusions are presented in section 4.

Theoretical model
Consider the electric fields of the pump, seed and plasma waves propagating through the plasma, to be given by where j = 0, s, p represents the pump, seed and plasma waves, respectively, and k e , , j j ĵ w and E j are the propagation vector, frequency, polarization vector and the slowly varying amplitude of the three waves. The wave equation in plasma is is the total electric field vector, with the individual fields satisfying is the plasma current carried by the electrons moving The cold relativistic fluid equations are given by is the relativistic factor. Using the cold relativistic fluid equations, the rhs of equation (1) to third order in the fields is Equation (2) gives the linear and nonlinear current densities due to the three waves. Since the velocities of the plasma electrons under the influence of electric fields E j is v r t ee m , 2 for linearly polarized pump and seed, assuming that the plasma density perturbations obey Poisson's equation we can substitute equation (2) into equation (1), to obtain the time dependent [10,12,34] nonlinear mode-coupled equations for the slowly varying amplitudes E j as d in the last term on the rhs of equation (2), and account for mode coupling between 2 p w (arising from v n p · ▿d ) and p w plasma wave modes, while the remaining terms are due to relativistic effects. For electric field amplitudes of the where j e and j a are the real amplitudes and phases respectively, the set of equations (3) describes the nonlinear coupled time dependent evolution of amplitudes and phases of the three waves, including mode coupling and relativistic effects, but disregarding any propagation effects.

Numerical and simulation results
To demonstrate how the presence of plasma wave harmonics leads to mode coupling nonlinearities as predicted by the set of equations (3), we first present 1D simulations using the code aPIC [35], where an infinitely long pump beam interacts with a short duration counter-propagating Gaussian seed pulse with full width at half maximum (FWHM) duration of 20 fs. It should be noted that the initial short seed pulse will broaden as it gets Raman amplified [1,3,5], leading to seed FWHM durations significantly larger than the plasma wave period. This ensures the applicability of the three wave coupling model in the given parameter regime, as long as the comparison with the simulations is restricted to the leading edge region (LER) of the probe, since propagation effects are not included in the model. Figure 1(a) presents a snapshot from the simulation showing resonant interaction between the pump (dashed-red curve) and seed (solid-blue curve) pulses along with the generated plasma wave (solid-green curve). Figure 1(b) shows the Fourier transform of E p in the seed LER extending from front (right edge) of the simulation window to the location of the seed leading edge peak a s in figure 1(a). The plasma wave harmonic modes are clearly evident in the LER of the seed where there is resonant interaction between the three waves. Furthermore, in order to demonstrate the existence of plasma wave oscillations at 2 , p w we also present the temporal Fourier transform (for the interval between points (f) and (g) in figure 1(c)) of E p in figure 1(d). The figure clearly shows plasma wave oscillations at harmonic frequencies. Existence of the plasma wave harmonic modes has been demonstrated in early experiments by Umstadter et al [30].  The set of figures 2(a)-(g) show the evolution of the fundamental and harmonic plasma wave modes corresponding to figure 1. The figures show that during the initial stage (figure 2(a)) of interaction only the fundamental plasma wave mode is observed. However, as the plasma wave grows in amplitude with time (along with the seed in figure 1(c)), the harmonics of plasma wave appear (figure 2(b)) and increase in magnitude along with the fundamental till they become almost constant close to saturation (figure 1(c)) and then decrease. Also, it can be observed from the figures 1(b) and 2(g) that the restricted window for the Fourier transform of the seed LER region ( figure 1(b)) leads to lower magnitude fundamental and harmonic plasma wave modes compared with the case where the entire moving window is used ( figure 2(g)). Numerical solutions of equations (3), for the parameters used in the simulations, are presented in figure 3. The initial amplitudes and phases in figure 3 are, respectively, a mc e a mc e , , Wavebreaking is defined as 'the loss of periodicity in at least one of the macroscopically observable quantities', such as the plasma wave electric field amplitude E p . In terms of particle trapping, this occurs when a significant population of background plasma electrons are trapped [20]. Since the set of equation ( is the Lorentz factor relating to the phase velocity of the plasma wave [22]. Figure 1(a) shows that the peak plasma wave amplitude (green curve) of the seed LER is 88 » GV m −1 , which is consistent with the theoretical predictions of 80 and 90 GV m −1 in the presence and absence of the nonlinearities, respectively, at saturation in figure 3(a). As these values are significantly lower than the wavebreaking limit, for the two cases presented in figure 3, wavebreaking can be ruled out as a mechanism for saturation. This is clear in the phase space analysis of plasma electrons presented in figure 4, at the same time step as figure 1(a). The blue-dotted curves in the figure show that in the seed LER (figure 4(a)) there is no particle trapping. However, in the mid ( figure 4(b)) and tail ( figure 4(c)) sections of the phase space, background plasma electrons are trapped by the plasma wave, and their number increases from the seed pulse leading edge to the tail region. For the given plasma density and approximately an order of magnitude lower pump and seed intensities (non-relativistic) 1D PIC simulations presented in [36] show similar trapping for thermal plasma with a temperature of ≈50 eV. The solid-red curves show that, although a small fraction of background plasma electrons are trapped, the corresponding collective plasma electron oscillations (leading to E p ) preserve their periodicity in the three regions, thus clearly proving that wavebreaking does not occur. Figure 3(a) shows that in presence of nonlinearities due to coupling between p w and 2 p w plasma wave modes and relativistic effects (crosses), RBS saturates at early times compared with the case where nonlinearities are neglected (solid curves). Figure 3(b) shows that there is a significant shift in the plasma wave phase p a due to the mode coupling nonlinearities during saturation, whereas without them the wave phases do not change (solid-black curve). Figure 3(a) shows that switching off the mode coupling and relativistic nonlinearities leads to complete pump depletion (solid-green curve), which is reduced to 65% » (cross-green curve) when the nonlinearities are included. This leads to a reduction of the maximum seed and plasma wave amplitudes by 20% and 35% when nonlinearities are present, compared with the case where they are absent. This is due to the plasma wave phase shift, which detunes the resonant three wave system and reduces backscattering of the pump (lower pump depletion) into the seed, while in the absence of the nonlinearities the system is perfectly tuned and leads to complete pump depletion. Further, as will be shown later, among the two nonlinearities, relativistic effects are not significant for these cases. An important new observation in the present regime is that the saturation time is significantly smaller than observed in previous studies [10,12,24] that included a noise source a s and considerably higher a 0 . This can be explained, in the present regime, because the initial seed is sufficiently intense, even for low pump a 0 (compared with the previous studies) the ponderomotive beat between a 0 and a s is strong, which drives a large amplitude plasma wave with 2 p w mode that couples with p w mode causing the plasma wave phase shift, which saturates RBS at an earlier time. It should be noted that in the present case, for amplified seed a s and a 0 smaller than in the previous studies, relativistic effects are weak (as shown below). For the case of RBS seeded by noise, even though a 0 is large the ponderomotive beat is weak until the time when the seed is significantly amplified (which delays saturation). However, when a s has been amplified significantly, because of the higher a 0 (compared with the current regime) relativistic effects become important in causing the phase shift. Also, recent studies [13,14] have shown that for pump intensities much lower than that used in the present and previous [10,12,24] studies, the seed can be amplified to ultra-high intensity ( 10 18 » W cm −2 ), which is an order of magnitude higher leading to saturation due to relativistic effects. Figures 5(a) and (b), respectively, show the evolution of the plasma wave amplitude and phase (from equation (3)) for different plasma densities (all other parameters are the same as those used for figure 1) in the presence of mode coupling and relativistic nonlinearities. The figures show that for dense plasma ( 2.0 0 p w w = ), where the wavebreaking limit is high, and for the same initial pump and seed amplitudes, the plasma wave phase shift ( figure 5(b)) is reduced thus taking a longer time to saturate. This allows the plasma wave to grow to a large amplitude (figure 5(a) solid curve) before saturation. Lowering the plasma density ( 4 0 p w w = and 7) leads to an enhanced plasma wave phase shift, which causes saturation at an earlier time, leading to lower amplitude plasma waves. However, for the 7, 0 p w w = figure 5(a) shows that the plasma wave amplitude at saturation is 37 GV m −1 , which exceeds the wavebreaking limit of 34 GV m −1 . Wavebreaking is therefore the main mechanism of saturation, which occurs more quickly compared to saturation due to the plasma wave phase shift [9]. This can be verified from figure 6, where the dotted-blue and dashed-red curves in figure 6(a) show, respectively, the electron phase space for 7 0 p w w = and the location of seed pulse in the simulation box. It is observed from the dotted-blue curves in the figure, that a significant fraction of background plasma electrons are trapped in the seed LER region ( figure 6(b)), and significantly enhanced in the mid (figure 6(c)) region, while almost all plasma electrons are trapped in the tail ( figure 6(d)). The corresponding solid-red curves in the figures show that moving from the leading edge towards the tail region, increased plasma electron trapping leads to disruption of the periodicity in the collective plasma electron oscillations (observed as E p ), which eventually collapses and loses its structure in the tail region, which confirms the presence of wavebreaking. It should be noted that in such a strong particle trapping regime frequency shifts due to trapped particles [37] can also become important for non-relativistic seed and pump in warm plasmas and it can increase with the increase in temperature [36]. In figure 7 plasma wave phase evolution is plotted with (dashed-red curve) and without (solid curve) relativistic effects (using equation (3)), for a a 0.08 0 s = = (blue curves) and 0.04 (green curves). The figure shows that the mode coupling nonlinearity is the dominant physical process leading to a shift in the plasma wave phase in this parameter regime, where the pump and seed intensities are low and relativistic effects [10,12,24] (dashed curves) only mildly reduce the plasma wave phase shift (by 6.6% and 7.3% respectively for a a 0.08 0 s = = and 0.04).