2 Theoretical analysis

The wave number of the IAW excited by backward SBS can be calculated by

(1) $$\begin{eqnarray}k_{A}\unicode[STIX]{x1D706}_{De}\simeq 2\frac{v_{te}}{c}\sqrt{n_{c}/n_{e}-1},\end{eqnarray}$$

where $v_{te}=\sqrt{T_{e}/m_{e}}$ is the electron thermal velocity, and $n_{e},T_{e}$ and $m_{e}$ are the density, temperature and mass of the electron, respectively. Considering fully ionized, neutral, unmagnetized plasmas with the same temperature of all ion species ( $T_{H}=T_{C}=T_{i}$ ), the linear dispersion relation of the IAW in multi-ion species plasmas is given by^{[Reference Williams, Berger, Drake, Rubenchik, Bauer, Meyerhofer, Gaeris and Johnston20–Reference Feng, Xiao, Wang, Zheng, Liu, Cao and He22]}

(2) $$\begin{eqnarray}\unicode[STIX]{x1D716}(\unicode[STIX]{x1D714},k=k_{A})=1+\unicode[STIX]{x1D712}_{e}+\mathop{\sum }_{\unicode[STIX]{x1D6FD}}\unicode[STIX]{x1D712}_{i\unicode[STIX]{x1D6FD}}=0,\end{eqnarray}$$

where $\unicode[STIX]{x1D712}_{e}$ is the susceptibility of the electron and $\unicode[STIX]{x1D712}_{i\unicode[STIX]{x1D6FD}}$ is the susceptibility of the ion $\unicode[STIX]{x1D6FD}$ .

When strong pump light interacts with plasmas, the strong-coupling mode is generated, which grows with time and is not damped, and is called the quasimode. The strong-coupling regime of SBS is characterized by $(v_{os}/v_{te})^{2}>4k_{0}c_{s}\unicode[STIX]{x1D714}_{0}/\unicode[STIX]{x1D714}_{pe}^{2}$ , where $v_{te}=\sqrt{T_{e}/m_{e}},\unicode[STIX]{x1D714}_{pe}=\sqrt{4\unicode[STIX]{x1D70B}n_{e}e^{2}/m_{e}}$ and $v_{os}=eE_{0}/m_{e}\unicode[STIX]{x1D714}_{0}$ are the electron thermal velocity, electron plasma frequency and electron quiver velocity, respectively. In this paper, strong laser refers to that with an intensity satisfying the condition of the strong-coupling regime of SBS. The dispersion relation of the quasimode in strong pump light is given by^{[Reference Drake, Kaw, Lee, Schmid, Liu and Rosenbluth23, Reference Liu, He, Zheng and Wang24]}

(3) $$\begin{eqnarray}\unicode[STIX]{x1D716}(\unicode[STIX]{x1D714}_{A},k_{A})=\frac{k_{A}^{2}v_{os}^{2}}{4}\unicode[STIX]{x1D712}_{e}\left(1+\mathop{\sum }_{\unicode[STIX]{x1D6FD}}\unicode[STIX]{x1D712}_{i\unicode[STIX]{x1D6FD}}\right)\left(\frac{1}{D^{-}}+\frac{1}{D^{+}}\right),\end{eqnarray}$$

where $D^{-}\equiv D(\unicode[STIX]{x1D714}_{A}-\unicode[STIX]{x1D714}_{0},k_{A}-k_{0})$ , $D^{+}\equiv D(\unicode[STIX]{x1D714}_{A}+\unicode[STIX]{x1D714}_{0},k_{A}+k_{0})$ and $D(\unicode[STIX]{x1D714},k)=k^{2}c^{2}+\unicode[STIX]{x1D714}_{pe}^{2}-\unicode[STIX]{x1D714}^{2}$ .

Under the condition $T_{e}=5~\text{keV}$ and $n_{e}=0.3n_{c}$ , one can obtain the wave number of the IAW, $k_{A}\unicode[STIX]{x1D706}_{De}=0.3$ , from Equation (1). By solving Equations (2) and (3), the contours of solutions of the IAW without pump light and with strong pump light are shown in Figure 1 under the condition $T_{i}=0.2T_{e}$ and $k_{A}\unicode[STIX]{x1D706}_{De}=0.3$ . Although there are infinite solutions of the IAW, as shown in Figure 1, the mode with the least Landau damping ( $|\text{Im}(\unicode[STIX]{x1D714})|$ ) is preferentially excited in SBS. There exist two groups of modes, called the ‘fast mode’ and the ‘slow mode’, in multi-ion species plasmas, as shown in Figure 1(a). Here, the fast mode and the slow mode refer to the least damped mode belonging to each class of mode. In the condition $T_{e}=5~\text{keV},T_{i}=0.2T_{e}$ and $n_{e}=0.3n_{c}$ in $\text{C}_{2}\text{H}$ plasmas, the laser intensity in the strong-coupling regime is $I_{0}>I_{sc}=1.98\times 10^{15}~\text{W}/\text{cm}^{2}$ . When strong pump light with an intensity greater than $I_{sc}=1.98\times 10^{15}~\text{W}/\text{cm}^{2}$ interacts with plasmas, the quasimode is excited. As shown in Figure 1(b), the imaginary part of the quasimode frequency is positive, which illustrates that the quasimode is a growing mode.

Figure 1.

Contours of solutions to the dispersion relations of (a) the fast IAW mode and the slow IAW mode without pump light and (b) the quasimode with strong pump light $I_{0}=1\times 10^{16}~\text{W}/\text{cm}^{2}$ . The red line is $\text{Re}[\unicode[STIX]{x1D716}]=0$ and the blue line is $\text{Im}[\unicode[STIX]{x1D716}]=0$ . The conditions are $T_{e}=5~\text{keV}$ , $T_{i}=0.2T_{e}$ , $n_{e}=0.3n_{c}$ and $k_{A}\unicode[STIX]{x1D706}_{De}=0.3$ in a $\text{C}_{2}\text{H}$ plasma.

3 Numerical simulation

A one dimension in space and one dimension in velocity (1D1V) Vlasov–Maxwell code^{[Reference Liu, Zhu, Cao, Zheng, He and Wang25]} is used to research the quasimode excited by SBS in multi-ion species plasmas. We choose the high-temperature and high-density region as an example: the electron temperature and electron density are $T_{e}=5~\text{keV}$ and $n_{e}=0.3n_{c}$ , where $n_{c}$ is the critical density for the incident laser. The electron density is taken to be higher than $0.25n_{c}$ – thus stimulated Raman scattering^{[Reference Feng, Zheng, Liu, Cao, Wang, Xiao and He26–Reference Weber and Riconda28]} and two-plasmon decay instability^{[Reference Xiao, Liu, Wu, Zheng and He29, Reference Xiao, Liu, Zheng and He30]} are excluded. We take C, $\text{C}_{2}\text{H}$ , CH, and H plasmas as typical examples since they are common in ICF^{[Reference He, Li, Fan, Wang, Liu, Lan, Wu and Ye1, Reference Glenzer, Froula, Divol, Dorr, Berger, Dixit, Hammel, Haynam, Hittinger, Holder, Jones, Kalantar, Landen, Langdon, Langer, MacGowan, Mackinnon, Meezan, Moses, Niemann, Still, Suter, Wallace, Williams and Young3]}. The ion temperature is $T_{i}=0.2T_{e}$ . The linearly polarized laser intensity is $I_{0}=1\times 10^{16}~\text{W}/\text{cm}^{2}$ , with wavelength $\unicode[STIX]{x1D706}_{0}=0.351~\unicode[STIX]{x03BC}\text{m}$ , and the intensity of the seed light from the left boundary is $I_{s}=10^{-6}I_{0}=1\times 10^{10}~\text{W}/\text{cm}^{2}$ with frequency $\unicode[STIX]{x1D714}_{s}=0.994\unicode[STIX]{x1D714}_{0}$ . The spatial scale is [0, $L_{x}$ ], discretized with $N_{x}=5000$ spatial grid points and a spatial step $\text{d}x=0.2c/\unicode[STIX]{x1D714}_{0}$ . The spatial length is $L_{x}=1000c/\unicode[STIX]{x1D714}_{0}\simeq 160\unicode[STIX]{x1D706}_{0}$ with $2\times 5\%L_{x}$ vacuum layers and $2\,\times \,5\%L_{x}$ collision layers on the two sides of the plasma boundaries. The plasmas located at the center, with a density scale length $L=0.8L_{x}$ , are collisionless. The incident laser propagates along the $x$ axis from left to right with outgoing boundary conditions. Strong collision damping layers are added into the two sides of the plasma boundaries ( $2\times 5\%L_{x}$ ) to damp electrostatic waves, such as IAWs at the boundaries, and decrease the effect of the sheath field. The electron velocity scale $[-0.8c,0.8c]$ and the ion velocity scale $[-0.03c,0.03c]$ are discretized with $2N_{v}+1$ ( $N_{v}=512$ ) grid points. The total simulation time is $t_{\text{end}}=1\times 10^{5}\unicode[STIX]{x1D714}_{0}^{-1}$ , discretized with $N_{t}=5\times 10^{5}$ and a time step $\text{d}t=0.2\unicode[STIX]{x1D714}_{0}^{-1}$ .

Figure 2.

Frequency spectrum of $E_{y}$ with the time range $t\in [0,1\times 10^{5}]\unicode[STIX]{x1D714}_{0}^{-1}$ at $x_{0}=25c/\unicode[STIX]{x1D714}_{0}$ . The parameters are $n_{e}=0.3n_{c},T_{e}=5~\text{keV},T_{i}=0.2T_{e}$ and $I_{0}=1\times 10^{16}~\text{W}/\text{cm}^{2}$ in a $\text{C}_{2}\text{H}$ plasma, the same as in Figure 1(b).

Figure 2 shows the spectra of the SBS scattered light in the case of $\text{C}_{2}\text{H}$ plasmas. Calculated by the Vlasov simulation, the two peak frequencies of scattered light in $\text{C}_{2}\text{H}$ plasmas are $\unicode[STIX]{x1D714}_{s}=0.9972\unicode[STIX]{x1D714}_{0},0.9966\unicode[STIX]{x1D714}_{0}$ . Calculated from Equation (2) and Equation (3), the frequency of the fast mode is $\unicode[STIX]{x1D714}_{f}=4.98\times 10^{-3}\unicode[STIX]{x1D714}_{pe}=2.73\times 10^{-3}\unicode[STIX]{x1D714}_{0}$ , the frequency of the slow mode is $\unicode[STIX]{x1D714}_{s}=3.49\times 10^{-3}\unicode[STIX]{x1D714}_{pe}=1.9\times 10^{-3}\unicode[STIX]{x1D714}_{0}$ , and the frequency of the quasimode with pump light $I_{0}=1\times 10^{16}~\text{W}/\text{cm}^{2}$ is $\unicode[STIX]{x1D714}_{q}=6.98\times 10^{-3}\unicode[STIX]{x1D714}_{pe}=3.82\times 10^{-3}\unicode[STIX]{x1D714}_{0}$ . Therefore, the corresponding theoretical frequencies of SBS scattered light of the fast mode, the slow mode and the quasimode are $\unicode[STIX]{x1D714}_{s}=0.9973\unicode[STIX]{x1D714}_{0},0.9981\unicode[STIX]{x1D714}_{0},0.9962\unicode[STIX]{x1D714}_{0}$ . The theoretical and simulation results are shown in Table 1. Compared to the theoretical frequencies, we can see that the two peaks of the SBS scattered light in Figure 2 are indeed the SBS of the fast mode and the SBS of the quasimode.

Figure 3.

(a) Evolution of the SBS reflectivities of different modes with time, where SBS is the total SBS with the frequency range $\unicode[STIX]{x1D714}\in [0.9\unicode[STIX]{x1D714}_{0},0.999\unicode[STIX]{x1D714}_{0}]$ , SBS of the fast mode with range $\unicode[STIX]{x1D714}\in [0.9968\unicode[STIX]{x1D714}_{0},0.9977\unicode[STIX]{x1D714}_{0}]$ and SBS of the quasimode with range $\unicode[STIX]{x1D714}\in [0.9960\unicode[STIX]{x1D714}_{0},0.9968\unicode[STIX]{x1D714}_{0}]$ . (b) Reflectivity and transmissivity of the total SBS. The condition is the same as in Figure 2.

Table 1.

Frequencies of different modes and the corresponding scattered light. The conditions are $T_{e}=5~\text{keV}$ , $T_{i}=0.2T_{e}$ , $n_{e}=0.3n_{c}$ , $k_{A}\unicode[STIX]{x1D706}_{De}=0.3$ and $I_{0}=1\times 10^{16}~\text{W}/\text{cm}^{2}$ in $\text{C}_{2}\text{H}$ plasmas.

Figure 4.

Evolution of the SBS reflectivities in different species plasmas with time.

Figure 3 gives the evolution of the SBS reflectivities of the fast mode and the quasimode with time. As shown in Figure 3(a), the SBS of the fast mode develops quickly before ${\sim}1\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ and saturates during $t\in [1\times 10^{4},3\times 10^{4}]\unicode[STIX]{x1D714}_{0}^{-1}$ . In this stage, the entire SBS reflectivity (blue line in Figure 3(a)) is from the SBS of the fast mode. After $3\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ , SBS of the quasimode develops and competes with the SBS of the fast mode, which leads to the burst behavior of the total SBS reflectivity. The SBS of the fast mode exhibits an obvious decrease when the SBS of the quasimode develops to the peak amplitude during $[3\times 10^{4},6\times 10^{4}]\unicode[STIX]{x1D714}_{0}^{-1}$ . After $6\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ , there is a contradictory trend between the SBS of the fast mode and the SBS of the quasimode, which illustrates the competition between the corresponding SBS of the two modes. The recurrence period of the total SBS reflectivity is $T_{R}=7.501\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}-6.494\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}=1.007\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ , as calculated from the two peaks of the SBS reflectivity shown in Figure 3(a); thus the recurrence frequency of SBS reflectivity is $\unicode[STIX]{x1D714}_{R}=1/T_{R}=9.93\times 10^{-5}\unicode[STIX]{x1D714}_{0}^{-1}$ . In the same way, the period of the SBS of the fast mode is $T_{f}=2.031\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ and the period of the SBS of the quasimode is $T_{q}=2.126\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ . Thus, the sum of the two periods is $\unicode[STIX]{x1D714}_{f}+\unicode[STIX]{x1D714}_{q}=1/T_{f}+1/T_{q}=9.63\times 10^{-5}\unicode[STIX]{x1D714}_{0}^{-1}$ , which is close to the recurrence period of the SBS reflectivity $T_{R}$ . This illustrates that the burst behavior of the SBS reflectivity is a result of competition between the SBS of the fast mode and the SBS of the quasimode, and the recurrence period of the SBS reflectivity burst comes from the sum of the periods of the SBS of the fast mode and the SBS of the quasimode. Figure 3(b) demonstrates the total reflectivity and transmissivity of the SBS. Since the burst behavior of the total reflectivity of SBS occurs due to competition between SBS of the fast mode and SBS of the quasimode, the transmissivity demonstrates the same burst behavior, which is complementary to the SBS reflectivity after $3\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ . Therefore, competition between SBS of the fast mode and SBS of the quasimode is an important factor in the burst behavior of SBS and the nonlinear saturation of SBS.

Figure 4 demonstrates the SBS reflectivities in different species plasmas. In stage I, the SBS reflectivities develop and then saturate at a nearly fixed level. The saturation time of the SBS of the electrostatic IAW mode can be calculated from $t_{\text{sat}}=\text{ln}(I_{\text{sat}}/I_{s})/\unicode[STIX]{x1D6FE}_{0B}$ , where $\unicode[STIX]{x1D6FE}_{0B}$ is the SBS maximum temporal growth rate, which is defined in Equation (7). Taking $\text{C}_{2}\text{H}$ as an example, $\unicode[STIX]{x1D6FE}_{0B}=2.2\times 10^{-3}\unicode[STIX]{x1D714}_{0}$ from Equation (7); thus, the saturation time of the scattering from $I_{s}=1\times 10^{-6}I_{0}$ to $I_{\text{sat}}\simeq I_{0}$ is $t_{\text{sat}}=\text{ln}(I_{\text{sat}}/I_{s})/\unicode[STIX]{x1D6FE}_{0B}\simeq \text{ln}(1/10^{-6})/(2.2\times 10^{-3}\unicode[STIX]{x1D714}_{0})=6.28\times 10^{3}\unicode[STIX]{x1D714}_{0}^{-1}$ . From Figure 3(a), the saturation time of the SBS of the fast mode with a reflectivity of approximately 1 is $t_{\text{sat}}\simeq 6.72\times 10^{3}\unicode[STIX]{x1D714}_{0}^{-1}$ , which is close to the prediction of the linear theoretical model. In stage II, when different modes, such as the fast IAW mode and the quasimode, are excited in SBS in our simulation conditions, other modes such as the slow IAW mode and the ion bulk (IBk) mode^{[Reference Feng, Zheng, Liu, Cao, Xiao, Wang, Zhang and He31]} could be excited in SBS in some conditions, and the competition among SBS of different modes leads to the low-frequency burst behavior of SBS reflectivity. As the ratio of C to H increases, from Figure 4(a) to Figure 4(d), the Landau damping of the IAW decreases and the SBS growth rate and SBS reflectivity increase. The saturation time $t_{\text{sat}}$ becomes shorter. When the SBS of the electrostatic IAW saturates, the IAW amplitude no longer increases and the quasimode starts to develop, as shown in Figure 3(a), which leads to competition between the quasimode and the IAW and thereby the burst behavior of SBS reflectivity. The earlier the SBS of the electrostatic IAW saturates, the earlier the burst behavior of SBS reflectivity occurs. As a result, the initial time of stage II will be earlier.

The burst period of SBS reflectivity $T_{R}$ is related to the maximum electrostatic field $E_{x}^{m}$ and the growth rate of SBS $\unicode[STIX]{x1D6FE}_{0B}$ . The electrostatic field grows from the initial noise level $E_{x}^{0}$ to the maximum electrostatic field $E_{x}^{m}$ in about half the burst period $T_{R}/2$ , satisfying the relation

(4) $$\begin{eqnarray}|E_{x}^{m}|^{2}=|E_{x}^{0}|^{2}\times \text{exp}(\unicode[STIX]{x1D6FE}_{0B}\times T_{R}/2).\end{eqnarray}$$

Thus, the approximate relation of the burst period can be expressed as

(5) $$\begin{eqnarray}T_{R}=\frac{4}{\unicode[STIX]{x1D6FE}_{0B}}\ln \frac{E_{x}^{m}}{E_{x}^{0}},\end{eqnarray}$$

where the maximum electrostatic field $E_{x}^{m}$ is decided by the wave breaking of the IAW. Taking $\text{C}_{2}\text{H}$ plasma as an example, $eE_{x}^{m}/(m_{e}\unicode[STIX]{x1D714}_{0}c)\simeq 4\times 10^{-3}$ can be obtained from the simulation. The noise electrostatic field $e\unicode[STIX]{x1D719}/T_{e}=\unicode[STIX]{x1D6FF}n/n=(I_{s}/I_{0})^{1/2}=10^{-3}$ ; thus, $eE_{x}^{0}/(m_{e}\unicode[STIX]{x1D714}_{pe}v_{te})=e\unicode[STIX]{x1D719}/T_{e}\times (k\unicode[STIX]{x1D706}_{De})=3\times 10^{-4}$ and $eE_{x}^{0}/(m_{e}\unicode[STIX]{x1D714}_{0}c)=eE_{x}^{0}/(m_{e}\unicode[STIX]{x1D714}_{pe}v_{te})\times (\unicode[STIX]{x1D714}_{pe}/\unicode[STIX]{x1D714}_{0})\times (v_{te}/c)=1.625\times 10^{-5}$ . As a result, one can calculate $T_{R}\simeq 1.001\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ from Equation (5), which is close to the burst period $T_{R}=1.007\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ from the Vlasov simulation. The maximum electrostatic field $E_{x}^{m}$ decided by the wave breaking of the IAW is related to the phase velocity of the IAW. In our simulation, the dominant IAW is the fast mode. The phase velocity increases and the maximum electrostatic field $E_{x}^{m}$ decreases in the plasma sequence CH, $\text{C}_{2}\text{H}$ , C and H. Thus, from CH, $\text{C}_{2}\text{H}$ , C to H plasmas in sequence, the burst period $T_{R}$ decreases, which is consistent with the results shown in Figure 4.

However, for C plasmas, the Landau damping of the IAW is very low and the gain of the SBS is very high in the condition $I_{0}=1\times 10^{16}~\text{W}/\text{cm}^{2}$ , and the SBS cascade occurs in stage III^{[Reference Feng, Liu, Zheng, Xiao, Wang, Zhang, Cao and He32]}. The low-frequency burst in stage II is due to the competition among SBS of different modes, while the high-frequency burst in stage III is a result of SBS cascade. These results in Figure 3 give a clear explanation of the low-frequency burst behavior in SBS, which is a common phenomenon of SBS in high-intensity laser–plasma interactions.

Figure 5.

(a) Early linear stage of SBS in different species plasmas. (b) Relation between the SBS reflectivity and SBS gain in different species plasmas, where the gains in multi-ion species plasmas, such as CH and $\text{C}_{2}\text{H}$ plasmas, are calculated by the kinetic theory, and the gains in single-ion species plasmas, such as H and C plasmas, are calculated by the fluid theory. The SBS reflectivities by the Vlasov simulation take the values at $t=1.3\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ .

Figure 5(a) shows the linear growth and saturation process of the SBS in the early stage, in which the conditions are the same as in Figure 4. We can see that the SBS in C plasmas increases the most quickly among the four cases, while the growth rate of the SBS in H plasmas is the slowest and the saturation level is the lowest. As the ratio of H to C in the plasma increases, the growth rate and saturation level clearly decrease, because of the increase in the IAW Landau damping. According to Figure 5(a), the SBS reflectivity values at the time $t=1.3\times 10^{4}\unicode[STIX]{x1D714}_{0}^{-1}$ are chosen as the saturation values, as shown in Figure 5(b), because they are stable at this time. The linear gain of SBS by fluid theory is given by

(6) $$\begin{eqnarray}G=2\frac{\unicode[STIX]{x1D6FE}_{0B}^{2}}{\unicode[STIX]{x1D708}_{A}v_{gs}}L,\end{eqnarray}$$

where

(7) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{0B}=\frac{1}{4}\sqrt{\frac{n_{e}}{n_{c}}}\frac{v_{0}}{v_{te}}\sqrt{\unicode[STIX]{x1D714}_{0}\unicode[STIX]{x1D714}_{A}}\end{eqnarray}$$

is the maximum temporal growth rate of SBS^{[Reference Berger, Suter, Divol, London, Chapman, Froula, Meezan, Neumayer and Glenzer33, Reference Lindl, Amendt, Berger, Glendinning, Glenzer, Haan, Kauffman, Landen and Suter34]} and $v_{0}=eA_{0}/m_{e}c$ is the electron quiver velocity. $v_{gs}=c^{2}k_{s}/\unicode[STIX]{x1D714}_{s}$ is the group velocity of the SBS scattered light and $L$ is the plasma density scale length, while $\unicode[STIX]{x1D714}_{A}\equiv \text{Re}(\unicode[STIX]{x1D714}_{A})$ and $\unicode[STIX]{x1D708}_{A}\equiv |\text{Im}(\unicode[STIX]{x1D714}_{A})|$ are the frequency and Landau damping of the IAW. The collision damping of the IAW can be neglected since the electron temperature is as high as $T_{e}=5~\text{keV}$ in our simulation; thus only the Landau damping of the IAW is considered. For single-ion species plasmas, such as H or C plasmas, the gain of SBS can be calculated by means of the fluid theory. However, for multi-ion species homogeneous plasmas, the SBS gain can be calculated using the kinetic theory^{[Reference Lindl, Amendt, Berger, Glendinning, Glenzer, Haan, Kauffman, Landen and Suter34]}:

(8) $$\begin{eqnarray}G(\unicode[STIX]{x1D714}_{s})=\frac{1}{4}\frac{k_{A}^{2}v_{0}^{2}L}{v_{gs}\unicode[STIX]{x1D714}_{s}}\text{Im}\left[\frac{\unicode[STIX]{x1D712}_{e}\left(1+\mathop{\sum }_{i}\unicode[STIX]{x1D712}_{i}\right)}{\unicode[STIX]{x1D716}(k_{s}-k_{0},\unicode[STIX]{x1D714}_{s}-\unicode[STIX]{x1D714}_{0})}\right],\end{eqnarray}$$

which is more precise than fluid theory. Here, subscripts $0,s,A$ represent the pump light, the SBS scattered light and the IAW. Under the strong damping condition $\unicode[STIX]{x1D708}_{A}/\unicode[STIX]{x1D6FE}_{0B}\times \sqrt{v_{gs}/v_{gA}}\gg 1$ ^{[Reference Forslund, Kindel and Lindman35]}, one obtains the Tang model^{[Reference Tang36]}:

(9) $$\begin{eqnarray}R(1-R)=\unicode[STIX]{x1D700}\{\text{exp}[G(1-R)]-R\},\end{eqnarray}$$

where $R$ is the reflectivity of SBS at the left boundary and $\unicode[STIX]{x1D700}$ is the seed light at the right boundary. If $R\ll 1$ , the Tang model can be approximated by the seed amplification equation:

(10) $$\begin{eqnarray}R=\unicode[STIX]{x1D700}\times \text{exp}(G).\end{eqnarray}$$

Although the strong damping condition in C plasmas or $\text{C}_{2}\text{H}$ plasmas is not satisfied, the Tang model can also be applied in predicting the SBS reflectivity in the linear saturation stage. As shown in Figure 5(b), the SBS reflectivities in different species plasmas are close to the Tang model. Also, the seed amplification equation cannot be applied to predict the SBS reflectivities in the condition of strong pump light, since the SBS reflectivity $R$ is not much lower than 1. Figure 5 shows the linear process, including a linear growth and saturation process, in which a linear theory such as the Tang model is applicable. However, in a nonlinear process, such as stage II and stage III in Figure 4, burst behavior occurs with the nonlinear saturation mechanism of SBS that has been explained in this paper.

5 Summary

In conclusion, the quasimode is excited by SBS in high-intensity laser–plasma interactions, which competes with the IAW excited by SBS. The competition between SBS of the quasimode and SBS of the fast mode in $\text{C}_{2}\text{H}$ plasmas demonstrates that this competition mechanism is the cause of the low-frequency burst behavior of SBS. In different species plasmas, the same low-frequency burst behavior will occur, thus illustrating that the competition mechanism is common in high-intensity laser–plasma interactions no matter what the plasmas are. These results also give a good explanation of the intermediate low-frequency burst and saturation process in C plasmas.