Mitigation of laser plasma parametric instabilities with broadband lasers

Laser plasma instabilities (LPIs) cause laser energy loss, asymmetric and insufficient compression, and target preheating, thus are assumed to be among the major concerns of inertial confinement fusion research. Mitigation of LPIs can enhance the laser–target coupling efficiency and optimize the target compression dynamics, which is critical for the realization of robust and high-efficiency fusion ignition. Broadband lasers with polychromatic components or random phases have been investigated for decades as an effective alternative to mitigate LPIs. Here, we present a brief overview on the progress of broadband LPIs, including the models of broadband lasers, the involved physics, the conditions for effective suppression of LPIs, and some schemes to produce broadband lasers.


Introduction
Large-scale plasmas are produced by laser irradiation in most of the inertial confinement fusion (ICF) schemes (Froula et al. 2010;Pesme et al. 2002;Hinkel et al. 2004;Tikhonchuk et al. 2019), where a target containing deuterium and tritium fuel is directly (direct drive) or indirectly (indirect drive) compressed by high-power laser beams (Craxton et al. 2015;Lindl et al. 2014;Hurricane et al. 2014;Betti and Hurricane 2016;He et al. 2016;Lan et al. 2017;Zhang et al. 2020a, b).During the interactions of high-power lasers with large-scale plasmas, laser plasma instabilities (LPIs) can be fully developed at the 10 picosecond time scale (Kruer 1988).The LPIs, such as stimulated Raman scattering (SRS), stimulated Brillouin scattering (SBS), and two-plasmon decay (TPD), are among the major obstacles to achieving fusion ignition, as they can cause significant laser energy loss (Ping et al. 2019;Rosenberg et al. 2018), asymmetric and insufficient compression (Igumenshchev et al. 2012;Strozzi et al. 2017;Moody et al. 2012), and target preheating (Dewald et al. 2010;Smalyuk et al. 2008;Christopherson et al. 2021;Montgomery 2016).Cross-beam energy transfer (CBET) may significantly affect the implosion velocity and ablation pressure in ICF typically involving many laser beams (Radha et al. 2016), where one of the pump beams is a seed light to drive SBS via beating with the other one.Even though the nonlinear processes of LPIs and their mitigation have been investigated for five decades (Kaw 2017), these problems have not yet been completely understood and solved, partially due to the parameter sensitivity and complicated coupling mechanism at large temporal and spatial scales.For the mitigation of LPIs, a few ideas have been proposed, such as the use of spike trains of uneven duration and delay (Albright et al. 2014), second-and third-harmonic mixing (Feng et al. 2020), application of strong DC magnetic fields (Liu et al. 2018;Zhou et al. 2021), polarization rotation (Barth and Fisch 2016;Zhong et al. 2021), as well as broadband lasers (Zhao et al. 2017b;Follett et al. 2019).
The recent achievement of the fusion energy yield of 1.3MJ at the National Ignition Facility (NIF) (Abu-Shawareb et al. 2022) indicates that the enhancement of the laser-target coupling efficiency remains a critical issue for achieving robust high-efficiency ignition.Compared to the compromise scheme of nearvacuum hohlraum (Hopkins et al. 2015), both the energy loss and the dynamic behavior of gas-filled hohlraum could be controlled using the broadband lasers, which may improve the system stability.The experiment result about shock ignition also indicates that a good control of both filamentation and SRS remains a key factor for the success of the shock ignition scheme (Cristoforetti et al. 2021).The basic physics associated with the mitigation of LPIs is to break the phasematching conditions so that the instabilities effectively grow within limited time and space.Currently, broadband lasers are considered to be promising for future ICF drivers.In this paper, we present an overview on relevant progress made on this topic.
Mitigation of LPIs via broadband lasers has been investigated from the early 1970s.The early broadband laser model was based upon the Kubo-Anderson process, where the laser phase jumps randomly within a range [0, 2 ] following a Poisson process (Laval et al. 1977;Lu 1988).The mean jump time Δt c (coher- ence time) of the random phase leads to a bandwidth ∼ 1∕Δt c .A simplification of the random-phase equations was usually made according to the pump frame transformation, where the group velocity of incident laser v 0 is assumed to be infinite while comparing to that of daughter waves near the quarter-critical density or other special configurations (Lu 1989;Divol 2007).The solutions of linear stochastic differential equations can be expressed by Greens function, and the growth rate is obtained from the imaginary part of the pole (Brissaud and Frisch 1974;Thomson 1975;Laval et al. 1977).The random-phase model was also used in the description of LPIs driven by smoothing lasers which introduce a small bandwidth though frequency modulation (Guzdar et al. 1991a;Divol 2007).Different from the random-phase model, low-coherence lasers are generated by a fundamental seed pulse with a large bandwidth (such as 50 nm), which have a broad spectrum at any time within the pulse laser duration (Gao et al. 2020).These kinds of lasers are often considered to be composed of many differentcolor components in numerical simulations (Estabrook and Kruer 1983;Zhao et al. 2017a).A sunlight-like model was proposed recently, which can describe a broadband laser with continuous frequency spectrum.Such lasers satisfy the Maxwells equations and can be easily implemented in simulation codes (Ma et al. 2021).Thomson and Karush (1974) investigated the mitigation of LPIs via a broadband laser with random phases in homogeneous plasmas, and found that the parametric instabilities can be reduced when the laser bandwidth is much larger than the linear growth rate.It was also found that the threshold of convective modes in inhomogeneous plasmas is enhanced by a large bandwidth (Thomson 1975).The reflectivity and nonlinear evolution of instabilities driven by broadband pumps were investigated via particle simulations, which show an enhancement of forward SRS via wave beating of different-color components in tenuous plasmas (Kruer et al. 1973;Estabrook et al. 1980).A considerable reduction of reflection was found when the bandwidth is ten times larger than the growth rate (Estabrook et al. 1981;Estabrook and Kruer 1983).The absolute and convective instability thresholds for incoherent pump were studied based on the envelop equations (Laval et al. 1976(Laval et al. , 1977)), and later the validity regions for the statistical description were established (Pesme et al. 2007).Bonnaud and Reisse (1986) found that the forward SRS is less sensitive to laser bandwidth than backward SRS, and the former is weakened only at large detunings.An effective Hamiltonian method was developed to describe the LPIs in the presence of space-time fluctuations in homogeneous and inhomogeneous plasmas, whose applications in the absolute instabilities near the quarter-critical density were presented in detail (Lu 1988(Lu , 1989)).Guzdar et al. (1991b) found that the interaction region of the daughter waves is extended by the bandwidth in inhomogeneous plasmas, which can compensate the reduction of temporal growth rate.DuBois et al. (1992) introduced the finite bandwidth description into the multibeam LPIs model.Control of SRS and hot-electron generation in the interactions of chirped pulse and low-density plasmas was found at the time scale less than one picosecond (Dodd and Umstadter 2001).However, the mitigation of SRS via chirped pulse fails at the 1 Page 4 of 36 nanosecond-scale LPIs (Zhong et al. 2021;Wu et al. 2001).The general dispersion relations of partially coherent light in a uniform plasma were obtained according to a generalized Wigner-Moyal statistical theory, which were then used to analyze the broadband effect on SRS and SBS (Santos et al. 2007;Brandão et al. 2021).Follett et al. (2018) concluded from three-dimensional (3D) simulations that 0.7% frequency detuning is sufficient to eliminate the hot electrons generated via TPD under the OMEGA laser conditions.They later evaluated the absolute instability threshold for SRS and TPD driven by a broadband laser numerically (Follett et al. 2019).Palastro et al. (2018) found that laser bandwidth has little effect on resonance absorption in the linear regime.Numerical simulations showed a significant reduction in CBET with bandwidth of 0.2-0.6%(Bates et al. 2018), which was confirmed in the experiments (Marozas et al. 2018).Control of CBET can significantly enhance the beam-target coupling efficiency (Döppner et al. 2020).Zhao et al. (2017bZhao et al. ( , 2019Zhao et al. ( , 2021) ) proposed the decoupling conditions of polychromatic lights in both homogeneous and inhomogeneous plasmas under single and multibeam configurations.Such polychromatic beams behave well both in the linear and nonlinear regimes of LPIs under weak decoupling conditions.
Some early experiments with low-intensity drive laser were reported to be consistent with the random-phase theoretical predictions (Obenschain et al. 1976).Mitigation of LPIs via temporal smoothing lasers or combining with polarization smoothing (PS) has been extensively studied in the experiments.The laser bandwidth introduced via smoothing by spectral dispersion (SSD) causes the speckle structure at the focal plane to vary in time due to the angular dispersion.For indirect-drive scheme, the experiments using 527 nm drive laser with a bandwidth of 0.1% indicate that the SRS levels in low-density region can be dramatically reduced (Montgomery et al. 1996).Recently, Seaton et al. (2022) found with particle simulation that the nonlinear ion behaviors can be suppressed by the introduced spectral dispersion.Beam smoothing suppresses the filamentation instability, which indirectly reduces the SRS backscatter.High-efficient laser coupling was found for a 527 nm drive beam smoothed with combined SSD and PS at laser intensity I 0 ≤ 7 × 10 14 W/cm 2 .Without SSD and PS the coupling is efficient only at lower intensities I 0 ≤ 3 × 10 14 W/cm 2 (Moody et al. 2009).The smoothing time scale of SSD is the decorrelation time of a speckle, whereas it is instantaneous for PS.Therefore, application of both 250 GHz of SSD and PS in combination produces a significantly reduction in the total backscattering from a long scale-length gasbag plasma (Moody et al. 2001).For laser and plasma scales relevant to NIF experiments, it was found that PS or SSD alone does not always effectively control LPIs.Theoretical work indicated that a small amount of laser bandwidth can strongly reduce the reflectivity from a spatially smoothed laser beam, and the effective bandwidth is doubled using PS (Divol 2007).Mitigation of LPIs via SSD was also found for direct drive.Shock ignition experiments performed on LMJ-PETAL showed the considerable decrease of SBS and SRS as well as hot-electron generation via SSD (Baton et al. 2020).Stimulated Raman side-scatter was found to be reduced by introducing one-dimensional (1D) SSD (Kang et al. 2021).
Based on the mitigation effect of SSD, a broadband laser driver scheme StarDriver containing a large number of beamlets was proposed (Eimerl et al. 2016).Each beamlet has a two-frequency two-dimensional (2D) SSD, with a total phase bandwidth up to 1.5 THz at 351 nm wavelength.As a result, the ensemble lasers have a bandwidth 2-10%.However, it was proposed that as many as 15,000 monochromatic beamlets are required to reach 1% asymptotic smoothness, which makes the system complicated.Some new smoothing schemes were proposed to reduce the amplitude modulation of SSD using multi-color spectral distribution (Fusaro et al. 2021) and to generate the light field with time-dependent polarization rotation (Zhong et al. 2021).Low-coherence green laser with energy over 400 J and bandwidth ∼ 0.6% 0 has been produced in KUNWU (Gao et al. 2020).The fourth generation laser for ultrabroadband experiments (FLUX) is under construction, which aims to generate 150 J and 1.5 ns ultraviolet incoherent pulses with fractional bandwidth over 1.5% (Dorrer and Spilatro 2022).
Significant efforts have been paid to produce the ultraviolet broadband lasers, and many new technologies have been proposed, such as stimulated rotational Raman scattering (Weaver et al. 2017), excimer lasers (Obenschain et al. 2015) and sum-frequency generation scheme (Dorrer et al. 2021;Dorrer and Spilatro 2022).However, no practical scheme has been demonstrated up to now for efficient third-harmonic generation (THG) of broadband beams (Dorrer et al. 2021).Similar to the control of LPIs, the phase-matching conditions required for the frequency conversion of broadband lasers are also hard to be satisfied, which lead to the inefficient conversion of near-infrared broadband pulses to ultraviolet.Zhao et al. (2022) proposed the polychromatic driver scheme with significant improvement of the conversion efficiency, where high-power polychromatic beamlets are generated via optical parametric amplification (OPA) (Cui et al. 2016;Eimerl et al. 1995;Banks et al. 2002;Haynam et al. 2007;Baisden et al. 2016;Mima 2010), and each beamlet is triple-frequency converted independently.The polychromatic scheme are based on the matured technologies, and both single-beam and collective LPIs are effectively mitigated, thus may pave the way towards realization of robust and high-efficiency fusion ignition.
In this paper, we present a review on the recent progress of LPI mitigation with broadband lasers.First, some models on broadband lasers are introduced in Sect. 2. The mitigation of LPIs by broadband lasers in homogeneous and inhomogeneous plasmas is discussed in Sects.3 and 4, respectively.In Sect.5, we introduce some broadband technologies for high-power laser systems, especially the polychromatic driver scheme.A summary and some discussions on the interactions of broadband lasers with plasma are given in Sect.6.

Description of broadband lasers
Five kinds of broadband lasers are mainly discussed here, i.e., chirped pulse, modulated pulse, polychromatic light, low-coherence laser and sunlight-like laser.We first investigate the chirped pulse with central frequency 0 = 2 ∕ , which can be described as Here, A(t), , and T are the pulse envelop, light period in vacuum, chirp rate and pulse duration, respectively.The pulse envelop can be flattened by increasing the parameter m, and the full-width at half-maximum (FWHM) of spectrum is √ ln 2(1 + 2 )∕ T for m = 1 .Note that the frequency of the chirped pulse var- ies with time, i.e., its spectrum is an instantaneously narrow band.For example, the pulse bandwidth is very small ∼ 0.05% 0 within the LPI developing time scale ∼ 10 ps regarding for a nanosecond pulse with a whole bandwidth even up to Δ 0 = 5% 0 .Therefore, chirped pulse may mitigate the convective SBS in a large- scale plasma, but has little effect on the long-time behavior of SRS.
Sinusoidal frequency modulation laser is considered here as an example of modulated pulse, which can be described by The bandwidth of a f is approximately given by Δ 0 ∕ 0 ≈ (2m + 1)f m , and f m is the frequency difference between two adjacent frequencies.The temporal part of the laser phase (t, y) = 0 t + m sin( 0 f m t + y) at the final focusing lens via SSD is the same to Eq. ( 2), where y is related to the phase shift across the lens due to grating dispersion (Moody et al. 2001).Note that Eq. ( 2) can be expanded to where a j = a 0 J j (m) , and J j (m) is the j-th order Bessel function of the first kind.For a finite beam number, the general expression of Eq. ( 3) is the polychromatic light model given by where j is a random phase between [− , ] , j and k j denote the frequency and wavenumber of j-th beam, respectively.The spectrum is truncated by the number of frequencies N for a finite bandwidth Δ 0 = N − 1 .The polychromatic model Eq. ( 4) has been widely used in the studies of broadband effects on LPI mitigation.Different power spectra of Eq. ( 4) can be applied such as Gaussian, Lorentzian, and flat profiles (Zhao et al. 2017b;Follett et al. 2019;Luo et al. 2022).For the flat spectrum, the amplitude of each color component is a j = a 0 ∕ √ N under the same incident energy of a monochromatic laser (Zhao et al. 2017b(Zhao et al. , 2019)).Low-coherence lasers are generated by a fundamental seed pulse with a large bandwidth, which have a broad spectrum at any time within the pulse duration (Gao et al. 2020).Therefore, it is always described by the polychromatic model with a continuous spectrum.
The above polychromatic light model can be extended to a sunlight-like laser with continuous spectra (Ma et al. 2021).If one introduces the amplitude frequency spectrum of a broadband laser light f ( ) and the phase frequency spectrum ( ) , (1) (2) where f ( ) is a predetermined positive function and ( ) is a random function of the frequency changing within [0, 2 ] , then the overall amplitude-phase frequency spectrum can be written as The electric field E(t) of a broadband laser in the time domain can then be obtained by the inverse Fourier transform of the amplitude-phase frequency spectrum , where F −1 denotes the inverse Fourier transform.The band- width of the sunlight-like laser is determined by the frequency range of f ( ) ≠ 0 .According to Eq. ( 5), the amplitude-phase frequency spectrum E( ) changes with different random-phase spectra.Therefore, when two linearly polarized broadband beams with different amplitude or phase spectra are synthesized in two orthogonal polarization directions respectively, one can obtain a new kind of laser with randomly changed polarization in time (Ma et al. 2021).This kind of light is named as sunlight-like laser, which may be generated by dividing a single low-coherence laser evenly into two pulses in time and then synthesizing them in two orthogonal polarization directions.Comparing to the typical PS method, the sunlight-like laser actually realizes temporal random polarization.
Coherence time is a good metric to describe the broadband laser with continuous spectra, which is defined as (Goodman 2015;Follett et al. 2019).Here, the angular brackets represent the time average.The bandwidth scale is usually estimated by the reciprocal of autocorrelation or coherent time (i.e., Δ 0 ∼ 1∕Δt c ) in the broadband LPI models with random-phase modulation (Thomson and Karush 1974;Lu 1988).

Random phase model of broadband lasers
Absolute instabilities are the major concern in homogeneous plasmas, which grows with time exponentially.The fundamental requirements for the development of LPIs are the wave-vector matching condition ⃗ k 0 = ⃗ k p + ⃗ k s and the frequency match- ing condition 0 = p + s , where ⃗ k 0 , ⃗ k p and ⃗ k s are the wave vectors of incident light, plasma wave and scattering light, and p and s are the frequencies of plasma wave and scattering light, respectively.In general, broadband laser mitigates LPIs by detuning the matching conditions.
Based on the three-wave coupling equations, Thomson and Karush (1974) studied the dynamic behavior of LPIs driven by a broadband laser , where m(t) is a stochastic component, and its autocorrelation time Δt c ∼ 1∕Δ 0 is much shorter than the characteristic growth time of the instability mode Δt g = 1∕Γ .The linear growth rate Γ is the imaginary part of ( 5) 1 Page 8 of 36 plasma-wave frequency.When Δ 0 ≫ Γ , the instability only effectively grows within time Δt c .Therefore, the growth rate is reduced by the factor Δt c ∕Δt g = Γ∕Δ 0 to be The reduction of growth rate can affect the threshold of LPIs.Under the dipolar case v 0 = ∞ , Laval et al. (1976) obtained the modified threshold based on the envelop equations with stochastic phase.When the damping rates of plasma wave p and scattering light s are approximately equal p ∼ s , and Δ 0 >  p |v s ∕v p | , the thresh- old of backscattering is given by Γ 2 ∕Δ 0 >  p .Otherwise the small bandwidth Δ 0 |v p ∕v s | ≪  p is just a perturbation of the usual coherent threshold.The corre- sponding threshold is estimated by where v s and v p are the group velocities of scattering light and plasma wave, respectively.
Analogous to the mitigation condition of temporal mode growth, Montgomery et al. (1996) gave the condition to control the filamentation by adding beam smoothing bandwidth, which is estimated by taking the coherence length of laser less than the characteristic growth scale of filamentation where H is a factor to account for thermal effects, T e is the electron temperature, and is the light wavelength.The relation between normalized amplitude a 0 and laser intensity I 0 is a 0 = √ I 0 (W∕cm 2 )[ (μm)] 2 ∕1.37 × 10 18 .For example, a bandwidth of Δ 0 ∕ 0 ∼ 4 × 10 −3 is needed for H = 2 at the intensity ∼ 10 15 W∕cm 2 .The fur- ther simulation results suggest more optimistic effect than that estimated by Eq. ( 8).Therefore, the theoretical estimate for filamentation mitigation is Δ 0 ∕ 0 ∼ 10 −3 , which supports the experimental observation using 527 nm drive laser (Montgomery et al. 1996;Moody et al. 2009).Santos et al. (2007) and Brandão et al. (2021) obtained the general dispersion relations for SRS and SBS driven by a low-coherence pump according to a generalized photon kinetics model where D( is the zero-order photon distribution function.For SRS, p = pe and D p = 2 − 2 pe , and for SBS, Here, Z is the ion charge, k pi and c s are the wavenumber and velocity of ion-acoustic wave, m e and m i are the electron ( 6) mass and ion mass, respectively.In general, k 2 c 2 ∕D p ≫ 1 .The one-dimensional water bag distribution was used for the zero-order photon distribution function.It was found that the conclusion that forward SRS is less sensitive to broadband fields than backward one, and the range of unstable wave number is increased by the bandwidth (Santos et al. 2007).Analogously, the backscattering growth rate of SBS is found to decrease with the increase of bandwidth in the weak coupling regime (Brandão et al. 2021), and at the large bandwidth limit, it is due to the widened unstable wave number.Actually, Eq. ( 10) can be reduced to Γ sat SBBS = Γ 2 SBBS ∕(4k 0 c∕ ) with Γ SBBS being the growth rate of backward SBS driven by a monochromatic laser, which corresponds to an ultra-broad bandwidth according to the work of Thomson and Karush (1974).Kruer et al. (1973) validated the mitigation of reflected light by random phase lasers with finite bandwidths using particle simulations.Zhao et al. (2017a) performed a series of particle-in-cell (PIC) simulations with different bandwidth models as shown in Figs. 1 and 2. The electron temperature and plasma density respectively are T e = 100 eV and n e = 0.08n c , i.e., plasma is initially in the fluid regime, and finally is heated to the kinetic regime.The amplitude of the monochromatic laser is a 0 = 0.04 .It is found that the broadband laser only reduces the linear growth of SRS, but has little effect on the saturation level and electron heating in the nonlinear stage even though the bandwidth Δ 0 > 10% 0 .The comparison of backscat- tering lights between sinusoidal-frequency-modulation laser and polychromatic light (10) exhibited in Fig. 1a indicates that the mitigation effects can be improved via dispersing the energy of the peak spectrum.Figure 1b shows that random polarizations have no significant improvement in the SRS mitigation, owning to the small changes of ponderomotive force though an average of large-scale space and time.The wavenumber-frequency distributions of Langmuir wave (k L , L ) displayed in Fig. 2 show a much wider wavenumber-frequency spectrum of the broadband light than that of the monochromatic laser, which indicates a broad range of the Langmuir-wave velocity for the broadband case.
The linear bandwidth theory is based on the fluid model, which cannot be directly used to describe the nonlinear process, since the latter involves the interaction between electrostatic fields and particles.Zhou et al. (2018) numerically studied the coupling of broadband lasers with the frequency shifted plasma wave in the deep nonlinear stage.According to the simulation results presented in Fig. 3, they found that much more hot electrons are produced by the broadband laser than the monochromatic laser at t = 2750 due to more intense frequency shift of Langmuir wave of the former.The inflationary spectrum coupling is also found in the broadband laser interaction with inhomogeneous plasmas, owning to the natural broadband of the Langmuir waves in inhomogeneous plasmas.The broad range of the plasma wavenumber may also lead to a significant level of wave turbulence (Zhao et al. 2017a;Brandão et al. 2021]).Actually, the amplitude modulation of low-coherence laser always leads to a peak intensity around twice higher of the average (Gao et al. 2020), which may eliminate the reduction effect of bandwidth by enhancing the temporal growth rate and developing nonlinear intermittent SRS (Liu et al. 2022).
The results shown above suggest that the broad bandwidth laser can still develop intense LPIs due to the coupling of different-color components, which indicate that the phase-matching conditions are satisfied in a finite spectrum range.The coupling of different-color components makes the LPI mitigation with broadband lasers ineffective at large scales.Therefore, the decoupling conditions of beamlets with different frequencies are critical to the effective applications of broadband lasers.

Decoupling conditions for the polychromatic light
For polychromatic lights, Estabrook et al. (1980Estabrook et al. ( , 1981)); Estabrook and Kruer (1983) found that each spectrum line acts independently if the frequency difference 0 = Δ 0 ∕(N − 1) of adjacent lines is much larger than the growth rate.The spa- tial gain was found to be proportional to Γ∕Δ 0 for Δ 0 ∕Γ > 10 .A large laser band- width Δ 0 ∕ 0 = 10-36% was used in their simulations, and filamentation and SRS are observed to be enhanced by the beating of discrete lines when the beat frequency is equal to pe .Bonnaud and Reisse (1986) found that forward SRS is insensitive to Δ 0 , and it is weakened only at a large frequency difference.However, the ambiguity condition  0 ≫ Γ fails to give an effective bandwidth model for the suppression of LPIs.Zhao et al. (2017b) first obtained the exact decoupling threshold based on the two-color laser model, which illustrates a physical picture of bandwidth coupling mechanism.As an example, we consider the development of SRS by two-color components (Zhao et al. 2021), which can be described by the following fluid equations where a 0 = a 1 + a 2 is the amplitude of incident light composed of two-frequency components 1 and 2 , v th is the electron thermal velocity, ⃗ Ãs = ⃗ Ãs1 + ⃗ Ãs2 and ñe are the perturbations of scattering light and electron density, respectively.In the strongly coupled regime, the scattering light developed by one of the incident lights can be shared by the other pump beam.Without loss of generality, the frequency coupling (11) Therefore, the Fourier analysis of Eqs. ( 11) and ( 12) in the strongly coupled regime leads to the dispersion relation where There- fore, the two beams are strongly coupled under a small frequency difference ) is the dominant term of D e2− for a large frequency difference  0 ≫ Γ , and correspondingly the two beams are decou- pled.Here, one comes to the same conclusion of Estabrook et al. (1981); Estabrook and Kruer (1983).
Based on the photon kinetics model (Santos et al. 2007), a similar dispersion relation can be obtained from Eq. ( 9) for two-photon beams  0 ( ⃗ where the two beams are strongly coupled.Therefore, Eq. ( 9) describes the dispersion relation driven by a finite-bandwidth laser with a continuous and slowly varying spectrum.
In the decoupling regime  0 ≫ Γ , different-color beams develop Langmuir waves without sharing the scattering light.Therefore, the off-resonant terms a 1 ⋅ ⃗ Ãs2 and a 2 ⋅ ⃗ Ãs1 are neglected, and the dispersion relation is reduced to For a broadband laser, the different-color beamlets propagate in the same direction.Therefore, Eqs. ( 13) and ( 14) can be solved in 1D regarding to the backward SRS.As exhibited in Fig. 4a, the two-color beams act as one monochromatic beam when their frequency difference is small enough 0 = 0.15% 0 .With the band- width increasing to 0 = 1% 0 , the instability regions are separated from each other, and each beam develops SRS independently.The numerical results are validated by the PIC simulations as illustrated in Fig. 4b.Note that the intensity of Langmuir wave is significantly reduced by the phase decoupling.Based on the phase (13) plots, the major parameter for the coupling of different-color beams is the width of instability region, which also indicates a finite valid range for the phase-matching conditions.One notes that the LPI intensity of long-wavelength pulse is stronger than that of short one due to the relatively higher density of the former.Equation ( 14) can also describe the collective SRS driven by decoupled multi-color beams with a certain intersecting angle (Zhao et al. 2021).In 2D geometry, both lasers with s-polarization (electric field components perpendicular to the incident plane) are mainly considered, which can develop more intense SRS than the p-polarized case (electric field components within the incident plane).According to the 2D phase plot displayed in Fig. 5a, the width of instability region for the j-th beam gradually decreases from the backscattering mode to the forward-scattering mode along the linear dispersion relation where = 3v 2 th ∕c 2 (1 − j ∕ pe ) is the thermal perturbation correction.Figure 5b indicates that the scattering light is decoupled from the other beam when the wavenumber difference of scattering lights k s = |k s1 − k s2 | is larger than the correspond- ing instability width.
One finds from Fig. 4c, d that Eq. ( 15) is validated at n e ≤ 0.23n c , since SRS backscatter is coupled with the forward scattering at n e ≈ 0.25n c .As an example, the instability width of backward SRS Δk LB can be obtained by substituting 15).The frequency detuning shifts the backward SRS wavenumber as . Therefore, the decoupling condition for SRS backscattering is  0 (dk L ∕d 0 ) > √ 2Δk LB , where their instability regions are no longer over- lapped in the phase space.In the underdense plasma with  pe ≪  0 , the decou- pling threshold of backward SRS can be reduced to where Γ SBRS = a 0 √ 0 pe ∕2 is the growth rate of backward SRS developed by a monochromatic laser.
Note that the frequency difference not only reduces the radius of the wavenumber circle r j = ( j − pe ) 0 ∕c √ 2 j − 2 j pe , but also shifts the center to the left ⃗ k j (1 −  j  0 ∕k 2 j c 2 ) as shown in Fig. 5a.Therefore, the mitigation condition of SRS forward scattering r j −  j  0 ∕k j c 2 > Δk LF is larger than that of SRS back- scattering Eq. ( 17), where Δk LF can be obtained by taking the wavenumber of ( 16) forward SRS mode into Eq.( 15).Even though forward SRS is insensitive to a small bandwidth (Bonnaud and Reisse 1986;Santos et al. 2007), frequency difference  0 ≳ 1% 0 is sufficient to mitigate forward SRS and decouple the scat- tering lights under a j < 0.01.
Considering for the filled gas in hohlraum, the developing conditions of SBS always satisfy a 0 ∼ 10 −2 and n e ≲ 0.1n c (Zhao et al. 2017b).In analogy to the SRS discussed above, the phase decoupling condition for backward SBS is where Γ SBBS = a 0 k 0 c pi ∕ √ 8 0 k 0 c s is the backward SBS growth rate developed by a monochromatic laser.In general, Γ SBRS ∕Γ SBBS > 2 under c s ∕c ∼ 10 −3 ; therefore, both SRS and SBS are effectively mitigated as long as Eq. ( 17) is satisfied.The speckles developed via beam smoothing, such as continuous phase plate (CPP), lead to the growth of transverse density perturbations that results in laser filamentation due to the nonuniform pressure.The growth of filamentation can be described by the same dispersion relation of SBS (Kruer 1988).Analogous to the two-dimensional dispersion relation of SRS discussed in the above, the decoupling condition of the laser bandwidth for filamentation is in the same order of that for backward SBS.
In the linear growth stage, the plasma instability level driven by the decoupled beamlets at different frequencies can be estimated as Δn j = N n exp(Γ j t) , where N is the beamlet number, n is the initial perturbation and Γ j is the growth rate of j-th beam.Analogously, the instability level developed by a single monochromatic beam is Δn 0 = n exp(Γ 0 t) .One always finds Δn j ∕Δn 0 > 1 at the time t < (1 + 1∕ √ N) ln(N)∕Γ 0 due to a j = a 0 ∕ √ N .For a large enough beamlet number N ≫ 1 , the instability level of decoupled beamlets is higher than that of the mono- chromatic beam if the latter is saturated before ln(N)∕Γ 0 but the former still grows.Therefore, there exists a maximum beamlet number N to mitigate LPIs in the fluid regime.However, the Langmuir wave is heavily damped under the ICF conditions.In hot plasmas without considering the damping of scattering light, the threshold for developing absolute instability is Γ j >  p √ |v s ∕v p |∕2 (Liu et al. 2019).Therefore, LPIs can be effectively mitigated by a polychromatic light when each beamlet intensity of the polychromatic light is below the developing threshold.According to the 1D simulation results shown in Fig. 6a, the polychromatic light can save at least 8.8% of the incident energy from LPIs in the hot plasmas compared to the monochromatic laser.A hot electron tail with 52.2 keV can be found for the monochromatic laser at 22.3 ps in Fig. 6b.However, the electrons have not been obviously heated at this time for the polychromatic light.

SRS mitigation by broadband sunlight-like lasers
The electric fields of two broadband laser ( E y and E z ) for constructing a sunlight- like beam are displayed in Fig. 7a, b, where both of them have the same amplitude  2021) performed a series of 1D PIC simulations.For a better comparison, three kinds of laser beams including monochromatic laser, conventional broadband laser, and sunlight-like laser were utilized in the simulations.Ions are immobile to evaluate the effect of SRS only in the simulations.The SRS reflectivity R shown in Fig. 8a indicates that the reflectivity steeply increases at a certain intensities and then slowly saturates (Montgomery et al. 2002;Yin et al. 2006).The intensity threshold of SRS is defined according to R = 1% .Under this crite- rion, one finds that the SRS intensity threshold is enhanced from 7.6 × 10 14 W/ cm 2 to 1.12 × 10 15 W/cm 2 for a conventional broadband laser with Δ 0 = 1% 0 .Note that the threshold is further enhanced to 1.42 × 10 15 W/cm 2 by the sunlight- like laser with the same bandwidth, which is about double of the monochromatic beam.The conventional broadband laser requires a nearly doubled bandwidth Δ 0 = 2% 0 to achieve the comparable mitigation effect of the sunlight-like laser.One finds from Fig. 8b that the sunlight-like laser as well as conventional broadband laser can delay the saturation time of the SRS, which is consistent with the studies of Zhao et al. (2015).
The instantaneous intensities of scattered lights developed by the conventional broadband laser and sunlight-like laser are found to be stronger than that of the monochromatic laser based on Fig. 8c.However, the average reflectivity R are 20.6%, 17.2% and 10.7% for the monochromatic, conventional broadband and sunlight-like laser, respectively.The electron energy distributions shown in Fig. 8d indicate that a significant reduction of hot electron generation can be found for the sunlight-like laser.Therefore, the sunlight-like laser may mitigate LPI further compared to the usual broadband laser due to the time-dependent random polarization, which will be experimentally investigated in the near future (Li et al. 2022).

Convective amplification of broadband lasers
Frequency bandwidth may partially compensate the wavenumber mismatch in inhomogeneous plasmas.Lu (1988) found that the inhomogeneous feature is significantly reduced at Δ 0 ∕K 0 ≫ 1 , where K 0 = k 0 − k p − k s is the wavenumber mis- match.Guzdar et al. (1991b) found that even though the temporal growth rate is reduced by the bandwidth the interaction time is extended where L is the density scale length.Note that the explicit value of the correlation function has been used to obtain Eq. ( 19) with an effective bandwidth Δ 0 (1 + |v s ∕v 0 |) , which can be reduced to Eq. ( 6) at |v s ∕v 0 | ≈ 1 .Finally, the overall amplification G ∼ Γ m Δt i becomes independent of bandwidth Δ 0 .As a special case of sinusoidal phase modulation with finite bandwidth, Wen et al. (2021) found that the resonance location has a nonzero velocity, and the scattered light is amplified continuously to achieve the maximum gain when its group velocity is equal to the velocity of the resonance location.Therefore, the maximum-gain condition is met when the indicator ≡ 2Δ 0 f m L∕ pe c satisfies = 1 .The simulation results dis- played in Fig. 9 demonstrate that the SRS gain peaks at = 1.16 , which is independ- ent of the density scale length.Different from the above random-phase model, Zhao et al. (2019) studied the polychromatic LPIs based on the Rosenbluth convective model (Rosenbluth 1972), where the linear wavenumber mismatch is introduced into the phase.Actually, the frequency difference 0 ∼ 10 −2 0 has little effect on the group velocity of scattering light.For example, v s ∕v s ∼ 10 −3 at n e = 0.1n c .When the frequency difference is much larger than the corresponding growth rate, LPIs driven by a two-color beam in a local region with phase coupling can be described by (21)  where the wavenumber mismatches are K , and the temporal effect of broadband lasers is included in the modified growth rate Γ m = Γ 2 ∕Δ 0 .The above Eqs.( 21) and ( 22) can be reduced to where K � 0 = dK 0 ∕dx and k � 0 = d k 0 ∕dx .Comparing to the monochromatic case, the drive terms of Eqs. ( 23) and ( 24) are reduced by a factor cos( k � 0 x 2 ∕4) .In the heavy damping plasmas, the saturation coefficient is obtained to the first order Equation ( 25) indicates that the two-color beams are coupled to develop convective instability in a common resonant region.However, the saturation level is decreased by the phase detuning as compared to the Rosenbluth gain saturation coefficient.Especially, for the SRS in an inhomogeneous plasma with density profile n e = n 0 (1 + x∕L) , Eq. ( 25) can be reduced to According to Eq. ( 26), the two-color beams are weakly coupled when As an example, for a plasma with n 0 = 0.08n c , T e = 3keV and L = 3000 , we have the threshold 0 ≈ 2.4% 0 .Therefore, for effective mitigation of SRS, the fre- quency difference between two adjacent components is in the order of 10 −2 0 .Note that the above linear model is suitable for describing the behavior of daughter waves in a finite density range.However, situations are different for the scattering light propagation in a large-scale plasma, because a scattering light produced by a certain color beam will be amplified again as a seed mode in the subsequent parametric excitation region where its frequency is equal to the scattering light developed by another low-frequency beam.An example is shown in Fig. 10, the scattering light developed by the laser with 1 = 1.02 0 at 0.1n c will be amplified again by the other color beam 2 = 0.98 0 as a seed mode at 0.078n c , since both of the scattering lights have an equal frequency s = 0.7 0 .The amplification coefficient is almost twice of Eq. ( 26), which counteracts the mitigation of broadband lasers.Therefore, the decoupling condition of SRS scattering lights in a finite inhomogeneous plasma with density range [n e (x 1 ), n e (x 2 )] (the corresponding Langmuir-wave frequency is , which is always larger than the thresh- old of Eq. ( 27) in a large-scale plasma.
The simulation result displayed in Fig. 11a shows that the polychromatic light with a finite bandwidth Δ 0 = 2% 0 can significantly reduce the strength of Lang- muir waves in inhomogeneous plasmas, and the saturation level is further reduced by increasing the bandwidth Δ 0 to 4% 0 .However, most of the beamlets are still coupled at the large beamlet number N = 20 according to Eq. ( 26).The reduc- tion of Langmuir wave leads to the suppression of hot-electron tail as shown in Fig. 11b.One notes that the frequency range of Langmuir wave is [0.33, 0.38] 0 in this simulation example.Considering the case with Δ 0 = 2% 0 and N = 2 , the spectra of the two scattering lights, respectively, are s1 = [0.61,0.66] 0 and s2 = [0.63,0.68] 0 , where an overlapping frequency range [0.63, 0.66] 0 can be found for them.The overlapping spectra are much higher than the sidebands since they are amplified twice as demonstrated in Fig. 11c.The sharing parts are shrunk with the increase of the bandwidth, and totally are separated at Δ 0 = 6% 0 .Based on the secondary amplification mechanism, the scattering lights are strongly coupled for N = 20 even at Δ 0 = 6% 0 as presented in Fig. 11d.
Different from SRS, the secondary amplification mechanism of scattering lights is the major concern for the mitigation of SBS and CBET (Zhao et al. 2019).The bandwidth of SBS scattering light can be estimated as Δ s ≲ 2Δ(k pi c s ) ∼ 0.004 0 with including the effect of plasma flow, where Δk pi ∼ 2 0 ∕c is the scale of ion- acoustic wavenumber, and c s ∼ 10 −3 c is the ion-acoustic velocity.In the linear regime, SBS and CBET can be effectively mitigated by the frequency detuning  0 ≳ 0.5% 0 between two-color beams.Numerical simulations performed by Sea- ton et al. (2022) indicate that laser bandwidth becomes highly effective in reducing the CBET gain exponent if it satisfies Δ 0 >  pi .Beams with sufficiently large bandwidth can induce both forwards (pump to seed) and reverse (seed to pump) 1 Page 22 of 36 energy transfers, and the net CBET is totally suppressed when the opposing transfers are balanced at a certain bandwidth.Beams with Gaussian and Lorentzian spectral profiles are close to optimal, producing a Δ −3 0 scaling of the CBET gain exponent, where CBET is insensitive to the damping rate of the ion-acoustic waves.However, the efficacy of bandwidth is not guaranteed in the nonlinear stage.(Bates et al. 2020) found that a bandwidth of about 8THz is required to eliminate CBET under the realistic laser features and the plasma conditions similar to those on the NIF.
In the multibeam configuration, considering for the case of two pumps with an incident angle v in vacuum, the common Langmuir wave shared by the monochromatic lasers at n e ∕n c ≤ cos 4 v ∕4 is divided into multiple weakly coupled modes by the polychromatic lights with a broad spectrum Δ 0 ≳ 3% 0 .The saturation coeffi- cient of Langmuir waves is reduced by the wave-vector detuning G m = G 0 [1 + cos( LB )] according to Eq. ( 26), where th n 0 ( 2 0 − 2 pe ) 3∕2 cos v for low-density plasmas (Zhao et al. 2021).Without loss of generality, assume that each polychromatic light has N beamlets, and every beamlet shares different common modes with all other polychromatic beamlets, i.e., there are N 2 backscattered common waves.For weakly coupled modes, the amplification coefficient for the polychromatic lights can be estimated as p = N 2 exp(1) , and that for the monochromatic lasers with the same energy is n = exp(N) based on a j = a 0 ∕ √ N .The amplification ratio p ∕ n is found to be less than one when N > 3 , and the energy fraction  2 p ∕ 2 n < 1 is always satisfied for any N > 1 .The scattering lights of different-color beamlets are decoupled from each other when Eq. ( 27) is satisfied, i.e.,  0 ≳ 1% 0 , and the secondary amplification is mitigated by the refraction of scattering light in 2D geometry.
One obtains the common wavevectors of forward and backward SRS as 16) at k Ly = 0 and v th = 0 .The wave-vector of beam 1 presented in Fig. 12a is O 1 (0.87, 0.17) 1 ∕c , and the common SRS backscatter can be found at | ⃗ k LB c| ≈ 1.17 1 , which is the most intensity mode at n e ∼ 0.188n c1 .The wave-vector of the scattering light developed by beam 1 is ⃗ , and the wave-vector of the other ≈ (1.17, −0.36) 1 are induced by the scattering light from the other beam.The intensities of the common Langmuir wave and the concomitant electrostatic field are found to be significantly reduced by the polychromatic lights with bandwidth Δ 0 = 3% 1 according to Fig. 12b.The polychromatic light is insufficient to drive intense plasma waves if the beamlets are weakly coupled, and therefore the nonlinear process is mitigated.The temporal growths of the electrostatic energies shown in Fig. 12c indicate that the mitigation of multibeam SRS by polychromatic lights can also be found in the nonlinear regime.The electron energy distributions in Fig. 12d also demonstrate that the polychromatic lights with a bandwidth Δ 0 = 3% 1 can effectively reduce the electron hot tail in the fully nonlinear regime.Follett et al. (2018) performed a series of 3D simulations with two different multicolor configurations, one is the polychromatic light and another is the tiled configuration (each beam is monochromatic, but the different beams have different frequencies).The incident beams were simulated with phase plates and polarization smoothing.The density scale length is L = 211 μm at 0.25 n c , and the electron and ion temperatures are T e = 2.6 keV and T i = 1 keV, respectively.They found that a polychromatic light with three colors and a 0.7% 0 (1.76 nm) detuning is sufficient to eliminate TPD-driven hot-electron generation at I 0 = 7 × 10 14 W/cm 2 .However, hot-electron generation cannot be obviously eliminated in the tiled configuration at the same laser intensity.The 2D simulations with normally incident plane waves 1 Page 24 of 36 indicated that the amplitude modulation of broadband lasers has a significant impact on the TPD threshold.When speckled beams are used, the threshold is always enhanced by increasing the beamlet number.Lu (1988) found an effective Hamiltonian H eff = H 0 + Δ 0 H 1 of the linear sto- chastic equation, where H 0 is the primary part describing the coherent effects, and Δ 0 H 1 is the finite bandwidth modification.The imaginary part of the eigenvalue of H eff gives the growth rate near 0.25n c in the pump frame, where v 0 = ∞ was assumed to eliminate the spatial variable x in the random phase term.A series of scaling results were found for SRS and TPD.As shown in Table 1, Follett et al. (2019) numerically evaluated the threshold of the absolute instabilities including SRS and TPD, and compared them to the theoretical predictions of Lu (1989).

Mitigation of absolute instabilities by broadband lasers
They finally found a qualitative match of the parameter scalings between them with significant quantitative differences.Lorentzian power spectrum of broadband lasers used in the simulations has an equal bandwidth to the assumption of the Kubo-Anderson process in the theoretical model, i.e., Δ 0 = 2∕Δt c .Approxi- mate scaling laws for the SRS and TPD thresholds with a modest bandwidth were obtained based on a combination of the simulation results and the analytical expressions At a large bandwidth, the scaling of ∕Δt c becomes linear.
For the polychromatic light, two absolute mechanisms are considered including the 90 • SRS side scattering and the degenerate modes around 0.25 n c .The disper- sion relation Eq. ( 14) provides the maximum decoupling condition for absolute side scattering in inhomogeneous plasmas due to the detuning of Langmuir waves.When the threshold Eq. ( 27) is satisfied, the convective couplings of sidescattered lights are detuned.Therefore, the effective mitigation of 90 • side-scattering lights can be found at  0 ≳ 1% 0 according to Eqs. ( 14) and ( 27).Under the multibeam con- figuration, the frequency difference of pump beams can shift the absolute region which introduces a wavenumber detuning into the SRS side scattering.
(28) a 2 SRS ≈ 0.067 Both of the absolute SRS and TPD are locally developed around 0.25n c , which is the critical factor for their decouplings.The absolute instability region of beam 1 can be obtained from the SRS and TPD dispersion relations in cold plasmas (Zhao et al. 2019(Zhao et al. , 2020)), which is the maximum valid range for the development of absolute instabilities in inhomogeneous plasmas.The polychromatic beamlets are detuned from each other at  1 ∕2 − a 1 k 1 c∕4 >  2 ∕2 + a 2 k 2 c∕4 , i.e., When the decoupling condition Eq. ( 31) is satisfied, different-color beamlets drive their respective absolute modes in non-overlapping density regions [( j ∕2 1 − a j k j c∕4 1 ) 2 n c1 , ( j ∕2 1 + a j j c∕4 1 ) 2 n c1 ] as shown in Fig. 13.For example, compared to the same energy of monochromatic laser with intensity I 0 ∼ 10 15 W∕cm 2 , absolute instabilities are developed independently by the polychromatic beamlets under beamlet number N = 5 and bandwidth Δ 0 = (N − 1) 0 = 1.7% 0 .According to the intensity threshold of absolute SRS in inhomogeneous plasmas, the whole level of absolute SRS is considerably mitigated by the decoupled beamlets when a j ≲ 0.323(L∕ j ) −2∕3 (Liu et al. 1974).For a large enough bandwidth, the above results indicate that the drive threshold is enhanced by the increase of beamlet number, which have been observed in the simulations performed by Follett et al. (2019).
As can be seen from Fig. 14a, b, the strength of TPD is reduced by the bandwidth Δ 0 = 2% 0 .However, TPD is still intense enough to produce abundant electrons as presented in Fig. 14d.When the decoupling condition is completely satisfied Δ 0 = 4% 0 > 2.7% 0 , TPD is almost totally suppressed as shown in Fig. 14c, d.
The temporal growth of the electrostatic energies exhibited in Fig. 15a validates that the saturation level and nonlinear evolution of absolute SRS can be reduced by a polychromatic light with large enough bandwidths.The electron energies are (31)  0 ≳ 0.433(a 1 + a 2 ).
diagnosed in the fully nonlinear regime t = 3000 , and the distributions are shown in Fig. 15b.The SRS intensity developed by two polychromatic beams with bandwidth Δ 0 = 3% 1 is comparable to that of one monochromatic laser.Single-beam SRS can be considerably suppressed by a polychromatic light with Δ 0 = 3% 1 , and the same mitigation effect is found at a larger bandwidth Δ 0 = 5% 1 for two polychromatic beams, where their distributions are almost overlapped without hotelectron generations.

Schemes for broadband laser drivers
There have been a few schemes for the construction of broadband laser drivers.Several proposed schemes are discussed briefly in this section.Different from the StarDriver scheme (Eimerl et al. 2016), which requires a large number of relatively narrowband lasers, the FLUX concept is aimed to produce broadband single beam.Sum-frequency generation of broadband optical pulses is use of a narrowband pulse at 526.5 nm, where the combination of angular dispersion and noncollinearity cancels out the wave-vector mismatch and its frequency deviation (Dorrer et al. 2021).The shaping approaches for the control of the spectral density and temporal shape of nanosecond spectrally incoherent pulses was proposed, which are applicable to high-energy facilities (Dorrer and Spilatro 2022).
Low-coherence green laser with energy over 400 J and bandwidth ∼0.6% 0 has been produced at the KUNWU (Gao et al. 2020), which may provide a platform to test the broadband theoretical model.Based on the broadband nature of second-harmonic generation (SHG) in Deuterated potassium dihydrogen phosphate (DKDP) crystals, low-coherence 529 nm lasers have been generated with conversion efficiency around 63% (Gao et al. 2020).The spectrum of the high-power broadband laser is a triangular profile with the FWHM ∼ 3.2 nm.An efficient smoothing effect, combining a CPP and the induced spatial incoherence (ISI), was found for the low-coherence pulse (Gao et al. 2020).The rms nonuniformity is reduced from 40.7 to 14.7%, and the smoothed focal spot has a shorter longitudinal speckle length.However, no practical scheme has been demonstrated up to now for efficient third-harmonic generation (THG) of broadband beams (Dorrer et al. 2021).None of the proposed technologies, such as stimulated rotational Raman scattering and excimer lasers, performs well in both of the sufficient bandwidth and high conversion efficiency (Weaver et al. 2017;Obenschain et al. 2015;Dorrer et al. 2021).Moreover, the commercial availability of diffraction gratings and nonlinear crystal lead to the extremely complicated and enormous cost of conventional broadband laser drivers.The amplitude modulation is a universal phenomenon of broadband lasers, which is one of the technical challenges for the optical damage threshold.The system simplification of broadband driver needs the innovations of optical technologies.The concept of polychromatic drivers is proposed as a feasible and effective way to mitigate the LPIs compared to the conventional broadband driver.The basic mechanism and preliminary design of decoupled polychromatic drive beams was first proposed in 2017 (Zhao et al. 2017b), and a scheme to realize such drivers is put forward recently (Zhao et al. 2022).
A typical polychromatic driver scheme was proposed based on the decoupling conditions of multi-color beams and the complexity of the project, where the full spectrum width is Δ 0 = 7.5% 0 with total 16 colors of monochromatic beamlets, and the differ- ence between every two adjacent frequencies is 0 = 0.5% 0 .The four-color beamlets with adjacent frequency difference 0 = 1% 0 are synthesized to form a single poly- chromatic beam with a spectrum width Δ 0 = 3% 0 .The frequency difference between any two beamlets of a polychromatic light is at least 1% 0 , which is sufficient to decouple the single-beam LPI submodes.The collective LPIs are also mitigated by the adjacent polychromatic beams due to the large enough beamlet number N = 4 and the minimum frequency detuning 0 = 0.5% 0 .With the same energy of a third-harmonic laser I 0 = 1.1 × 10 15 W/cm 2 , the beamlet intensity of a polychromatic light is significantly reduced to I j = 2.75 × 10 14 W/cm 2 .Therefore, the whole LPI intensity is effectively miti- gated when the beamlets are weakly coupled.
In general, a polychromatic driver is able to generate numerous polychromatic beamlets with whole spectrum width larger than 10% 0 (Cui et al. 2016).Therefore, the parameter range is large enough for design of the optimal scheme according to the drive configuration and facility cost.The two-color light generator, including an Nd:glass laser (1053 nm) and the other color signal light, is a unit of the polychromatic system as shown in Fig. 16.The second harmonic of Nd:glass laser with wavelength 526.5 nm is the pump source to amplify the signal light and simultaneously produce the idler beam via OPA.The generated twocolor beams are then second-harmonic and triple-frequency converted independently.The other polychromatic beamlets are produced in the same way, and parts of the beamlets are synthesized and focused via the dichroic mirror and off-axis parabolic mirror, respectively.
DKDP is one of the optimal nonlinear crystals currently used in the high-power laser systems due to its high optical quality for kilojoule pulses.The lengths of DKDP in the processes of OPA and SHG are determined by the intensity fraction of different-color beams.The total energy conversion efficiency of polychromatic ultraviolet light is about 60%, where additional 15% of energy is lost in the OPA process with comparing to the general monochromatic laser.The conversion efficiency of a fundamental laser to the third-harmonic laser with bandwidth 1.5% 0 as reported is only a third of the efficiency of the polychromatic light with a spectrum width 3% 0 .The suppression effects of such polychromatic lights have been investigated via large-scale PIC simulations, which indicate that more than 35% of the incident energy can be saved from the LPIs compared to monochromatic lasers for the direct-drive scheme, or high-density filled target for the indirect-drive scheme.The project cost and construction difficulty of the proposed scheme could be controlled, as moderate changes are required in the whole system.The conversion efficiency and beam collimation of the drive system could be enhanced by simplifying the beam synthesis, which will be investigated in the near future.

Discussion and summary
The interactions of broadband lasers with plasma at the nanosecond time and millimeter spatial scale may deserve further investigations, as such interactions may involve new physics associated with LPIs.For example, the self-adaptive coupling between the plasma evolution and broadband LPIs may have considerable impact on implosion dynamics and final yields.The introduced bandwidth can modify the properties of plasma in a long-time scale, including the energy deposition rate and hot electron generation, which in turn affect the subsequent LPI processes.
The production of hot electrons by broadband lasers has been comprehensively studied for n e ≤ 0.25n c under the same initial conditions of plasmas.The broadband lasers with coupled beamlets or continuous spectra can still produce a large number of hot electrons, especially near the quarter-critical density (Zhao et al. 2015(Zhao et al. , 2017a).The growth rates of forward and side-scatter SRS are found to be insensitive to the bandwidth (Santos et al. 2007;Zhao et al. 2022), which may produce a large number of hot electrons in the kinetic regime.The large frequency shift of Langmuir wave can compensate the three-wave mismatch introduced by the bandwidth (Zhou et al. 2018), which along with the broad range of plasma wavenumber could lead to a significant level of wave turbulence (Zhao et al. 2017a;Brandão et al. 2021).The amplitude modulation of broadband lasers has a significant impact on the LPI threshold, which may cause a series of nonlinearly collective bursts (Follett et al. 2018;Liu et al. 2022).However, the effect of amplitude modulation can be reduced by phase smoothing of incident lasers, and the TPD-driven hot-electron generation is greatly reduced by the smoothed laser with a small bandwidth 0.7% 0 (Follett et al. 2018).The decoupling threshold is critical to the effective suppression of hot electrons by broadband lasers.Actually, when the beamlets of a polychromatic light are decoupled, each beamlet is insufficient to drive intense plasma waves.Therefore, the nonlinear processes are suppressed effectively as shown in Figs. 6,11,12,13 and 14.The kinetic behaviors of broadband lasers near the critical density need to be explored further.Palastro et al. (2018) studied the linear and nonlinear resonance absorption of a broadband pulse with a fluid code.In the linear regime, the modified fractional absorption is f rm ≈ f r ( 0 ) 1 + (1 − 12q 3∕2 + 8q 3 ) Δ 0 ∕ 0 , where f r ( 0 ) is the fractional absorption of an electromagnetic wave with frequency 0 , and q = (2 L∕ ) 2∕3 sin 2 v .For a broad bandwidth Δ 0 = 10% 0 , the correction is only 2% near the peak absorption where q = 0.5 .Therefore, laser bandwidth fails to mitigate resonance absorption in the linear regime, where the ponderomotive response of the ions is neglected.In the nonlinear stage, the enhancement of absorption caused by the steep ion-density profile and electromagnetic decay instability can be suppressed by the broadband lasers.Broadband lasers can save considerable energies from SRS, TPD and SBS in the plasma n e < n c .The hot electrons generated by parametric decay instabilities around n c cannot be ignored in the case of the broadband lasers, since more energies could be deposited at the critical density than the monochromatic light.
In summary, mitigation of LPIs with broadband lasers can contribute to the enhancement of beam-target coupling and robust ignition performance.Decoupling conditions of polychromatic beams in homogeneous and inhomogeneous plasmas are identified and tested by numerical simulation.It is found that a laser bandwidth larger than 1% is typically required to mitigate LPIs effectively.This is applied to the design of broadband laser drivers.Up to now, no practical scheme has been demonstrated for efficient THG of broadband lasers.A new design of polychromatic drivers is proposed to effectively mitigate both the single and multibeam LPIs.It is based on the matured technologies and, thus, shows high feasibility and cost-effectiveness.

Fig. 1 a1
Fig. 1 a Temporal profile of the backscattered light waves.The polychromatic light and sinusoidal-frequency-modulation light are denoted as inc and sin, respectively.Both of them have the same bandwidth Δ 0 = 15% 0 .b Temporal evolutions of the electron temperature for the incident lights with different bandwidths or polarizations.Linear polarization and random polarization are denoted as lp and rp, respectively.The random polarization light is described by a rp = Σ N j=1 a j [cos(2 p j )e y + sin(2 p j )e z ] cos( j t) , where p j is a random number in [0, 1], and e y and e z are the unit vectors in the transverse directions [Reprinted from Zhao et al. (2017a), with the permission of AIP Publishing.]

Fig. 2
Fig. 2 Phase-space distributions of the Langmuir waves driven by a monochromatic light (a) and a polychromatic light (b), which are obtained for the time windows [301, 600] and [601, 1000] , respectively [Reprinted from Zhao et al. (2017a), with the permission of AIP Publishing.]

Fig. 3 a
Fig. 3 a Electron energy distribution at t = 2750 .b Spatial frequency spectrum of Langmuir wave at t = 2750 .The intensity of incident monochromatic laser is I 0 = 2.8 × 10 15 W/cm 2 .The plasma density is n e = 0.128n c with electron temperature T e = 3 keV [Reprinted from Zhou et al. (2018), with the permission of AIP Publishing.]

Fig. 4 a1
Fig. 4 a The growth rate Γ of backward SRS versus k L for n e = 0.08n c , and a 1 = a 2 = 0.02 under different 0 .b The corresponding phase plot of Langmuir wave at t = 340 from 1D PIC simulations, where the laser amplitudes and plasma density are the same to a. c, d show the comparison between analytical results of Eq. (15) and numerical solutions of SRS dispersion relation regarding the backscattering mode [Reprinted from Zhao et al. (2017b), with the permission of AIP Publishing.]

Fig. 5
Fig. 5 Distributions of the SRS growth rate Γ driven by two decoupled beams in the wavenumber space.The plasma density is n e = 0.2n c1 , and the electron temperature is T e = 2 keV, where n c1 is the critical density for beam 1.The amplitude of each beam is a 1 = a 2 = 0.01 with frequency difference 0 = 1 − 2 = 2% 1 .a The two beams propagate in the same direction, where the locations for stimulated forward (SFRS) and backward (SBRS) Raman scattering are shown.b The angle between the two beams is 16 • [From original figure in Ref. Zhao et al. (2022), © IOP Publishing.Reproduced with permission.All rights reserved.] Page 16 of 36 frequency spectrum but different phase frequency spectrum.The amplitude and phase frequency spectrum of E y are displayed in Fig. 7c.The sunlight-like laser, synthesized by E y and E z , has a random polarization with time, which is validated by the Stokes parameters shown in Fig. 7d.To some extent, a sunlight-like laser can be considered as a broadband elliptically polarized light whose polarization orientation, axial ratio and handedness changes randomly with time.For some relevant studies, it was found that the parametric backscattering can be reduced significantly when the polarization of incident laser light rotates slowly(Barth and Fisch 2016).Zhong et al. (2021) also found the time-dependent polarization rotation of incident laser can greatly suppress LPIs, since the ponderomotive force is reduced by the polarization difference between incident laser and scattering light in certain areas.To verify the advantages of sunlight-like lasers in mitigating LPIs,Ma et al.  (

Fig. 6 a
Fig. 6 a Transmission rate of incident energies at different plasma densities.b Electron energy distributions heated by monochromatic laser and polychromatic light at t = 15000 under n e = 0.07n c .Large-scale plasmas with a whole length of 1 mm are set in the simulations to include the convective effects.The polychromatic light contains four-color beamlets, and the electron temperature is 2.5 keV [From original figure in Ref. Zhao et al. (2022), © IOP Publishing.Reproduced with permission.All rights reserved.]

Fig. 7 1
Fig. 7 Temporal profiles of electric field components (a) E y and (b) E z .c The amplitude frequency spectrum f ( ) and phase frequency spectrum ( ) of E y .d Stokes parameters of the sunlight-like laser, where I denotes the intensity regardless of polarization, Q the linear polarization along the y ( + ) or z (−) axes, U the linear polarization at +45 • ( + ) or −45 • (−) from the y axis, V the right-handed ( + ) or lefthanded (−) circular polarization [Reprinted from Ma et al. (2021), with the permission of AIP Publishing.]

Fig. 8 a
Fig. 8 a The average SRS reflectivity within 6000 obtained from 1D PIC simulations as a function of laser intensity for monochromatic laser (MCL), broadband laser (BBL) with bandwidths Δ 0 = 1% 0 and 2% 0 , and sunlight-like laser (SLL) with bandwidth Δ 0 = 1% 0 , respectively.The purple dashed line denotes R = 1% .b The saturation time of the backward reflected light as a function of laser intensity for three kinds of laser beams.c The temporal evolution of the backscattering light.d The electron distributions for monochromatic laser, broadband laser and sunlight-like laser at 6000 .e, f show the temporal evolution of |E x | for broadband laser and sunlight-like laser, respectively.All the incident lasers have the same intensity scale ∼ 10 15 W/cm 2 , and the broadband laser and sunlight-like laser have the same bandwidth Δ 0 = 1% 0 .The plasma density is n e = 0.128n c and initial electron temperature is 3 keV [Reprinted from Ma et al. (2021), with the permission of AIP Publishing.]

Fig. 9 1
Fig. 9 The mean convective SRS gains, normalized by the Rosenbluth gain G r , as a function of the maximum-gain indicator for density scale length L = 200 μ m (squares), 330 μ m (circles), and 660 μ m (crosses).The plasma density linearly ranges from 0.09n c to 0.14n c with electron temperature T e = 0.1 keV.The intensity of monochromatic laser is I 0 = 5 × 10 14 W/cm 2 [Reprinted from Wen et al. (2021), with the permission of AIP Publishing.]

Fig. 10
Fig. 10 Schematic diagram as an example for the secondary amplification of backscattering light developed by two-color beamlets with frequencies 1 = 1.02 0 and 2 = 0.98 0 .The backscattering light developed at x 2 will be amplified again at x 1 [From original figure in Ref. Zhao et al. (2019), © IOP Publishing.Reproduced with permission.All rights reserved.]

Fig. 11 a
Fig. 11 a Temporal evolution of electrostatic energy driven by incident lights with different bandwidths.b Electron energy distributions for monochromatic laser and polychromatic light with different beam number N under the same energy at t = 5000 .Spectra of the backscattering lights for the polychromatic lights composed of c N = 2 and d N = 20 beamlets under the same energy.The density scaling length is L = 1000 , and the density range is [0.08, 0.12]n c .The initial electron temperature is T e = 2 keV.The incident monochromatic laser has a semi-infinite uniform amplitude a 0 = 0.014 (the corresponding intensity is I 0 = 2.5 × 10 15 W∕cm 2 with = 0.33 μm ) [From original figure in Ref. Zhao et al. (2019), © IOP Publishing.Reproduced with permission.All rights reserved.]

Fig. 12
Fig. 12 Wavenumber distributions of Langmuir wave driven by a two monochromatic lasers at = 1500 and b two polychromatic lights with bandwidth Δ 0 = 3% 1 at t = 3000 .c Temporal growth of electrostatic energies developed by pump beams with different bandwidths.d Electron energy distributions under different bandwidths at t = 6000 .The incident angles of beam 1 and beam 2 in vacuum are v1 = 10 • and v2 = −10 • , respectively.Plasma density of the beam overlapping region is n e ≈ 0.188n c , and the density scaling length is L = 1000 .The initial electron temperature is T e = 2 keV.The incident monochromatic beams have an equal amplitude a 1 = a 2 = 0.01 and same frequency 1 = 2 [From original figure in Ref. Zhao et al. (2021), © IOP Publishing.Reproduced with permission.All rights reserved.]

Fig. 14 1
Fig. 14 Spatial Fourier transform of the electric field at t = 800 under a Δ 0 = 0 , b Δ 0 = 2% 0 and c Δ 0 = 4% 0 , where N = 20 .The white line denotes the expected wavenumber distribution for plasmons occurring in the TPD instability.d Energy distributions of electrons for different bandwidth at t = 1200 .The density scaling length is L = 660 , and the initial electron temperature is T e = 2 keV.The amplitude of incident monochromatic laser is a 0 = 0.014 [From original figure in Ref. Zhao et al. (2019), © IOP Publishing.Reproduced with permission.All rights reserved.]

Fig. 15 a
Fig. 15 a Temporal growth of electrostatic energies developed by two pumps with different bandwidth.b Electron energy distributions under different bandwidth or beam number at t = 3000 .The incident angles of beam 1 and beam 2 in vacuum are v1 = 28 • and v2 = −28 • , respectively.Plasma density of the beam overlapping region is n e ≈ 0.243nc .The density scaling length is L = 1000 , and the initial electron temperature is T e = 3keV.The amplitude of incident monochromatic lasers have an equal amplitude a 1 = a 2 = 0.01 [From original figure in Ref. Zhao et al. (2021), © IOP Publishing.Reproduced with permission.All rights reserved.]

Fig. 16
Fig. 16 Optical design diagram of a two-color beamlet generator

Table 1
Follett et al. (2019)threshold power-law dependence on Δ 0 , L, , and T e from the analytical monochromatic result ( I thr,srs and I thr,tpd , respectively), the analytical result with bandwidth ( I Lu,srs and I Lu,tpd , respectively), and the fits to the result of laser-plasma simulation environment code ( I thr ) [Reprinted fromFollett et al. (2019), with the permission of AIPPublishing.]