The influences of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer (CBET) are examined using two- and three-dimensional particle-in-cell simulations in parameter regimes relevant to recent CBET experiments at the OMEGA laser facility. In mid-Z plasma irradiated by an intense pump beam and weaker probe beam, ion trapping, collisional de-trapping, and plasma flow induced by thermal effects are shown to affect the CBET gain. Ion trapping can enhance or detune the CBET resonance [Nguyen et al., Phys. Plasmas 28, 082705 (2021)]. Collisional de-trapping can affect the CBET gain at low seed beam intensity near the onset threshold for ion trapping. Thermal-effects-induced flow can also detune the CBET resonance at a level comparable to that from trapping at low seed beam intensity. As a consequence, the CBET gain is sensitive to collisions and dimensionality at low seed beam intensity where ion trapping is weak but is insensitive to collisions and dimensionality at high seed beam intensity where ion trapping is strong.
Laser-plasma instabilities (LPI) reduce the energy coupling to plasma and can affect both the direct- and indirect-drive approaches to laser-driven initial confinement fusion (ICF).1–6 Two major LPI processes are stimulated Raman scattering (SRS) and stimulated Brillouin scattering (SBS), where the laser light scatters off an electron-plasma wave (EPW) and an ion-acoustic wave (IAW), respectively.7 The cross-beam energy transfer (CBET), a special case of SBS, allows crossing laser beams to exchange energy through the excitation of IAW. An understanding of LPI dynamics and saturation in a variety of parameter regimes would enable the coupling of nonlinear LPI effects into ICF implosion modeling using multi-physics codes, improving their predictive capability.
To achieve this goal, dedicated LPI experimental platforms are valuable since they provide insight into LPI behavior and allow us to validate our LPI modeling approaches. Historically, these platforms have enabled an understanding of the essential nonlinear behavior of LPI in isolated laser speckles, how laser speckles interact with one another, and LPI processes in more complex settings.8–18
Recent CBET experiments at the OMEGA facility19,20 at the Laboratory for Laser Energetics (LLE) were designed to provide a simplified, well-diagnosed CBET platform to study the interaction of a solitary seed laser beam with up to four pump beams. By isolating the CBET interaction, this platform enables detailed comparison with fully self-consistent, nonlinear, particle-in-cell (PIC) modeling. A recent study of this type led to an improved understanding of nonlinear CBET phenomena.21 Using two-dimensional (2D) VPIC22–25 simulations, Nguyen et al. found that CBET can saturate through ion trapping-induced detuning, IAW nonlinearity, and collisional ion heating. In particular, ion trapping and heating were found to result in the peak of the plasma response function shifting to off-resonant frequencies, which reduced the CBET gain. More broadly, ion trapping was found to be one of the key CBET saturation mechanisms in indirect drive settings as well.26–29
Another nonlinearity that can affect CBET is laser-beam self-focusing. One form of self-focusing is ponderomotive self-focusing, arising when plasma density depressions form in high-intensity beams as a result of plasma hydrodynamics response to the ponderomotive pressure.7 The other is thermal self-focusing from plasma heating by the inverse bremsstrahlung.30–37 CBET in the presence of speckled-laser beam self-focusing has been examined in studies using a fluid model of the plasma.38,39 It was found that ponderomotive self-focusing in speckled beams has a significant impact on CBET resulting from beam bending and plasma-induced smoothing. However, the hydrodynamic equations employed an isothermal approximation and, thus, thermal effects37 in laser self-focusing were not included. In Ref. 21, ion–ion collisions were used to examine ion trapping and heating with a high intensity seed beam, but thermal effects were not included since the simulations were performed without electron collisions and, hence, the inverse bremsstrahlung process. Similarly, VPIC simulations of CBET driven by large-scale speckled beams28,29 captured ponderomotive effects but the parameters were relevant to indirect drive He+2 plasmas, where collisional effects and, thus, thermal effects are unimportant.
To date, the influence of ion trapping and thermal effects driven by the laser beam on CBET dynamics have not been modeled together. In this paper, we explore the nonlinear dynamics and saturation of CBET using two- and three-dimensional (3D) VPIC simulations of a seed beam interacting with a pump beam in parameter settings relevant to the OMEGA CBET experiments, namely, in a mid-Z plasma with an intense pump laser beam. To the best of our knowledge, these are the first ever multi-speckled 3D fully kinetic simulations of CBET. We focus here on the influences of ion trapping and laser-induced thermal effects. In Sec. II, we describe the PIC simulation setup for modeling the OMEGA CBET experiments, including the methods for modeling the super-Gaussian electron distribution, inter-particle collisions, and thermal effects. In Sec. III, the effects of ion trapping, collisional de-trapping, and plasma flow induced by thermal effects are discussed. Results at low seed beam intensity near the onset threshold for ion trapping are contrasted with those at high seed beam intensity with strong ion trapping. Sensitivity of simulation results to collisions and dimensionality are examined in detail by comparing collisionless, only ion–ion collision and full-collisions simulations in 2D and 3D. The main results from this work are summarized in Sec. IV.
II. MODELING OF OMEGA CBET EXPERIMENTS WITH PIC SIMULATIONS
A. Simulation setup
The plasma and laser conditions for the VPIC simulations were chosen to correspond to those of recent CBET experiments19,20 using the Tunable OMEGA Port-9 laser.40 The platform employs a gas-jet target comprising 45% nitrogen and 55% hydrogen in a plume heated by UV laser beams. To avoid overlapping the heater and pump/probe beams, the pump beams (with characteristic overlapping intensity W/cm2) and probe beam (with intensities ranging from 0.1 to W/cm2) were turned on 800 ps into the experiment. These beams had wavelengths of 351.11 and 351.40 nm and crossed at an angle over an interaction length of 220 μm. The pump and probe beams were linearly co-polarized. Using Thomson-scattering measurements, the plasma was observed to be spatially uniform over a 1.5 mm diameter plateau at initial temperatures Te = 450 eV, , and density cm−3, though the density was observed to decrease slowly to cm−3 over a period of 700 ps.
In the VPIC simulations, the plasma has electrons and pre-ionized hydrogen and nitrogen ion species H+ and N7+ (with realistic mass ratios) at 55% and 45% number density, respectively, with initial Maxwellian ion and non-Maxwellian electron (super-Gaussian of order three) distributions, matching the experiments.19,20 We model the experiments with a spatially uniform, representative density cm−3 with (where is the critical density) and electron-to-ion temperature ratio and Te = 600 eV. Refluxing boundary conditions are used for the particles: when a particle encounters a boundary, it is destroyed and a new particle is injected with velocity sampled randomly from a Maxwellian distribution for the ions and a super-Gaussian of order three for the electrons at initial temperatures Ti and Te.
The simulations are performed in two spatial dimensions (2D) in the (x, z) plane and in three spatial dimensions (3D) with cell sizes nearly equal to the Debye length λD (0.02 μm) and with time step where is the electron plasma frequency. We use 512 and 80 particles per cell per species in 2D and 3D, respectively, for collisionless simulations. In simulations with a binary collision operator, we also use 512 and 80 electron macro-particles per cell in 2D and 3D. The number of ion macro-particles per cell are chosen so that the statistical weight of the individual ion macro-particles is equal to that of the electron macro-particles, ensuring strict energy and momentum conservation in the collision operator.
The multi-speckled laser beams have an average intensity and a flat-top temporal pulse shape with rise time of 0.25 ps to . (This rise time gives the same instability growth as that from the use of a longer rise time.) is calculated on the injection boundary using the intensity values after the rise time. The laser pulse is launched from the boundary with a polarization in the y direction with the field Ey specified in a manner that approximates a Gaussian random field.41–44 In vacuum, this creates a random distribution of F/6.7 speckles with characteristic width μm and length μm. The 2D simulation geometry is shown in Fig. 1, illustrating the pump beam, seed beam, and IAW propagation directions. First-order Higdon absorbing boundary conditions45 are used for the fields.
The CBET behavior is explored in simulations with pump laser beam intensity W/cm2 (as in Ref. ) and either a low seed intensity W/cm2 or a high seed intensity W/cm2 with beam crossing angle . The beam diameter is 68 μm in 2D and 34 μm in 3D. Except for the reduced beam diameters, the laser and plasma conditions in simulations match those measured in the experiments.19,20 Under these laser and initial plasma conditions, the time evolution of CBET is modeled until the gain approaches a steady-state value (except for collisionless simulations at low seed intensity, which were terminated when CBET saturation was clearly seen but before the gain reached a steady-state value). The behavior of CBET gain as a function of the evolving plasma density measured in the experiments is not modeled because the simulations do not include the hydrodynamics of the plasma expansion.
B. Super-Gaussian electrons
C. Collisions and pump-beam-induced thermal effects
Since Coulomb collision rates scale dominantly as , where Zi and Zj are the charges of the particles, collisions are important for modeling the OMEGA CBET experiments where the plasma has mid-Z N7+ ions. The collisional simulations are performed using a binary collision model50 and can include all self- and cross-species collisions or only selected species collisions. The collision operator uses a fixed Coulomb logarithm for all collisions.51 To ensure accuracy of collision sampling,52 the collision operator is applied every N time steps, where N is defined by the condition , is the sub-cycling time step of the collision operator, dt is the simulation time step, and is the collision rate for species i and j (including self-collisions for i = j and cross-species collisions for ). The collision model has been validated by comparisons with classical transport theory53,54 and has been used successfully in prior studies of LPI experiments.19,21,41,55,56
To understand collisional effects on ion trapping and CBET dynamics, we first examine the collision rates for various scattering processes.51 The dominant rates are from N7+, H+, and electron scattering with a background of N7+. Figure 2 shows rates of frictional slowing down νs (blue curves), perpendicular diffusion (green), and parallel diffusion (red) as a function of test particle velocity for N7+ (a), H+ (b), and electrons (c) propagating in a background of N7+ (with temperature 150 eV and number density cm−3). The vertical dashed lines indicate the phase velocity of the resonant CBET IAW expressed in units of the test particle thermal speed. This provides a reference for the region of the velocity space distribution where wave-particle interactions occur in VPIC simulations of CBET. Ion trapping modifies the velocity space distribution, IAW damping rate, and CBET gain. Collisional scattering may significantly affect the velocity space distribution in regions around the IAW phase velocity and, as a result, the IAW damping rate and the CBET gain.
For test N7+ ions, shown in Fig. 2(a), the slowing-down rate is comparable to the perpendicular diffusion rate, and both rates are much higher than the parallel diffusion rate in the region around the IAW phase velocity. The estimated timescale for slowing-down and perpendicular diffusion is about 10 ps, shorter than a typical VPIC CBET simulation time ( ps). In CBET, IAW can grow to a large amplitude and trap ions, forming a trapped particle population and modifying the initial Maxwellian distribution function. Collisional scattering processes (mainly slowing-down and perpendicular diffusion) can modify the trapped population and restore the distribution toward a Maxwellian form. The effects of collisions on N7+ ions trapped by the IAW can lead to significant N7+ ion heating (see Fig. 5 in Ref. 21) as a result of the strong self-collisions.
For test H+ ions, shown in Fig. 2(b), the perpendicular diffusion rate is greater than the slowing-down rate, and both rates are much higher than the parallel diffusion rate in the region around the IAW phase velocity. The estimated timescale for perpendicular diffusion and slowing-down is shorter than 10 ps. The IAW energy is first converted to the directed kinetic energy of the trapped ions. Collisions convert this directed energy into ion thermal energy. However, the partition of IAW energy that gets transferred to each ion species cannot be determined from theory estimates. Nevertheless, we observe from simulations that the H+ ions heat less than the N7+ ions (also see Fig. 5 in Ref. 21).
As a result of collisions, the ion distributions in the VPIC simulations show a reduced population of trapped ions and increased ion temperatures as time increases. Both the simulated timescale and the amount of ion collisional heating are consistent with experimental observations.19–21 Since N7+ and H+ ion scattering rates on background electrons are much lower, the effects of collisions on ion trapping can be represented in VPIC simulations using only ion–ion collisions as done in Ref. 21.
The CBET IAW may interact with electrons that have a velocity much smaller than the electron thermal speed. For test electrons, shown in Fig. 2(c), the scattering rates increase sharply as the test electron velocity decreases. For the parameters considered here, the main effects that electron collisions have on CBET are fluctuations of both electron temperature and plasma flow.
To understand the thermal effects (i.e., transport effects resulting from electron collisions) on CBET gain, we first consider the dynamics of a high-intensity pump beam in the absence of a seed beam by exploring the density response, thermal fluctuations , and thermal-effects-induced flow in the high-intensity pump beam with and without collisions. In Fig. 3, we show the plasma density, flow vx, and temperature profiles across the pump beam from simulations at pump intensity W/cm2. The laser field is given in frame (a) and frames (b) to (e) show the profiles along x of charge density, flow velocity vx, and temperature from the simulation with full collisions for the electrons (black curves), H+ (blue), and N7+ (red) ions at an early time and a late time that is long compared with the collision times. Similarly, vx and temperature profiles from the collisionless simulation are shown in frames (f) and (h). The laser light drives perturbations in density, flow velocity, and electron temperature in both the collisional and collisionless simulations. The two dominant transverse gradient scales in these perturbations are the characteristic width of the speckles and the beam width. At early times, the flow perturbations at the speckle scale are comparable in the collisional and collisionless simulations. Later, these small-scale flow perturbations evolve to a larger scale comparable to the beam width. Contrasting the results in frames (c)–(e) with those in (f)–(h), we find enhanced flucturations from thermal effects, as shown by the differences in vx and Te with and without collisions. With full collisions, fluctuations of electron density ne and temperature Te are 10% and 25% of their initial values, respectively. Additionally, thermal effects induce a plasma flow away from the center of the pump beam with a peak vx value ∼5% (for the N7+ ions) to 10% (for the H+ ions) of the IAW phase velocity . In contrast, in the collisionless case, the fluctuations of electron temperature Te are smaller with peak values at 10% and the plasma flows are also reduced by 2× without thermal effects.
If we were to add a seed beam with a frequency initialized to the resonance condition for CBET in the absence of flow, these modifications to the plasma conditions could affect the CBET matching conditions as well as the CBET IAW damping rates through increased Te. The impact of thermal-effects-induced flow on the CBET gain can be understood from linear theory: The energy transfer is proportional to the imaginary part of the plasma response function.58,59 These flows can Doppler-shift the CBET IAW frequency by a few percent, changing the peak of the plasma response function from being on resonance to being off resonance. At the speckle scale, this leads to up-shifted frequencies in the left half of the speckles and down-shifted frequencies in the right half of the speckles. At the laser beam scale, this leads to up-shifted frequencies in the left half of the laser beam and down-shifted frequencies in the right half of the laser beam. In the presence of ion trapping, such a Doppler-shift away from resonance can be comparable to that caused by trapping-modified distributions, as reported in prior studies of CBET.21,29 In Sec. III below, we examine the combined detuning effects from ion trapping and flow on the nonlinear saturation of CBET.
III. NONLINEAR SATURATION OF CBET BY ION TRAPPING AND THERMAL-EFFECTS-INDUCED FLOW
CBET dynamics can be influenced by ion trapping and ion-heating-induced detuning,21 as well as by flow induced by thermal effects as explained above. Collisions may affect both processes: Collisional de-trapping is important near the onset threshold intensity for trapping; collisions also introduce thermal effects, enhancing the plasma flow. Since ion-electron collision rates are small compared with Ion–ion collision rates (see Sec. II C), only ion–ion collisions were included in prior CBET simulations21 to isolate the effects of ion trapping and heating from those of evolving plasma conditions. Here VPIC simulations were run with full collisions, ion–ion collisions, and no collisions to determine the effect of collisions on ion-trapping and flow modifications in CBET.
In this work, the CBET gain percentage of the seed beam is defined as , where and are the input (measured on the right boundary) and output power (measured on the left boundary) of the seed beam, respectively. The gain percentage can be converted to gain as used in linear theory.
FSRS may be unstable in the CBET amplified seed beam, as found in VPIC simulations relevant to indirect-drive.28,29 However, FSRS is found to be unimportant in this work because of the low plasma densities used in the OMEGA CBET experiments. In the absence of FSRS, ion trapping is the main nonlinear effect aside from detuning via flow fluctuations and pump depletion.
We start with a discussion of results from 2D simulations, followed by 3D simulations that confirm the main results.
A. 2D VPIC simulations
In Fig. 4, the CBET gain vs time from collisionless, ion–ion collision, and fully collisional 2D simulations are shown at two seed intensities: low seed intensity W/cm2 (red curves) and high seed intensity W/cm2 (black, blue, and green curves). Simulations using different speckle patterns show that the CBET gains are insensitive to speckle statistics in the cases examined. At low seed beam intensity, the CBET gain is sensitive to collisions as evidenced in the substantially different CBET gains shown by the three red curves. In contrast, at high seed intensity, the CBET gain is insensitive to collisions shown by the black, blue, and green curves.
Differences in the CBET gain can arise from the following: (1) Full collisions lead to changes in ne and Te and modify the plasma conditions for CBET resonance. An increase in Te can also lead to a reduction of IAW damping. (2) Full collisions enhance laser-induced flow, as shown in Fig. 3, which also modifies the CBET resonance. (3) Collisional de-trapping has a significant effect on the CBET gain when ion trapping is relatively weak.
First, we use the low seed intensity cases to explain the effects of collisional de-trapping of ions. In Fig. 5, at low seed intensity W/cm2, the characteristic ion 2D velocity space distributions are compared from collisionless, ion–ion collision, and fully collisional simulations shown in frames (a1) and (b1), (a2) and (b2), and (a3) and (b3), respectively. These distributions are sampled at a time (t = 90 ps) during CBET saturation but are generally representative of the entire durations of the simulations. The trapping tail, which indicates the IAW propagation direction (see Fig. 1), is the strongest for the collisionless case in Figs. 5(a1) and 5(b1). In simulations with collisions, collisional de-trapping and ion heating resulting from collisional scattering can effectively reduce the trapping tail in the velocity space distributions. As explained in Sec. II C, the N7+ and H+ ion scattering rates on background electrons are much lower than ion–ion scattering rates. Thus, the effects of collisions on ion trapping are similar in simulations with ion–ion collisions and with full collisions, as seen by comparing Figs. 5(a2) and 5(b2) to Figs. 5(a3) and 5(b3). The CBET gain reduction from the solid red curve to the dashed red curve in Fig. 4 results from ion–ion collisional de-trapping alone since the ponderomotive-effects-induced flow is the same in collisionless and ion–ion collision simulations. We note that while electron trapping in the IAW and its induced nonlinearity (i.e., a positive nonlinear frequency shift as opposed to the negative frequency shift associated with ion trapping) have been reported with large ,60,61 they do not present discernible effects in our CBET simulations.
In Fig. 6, at low seed intensity W/cm2, 2D simulation results for flow velocity vx and temperature are compared for collisionless [frames (a1)–(d1)], ion–ion collision [frames (a2)–(d2)], and fully collisional [frames (a3)–(d3)] cases. As we explained in Sec. II C, laser-driven flow at the speckle scale can begin to detune the CBET resonance at early times of the interaction. Here, we focus on flow induced by thermal effects at long time scales. Profiles of vx and temperature along x (at z = 0.0 μm) for electrons (black curves), H+ (blue) and N7+ (red) ions are shown at t = 50.0 ps (near the peak of CBET gain) and at 91.5 ps (during saturation). Note that in the presence of ion trapping, taking the first moment of the distribution function yields an apparent flow with all positive values in velocity even though there is no initial flow in the simulation. Furthermore, taking the second moment of the distribution function yields an effective temperature, which is higher than the initial ion temperature. Both arise from trapping as seen in Fig. 5. In the collisionless and ion–ion collision simulations, the ponderomotive-effects-induced flow is small and the ion velocity is predominantly from the modification of the velocity space distribution due to ion-trapping. The effects of de-trapping can be discerned from and the ion temperature profiles by comparing the collisionless and the ion–ion collision simulations: with ion–ion collisional de-trapping, the amplitudes of vx and ion temperatures are lower in the ion–ion collision case.
Next, we examine how flow induced by thermal effects accounts for the reduction in the CBET gain from the dashed red curve to the dotted red curve in Fig. 4. As shown in Fig. 3, thermal effects in the intense pump beam induce a plasma flow outward from the axis of the pump beam with a peak value ∼5% (for the N7+ ions) to 10% (for the H+ ions) of the IAW phase velocity . In the simulation with full collisions, as a result of the combined effects from ion trapping, collisional ion heating, and thermal-effects-induced flow in the pump beam, the profile of is asymmetrical about the center of x-domain (x = 69 μm), i.e., the magnitude of the flow is larger to the right of x = 69 μm than to the left. This asymmetry decreases in time as trapped ions in the tails of the distributions go through slowing-down, perpendicular diffusion, and thermalizing processes, leading to a flow resembling that in the pump-only simulation in Fig. 3. This induced flow further detunes CBET away from resonance in addition to the ion trapping/heating induced detuning, leading to the CBET gain reduction from the dashed red curve to the dotted red curve in Fig. 4. We note that in the simulation with full collisions, there is also heating in Te up to 15%.
At the high seed intensity W/cm2, in Fig. 7, 2D simulation results for flow velocity vx and temperature are compared for ion–ion collision [frames (a1) and (b1)] and full collision [frames (a2) and (b2)] cases. The ion trapping is strong and leads to nearly all positive values of and large ion temperature increases. The ion trapping induced detuning dominates over effects from collisional de-trapping and flow. These results explain why the CBET gain at high seed intensity is insensitive to collisional de-trapping and thermal-effects-induced flow, as shown by the lower black solid, blue dashed, and green dotted curves in Fig. 4. Ion trapping induced saturation leads to a reduction of the CBET gain to a level that is much lower than seen at low seed intensity, as shown in Fig. 4, or as predicted by linear theory including pump depletion.19,21
Note that for the 2D simulations, the average intensity for the speckled pump and seed beams is defined , as used in Omega experiments,62 where brackets denote spatial averages over the beam cross section at the laser launch surface. Converting the intensity to that based on the normal average intensity definition and using simple linear theory estimates,21 we find that the gain percentage from linear theory is % (using the interaction length where d is the beam diameter and θ is the beam crossing angle). This value is below those from simulations at low seed intensity but above those from simulations at high seed intensity.
B. 3D VPIC simulations
To determine the effects of dimensionality and to confirm the main findings from 2D simulations, additional 3D simulations are performed with and without collisions at low and high seed intensities. The simulation domain is 70.4 μm3 in x, y, z with periodic boundary conditions in y. The pump and probe laser beams on the x and z launch boundaries have rectangular cross sections of dimension μm2 in their respective launch planes. With more than 5 × 109 cells and 80 particles per cell, the simulations employ 211 600 MPI ranks and over 1.2 × 1012 particles, which reaches the limit of what is practical with the computing resources available to us. However, at reduced scale, for a random distribution of F/6.7 speckles with a characteristic width μm and length μm, these simulations represent the multi-speckled nature of the 3D CBET interactions. The 3D geometry allows for additional laser-driven flow and temperature fluctuations in y as well as modifications to ion trapping such as enhanced side loss compared with the 2D simulations presented in Sec. III A.
The time histories of the CBET gain from the 3D simulations at low seed intensity W/cm2 and high seed intensity W/cm2 are shown in Fig. 8(a). (Note that here, in contrast with 2D simulations, we use the more traditional definition of intensity , where brackets denote a spatial average over the beam cross section at the launch plane.) The red solid curve and the red dotted curve are at low seed intensity without collisions and with full collisions, respectively. Similar to the 2D simulation results shown in Fig. 4, the 3D simulations also show that at low seed beam intensity, the CBET gain is sensitive to collisions. This sensitivity, as explained using results from the 2D simulations, arises from collisional de-trapping and thermal-effects-induced flow, as well as from changes in the plasma conditions in the presence of full collisions.
To compare with the thermal effects in 2D, results from the 3D simulation with full collisions at low seed intensity are shown at a time t = 64.2 ps during CBET saturation in Figs. 8(b)–8(d). Profiles along x (averaged over 10 cells around ) for charge density (b), flow velocity vx (normalized to the IAW phase velocity ) (c), and temperature (normalized to the initial temperature for each species) (d) are shown for the electrons (black curves), H+ (blue) and N7+ (red) ions. Thermal effects in the intense pump beam are stronger in 3D, as evidenced by the stronger heating of Te in 3D with an increase in Te up to 30%, in contrast with that in Fig. 6(d3). The asymmetry in the profile of about the center of the x-domain (x = 35 μm) reflects the presence of induced flow superimposed onto an all-positive valued profile due to ion trapping alone. The net profile has a larger magnitude in the x > 35 μm region than in the x < 35 μm region.
For high seed intensity, the time histories of CBET gain from 3D simulations are shown in Fig. 8(a) by the black solid and green dotted curves without collisions and with full collisions, respectively. Profiles along x for flow velocity vx (e) and temperature (f) are also shown for the different species. The strong ion trapping effects dominate over those from thermal effects, and the values of are largely positive. In contrast to the low seed intensity case, at high seed intensity, the CBET gain is insensitive to collisions as well as dimensionality as a result of the dominant ion trapping effects.
Though 2D and 3D simulations of CBET exhibit similar behavior in many regards, there are also important differences. In 3D, the most intense laser speckles in the crossing beams, for which CBET IAW amplitudes are largest, are less likely to overlap than in 2D, where essentially every intense speckle in the pump and seed beams overlap. One would also expect weaker transverse coupling among neighboring speckles through the exchange of side-loss ions as compared with 2D. The side-loss ions originate from the trapped ions that escape the narrow IAW potential well63 at a rate . The side-loss ions can lower the Landau damping rate for IAW as they enter neighboring speckles, thereby increasing the wave coherence length and enhancing energy transfer.28 In 2D, every side-loss ion leaving a given speckle will encounter a neighboring speckle (except for speckles at the edge of a beam). In 3D, only a fraction of the side-loss ions leaving a given speckle encounter a neighboring speckle, which reduces the transverse coupling among neighboring speckles.
Another interesting result of our study is the apparent effects of boundary conditions on the CBET gain at low seed intensity. Comparing the red solid and dotted curves in Fig. 4 (with larger beam diameter) to their counterparts in Fig. 8(a) (with smaller beam diameter), we see that although the peak of the gain curve from the 3D collisionless simulation is lower than that from its 2D counterpart, the 3D gain saturates at a higher level than that of the 2D simulation; for fully collisional simulations, the gain from 3D is also higher than that from 2D. This is because in the 2D simulations, all physical boundaries are open to waves and particles. As a result, energy from the laser beams enters and exits the boundaries and particle and wave energy from the CBET processes also exit the boundaries, which act as an energy sink. In particular, ions accelerated by wave-particle interactions and electrons heated by the laser can propagate out of the boundaries and are replenished by Maxwellian and super-Gaussian ions and electrons at initial temperatures flowing into the system from the boundaries. The 3D simulations allow for physical processes, such as ion trapping and laser-driven flow and temperature fluctuations in the y dimension, but the boundary conditions are periodic in y. Ions accelerated by wave-particle interactions and electrons heated by the laser are wrapped around in the y domain, which is effectively infinite. It is unclear in this case what the correct interaction length is in estimating the gain percentage using simple linear theory.21
IV. SUMMARY AND CONCLUSIONS
Using 2D and 3D VPIC simulations for parameters relevant to recent CBET experiments at the OMEGA laser facility, the influences on CBET of ion trapping in ion acoustic waves and laser-driven flow and temperature fluctuations have been examined. In the mid-Z plasma with an intense pump laser beam, ion trapping, collisional de-trapping, and plasma flow induced by thermal effects all have important effects on CBET. Ion trapping can enhance or detune the CBET resonance. At low seed beam intensity near the onset threshold for ion trapping, collisional de-trapping has more significant effects on CBET and thermal effects can induce flows, which also detune the CBET resonance at a level comparable to that from trapping. At high seed beam intensity, ion trapping effects on CBET dominate over those from collisional de-trapping and the flow induced by thermal effects. The trapping dominated CBET saturation has been confirmed from simulations at seed intensities ranging from 3 to W/cm2. As a consequence, we found that CBET gain is sensitive to collisions and dimensionality at low seed beam intensity where ion trapping is weak but is insensitive to collisions and dimensionality at high seed beam intensity with strong ion trapping. In both the high- and low-intensity seed cases, the amount of energy gained by the seed beam represents a very small fraction of the pump and the effects of pump depletion are minimal.
We note that the pump intensity used in this study ( W/cm2) is higher than the intensity of the individual pump beams of the experiments, though matches the pump intensity used in our prior CBET modeling study.21 This choice allows us to compare the simulation results and isolate the thermal effects in the presence of ion trapping. These effects include the introduction of a spatially varying (over lengths comparable to the width of a laser speckle or the beam) transverse flow whose speed is a small fraction of the IAW speed. Such a flow may be difficult to measure experimentally, yet may nevertheless affect CBET gain. An additional pump-only simulation with full collisions at a lower pump intensity W/cm2 has also been performed and the plasma flow vx profiles are shown in Fig. 3(d) by the orange and green curves for N7+ and H+ ions, respectively, with the peak vx values ∼4% to 5% of the IAW phase velocity (compared to the ∼5% to 10% at W/cm2 that can Doppler-shift the CBET IAW frequency by a few percent and detune the peak of the plasma response function from being on resonance to being off resonance). The results suggest that the intensity threshold for thermal effects to affect CBET gain is near W/cm2. For these simulations, the average intensity for the speckled beam defined as are used. Converting the intensity to that based on the normal average intensity definition , the pump beam intensities for these simulations are and W/cm2.
The electron temperature increases at 25%–30% resulting from laser heating shown in Figs. 3 and 8 occur over a time interval that is long compared to collision times but short compared to the 50 ps temporal resolution of the Thomson-scattering system measurements.19 Together with the cooling associated with the slowly decreasing density observed in the experiments, the electron heating may be difficult to observe experimentally.
One of the main results of this work is the elucidation of important effects of ion trapping on CBET saturation. This has broad applicability to fusion-relevant laser-plasma configurations as discussed in prior work.26,27 Collisional de-trapping and plasma flow induced by thermal effects also have wide applicability to ICF and HED settings. The choice of parameter regime for the simulations (and the Top9 experiments) has relevance to other ICF and HED settings as well. The plasma composition and used in this work resembles that of CH target in direct-drive ICF implosion. Also, in the low plasma density regime, as in the case for low gas fill and near-vacuum hohlraums,66–70 FSRS and BSBS are not dominant instabilities, so ion trapping is the main nonlinear effect aside from detuning via flow fluctuations and pump depletion, and its effects can be examined in isolation. Key parameters determining the CBET gain include , IAW damping rate, and IAW (or IAW phase velocity normalized to the ion thermal velocity).29 Although in the regime of low density, the simulations in this work (and the Top9 experiments) can access different regimes of IAW and IAW damping rates by varying the crossing angle and density. The resulting values of the IAW damping rates are relevant to those in OMEGA and NIF scale direct-drive experiments.21 Moreover, this work reports on key processes in the nonlinear optics of plasmas, identified as a grand challenge problem in plasma and high energy density physics. As such, these results would be of broad interest to the plasma physics community. For example, results from recent beam combiner experiments indicate that the energy output in the seed beam may be increased further by a reduction in the plasma density64,65 if the nonlinear saturation of CBET were mitigated.
We further note that the kinetic simulations for our study were performed under initial conditions of no bulk plasma flow but with a frequency difference between the pump and seed beams giving rise to the excitation of CBET. In the presence of plasma flow, CBET can also occur between beams of the same frequency in regions of sonic flow.38,39,71 Under these conditions, thermal-effects-induced flow would be expected to have similar detuning effects on CBET, a point we intend to consider in future work.
This work was performed under the auspices of the U.S. Dept. of Energy by the Triad National Security, LLC Los Alamos National Laboratory and was supported by the LANL Office of Experimental Sciences Inertial Confinement Fusion program. LLE work was supported by the Department of Energy National Nuclear Security Administration under Award No. DE-NA0003856. The authors acknowledge valuable discussions with Dr. Harvey Rose. VPIC simulations were run on ASC Trinity supercomputer under ATCC and the Large Scale Calculations Initiative (LSCI).
Conflict of Interest
The authors have no conflicts to disclose.
Lin Yin: Conceptualization (lead); Data curation (lead); Formal analysis (lead); Investigation (lead); Methodology (lead); Software (equal); Visualization (lead); Writing – original draft (lead); Writing – review & editing (lead). K. L. Nguyen: Investigation (supporting); Writing – review & editing (supporting). Brian J. Albright: Conceptualization (supporting); Investigation (supporting); Methodology (equal); Software (equal); Writing – review & editing (supporting). Alexander G. Seaton: Investigation (supporting); Methodology (supporting); Software (equal); Writing – review & editing (supporting). Aaron Hansen: Investigation (supporting); Writing – review & editing (supporting). Dustin Froula: Investigation (supporting); Writing – review & editing (supporting). David Turnbull: Investigation (supporting); Writing – review & editing (supporting). John Patrick Palastro: Investigation (supporting); Writing – review & editing (supporting).
The data that support the findings of this study are available within the article.
APPENDIX: SUPER-GAUSSIAN ELECTRONS
Generate four samples p1, p2, p3, and p4 independently from a uniform random distribution from zero to one.
Pick a sample speed , obtained from inverting .
If , accept the sample and proceed to step 4. Otherwise, reject the sample and return to step 1.
- Generate the components of u by sampling uniformly over solid angle
Generate two samples p1 and p2 independently from a uniform random distribution from zero to one.
Pick a sample u value , obtained from inverting .
If , reject the sample and return to 1.
This prescription rejects approximately one third of samples and has proven to be acceptably efficient.
Generate three samples p1, p2, and p3 independently from a uniform random distribution from zero to one.
Pick a sample u value , obtained from inverting .
If , reject the sample and return to 1.
- Set and set the two tangential components to
This sampling rejects a fraction 0.41 of trials on average. For numerical stability, it is recommended that uz be rejected above uz = 4 prior to entering this portion of the sampling algorithm. Overall, for VPIC calculations on a variety of many-core supercomputers available to us, the Maxwellian reinjection of super-Gaussian distributions adds of order 10%–20% overhead to the runtime, making these runs practicable.