A Wave-Based Simulation Model for Cross-Beam Energy Transfer and Stimulated Brillouin Scattering in Laser-Plasma Systems

We present WEBS (WavE-Based Simulations), an efficient wave-based simulation model designed to investigate the dynamic interplay between cross-beam energy transfer (CBET) and stimulated Brillouin scattering (SBS) in laser-plasma systems. By employing a unified Schrodinger-type envelope formulation for the laser and ion-acoustic waves, our model enables the use of a single, unconditionally stable Du Fort-Frankel numerical scheme, which maintains excellent long-term energy conservation even with coarse spatial grids. This approach not only achieves high computational efficiency validated against particle-in-cell simulations but also allows the selective activation or suppression of CBET and SBS processes, offering a clear diagnostic of their mutual coupling. Our simulations reveal that at high laser intensities, CBET and SBS reach a coupled steady state, leading to significant deviations from classical fluid theory predictions. Specifically, CBET gain is suppressed due to enhanced SBS reflectivity, while strong asymmetry in SBS reflectivity emerges between the interacting beams. These findings highlight regimes where the two instabilities strongly influence each other, providing critical insights for inertial confinement fusion research and offering a practical numerical tool for instability control and scenario design.

💡 Research Summary

The paper introduces WEBS (Wave‑Based Simulations), a novel wave‑based numerical framework designed to study the coupled dynamics of cross‑beam energy transfer (CBET) and stimulated Brillouin scattering (SBS) in laser‑plasma interactions. Traditional particle‑in‑cell (PIC) methods capture kinetic effects but are computationally expensive and noisy, while fluid‑based models are fast but lack detailed wave‑phase information and nonlinear coupling. WEBS bridges this gap by representing both electromagnetic laser fields and ion‑acoustic waves with a unified Schrödinger‑type envelope formulation. The three coupled first‑order equations for the pump envelope a₀, seed envelope a₁, and ion‑acoustic density perturbation ~nₑ are derived from the full Maxwell‑fluid system under the slowly‑varying envelope approximation, retaining key physics such as plasma flow (Vₓ, V_y), Landau damping ν, and nonlinear frequency shift δ ω.

A single unconditionally stable Du Fort‑Frankel scheme integrates these equations in time and space. The scheme uses three time levels (n‑1, n, n+½) and second‑order accuracy O(Δt²/Δx²), guaranteeing excellent long‑term energy conservation even on coarse grids. Because each update only requires nearest‑neighbor values, the algorithm is naturally parallelizable and suitable for large‑scale 2‑D or 3‑D simulations.

The authors validate WEBS against the PIC code EPOCH for a canonical 2‑D CBET scenario: two Gaussian beams intersect at 31°, each with peak intensity 3 × 10¹⁴ W cm⁻², wavelength 0.351 µm, and a uniform He plasma (n = 0.05 n_c, Tₑ = 1 keV, Tᵢ = 0.2 keV). Simulations with 10 cells per wavelength reproduce the PIC results for pump attenuation and seed amplification with high fidelity. Grid‑convergence tests (10, 15, 20 cells/λ) show negligible differences, demonstrating robustness against spatial resolution. Additional runs varying the plasma flow velocity V_y (0.5 Cₛ, 0.9 Cₛ, Cₛ) confirm the resonant CBET condition at V_y ≈ Cₛ, where power transfer peaks.

When the laser intensity is high enough for both CBET and SBS to be active, the simulations reveal a self‑consistent steady state: CBET initially transfers energy from the pump to the seed, which simultaneously drives ion‑acoustic waves. The amplified ion‑acoustic field enhances SBS reflectivity, causing a substantial fraction of the pump energy to be back‑scattered. Consequently, the net CBET gain is reduced by up to 30 % compared with classical fluid CBET theory, and a pronounced asymmetry appears in the SBS reflectivity of the two intersecting beams. The asymmetry is strongest when Landau damping is weak and the flow velocity matches the ion‑acoustic phase speed, indicating that the two instabilities cannot be treated independently in such regimes.

Key advantages of WEBS are: (1) a single scheme that conserves energy over long simulation times, (2) insensitivity to grid coarseness, allowing significant computational savings, (3) the ability to toggle CBET and SBS on or off, providing a clean diagnostic of their mutual coupling, and (4) quantitative agreement with both PIC benchmarks and established fluid theory, establishing confidence in the model.

The paper concludes with a roadmap for future development: incorporation of stimulated Raman scattering (SRS) and two‑plasmon decay (TPD) into the same framework, extension to fully three‑dimensional geometries with non‑uniform density and flow profiles, and large‑scale runs on high‑performance computing platforms to generate design‑space databases for inertial confinement fusion (ICF) experiments. WEBS thus offers a powerful, efficient tool for exploring and controlling laser‑plasma instabilities in realistic ICF scenarios.