Resonant upconversion through a sixth order relativistic nonlinearity resulting in a unique resonance was recently proposed [Malkin and Fisch, Phys. Rev. E 108, 045208 (2023)]. The high order resonance is a unique non-integer multiple of a driving pump frequency resulting in a frequency upshift by a factor of . We demonstrate the presence, unique requirements, and growth of this mode numerically. Through tuning waves to high amplitude, in a mildly underdense plasma, the six-photon process may grow more than other non-resonant but lower order processes. The growth of the high frequency mode remains below the nonlinear growth regime. However, extending current numerical results to more strongly coupled resonances with longer pulse propagation distances suggests a pathway to significant upconversion.
I. INTRODUCTION
The falloff in available laser power below the optical is severe,1 and producing high intensity UV, XUV, or x-ray pulses remains a challenge with current electron beam driven2,3 or laser driven sources.1,4–6 Recent proposals detail using resonant interactions in plasma to upconvert high fluence sources, such as the National Ignition Facility,7 to shorter wavelength.8–11 This paper details the first numerical test of higher order plasma mediated wave mixing for upconversion.12
Wave mixing builds off of a long history of using plasmas components to replicate the behavior of well-known solid state parts, but which operate at higher laser intensity. Pulse amplification can be done with electron plasma waves,13–19 ion acoustic waves,20–22 or magnetized waves.23 Ion acoustic waves might also enable beam combination.24–26 Plasmas could serve as polarizers,27,28 q-plates,18 lenses,29 or gratings.30
Mildly relativistic wave mixing to upconvert optical light with high efficiency to much shorter wavelengths is a more recent proposal. Pump photons are combined into a higher frequency output through nonlinear coupling. The nonlinear coupling arises from corrections to the plasma frequency at significant electron Lorentz factor.8 This coupling contrasts with previous work in that it can be resonant, without requiring a modulated plasma density.31 The lowest order nonlinearity results in four photon coupling. A simple colinear alignment between pump and seed beams requires a complex set of resonance conditions,9,10 and non-colinear arrangements result in difficulties for amplification.11 A higher order, six-photon, process addresses both of these problems.12 This six-photon scattering process consumes four input photons with frequency ω to produce a pair of photons at frequencies of approximately and . This unique non-integer multiple frequency shift from the six-photon process could allow for large changes in wavelength if it were to be cascaded multiple times.
We present a numerical investigation, which demonstrates the unique six-photon resonance. In contrast to previous work, here the resonance is not assumed in a reduced (slowly varying envelope) and expanded model. Instead, we consider a model, equivalent to a cold relativistic electron fluid, in which the full wave dynamics are evolved without the resonance embedded. As the resonance is of higher order, coupling must be maximized to generate observable effects. Thus, we consider the dynamics without expanding in either normalized plasma density, , for laser frequency ω and plasma frequency , or normalized wave amplitude, , for laser vector potential A. When both parameters are not , which is required to achieve significant coupling, there may be large errors in any expansion of nonlinear terms to finite order. Numerical experimentation is used to examine the fundamental physics at play and validate the viability of relativistic six-wave mixing for upconversion.
The exposition of our numerical results on relativistic six-wave mixing in plasmas is outlined as follows. In Sec. II, we review the dynamics and resonances available for the coupling between mildly relativistic waves, a < 1, in an underdense plasma. In Sec. III, we describe the choice of the particular numerical tools used to examine this regime. In Sec. IV, we demonstrate the unique resonance proposed in by Malkin and Fisch. In Sec. V, we discuss the applicability of this work to the aim of achieving significant upconversion, review our work, and discuss possible avenues for continued exploration.
II. BACKGROUND
The dynamical equations contain multiple sources of scattering. Raman scattering, where there is coupling between two electromagnetic waves and an electron plasma wave, is the lowest order process. When terms are limited to the second order, fluctuations of density are captured by the electrostatic potential in Eq. (1), which are driven by the ponderomotive force, derived from combining Eqs. (2) and (3) at the lowest order. For the purpose of this paper, Raman scattering is purely parasitic, and forward Raman scattering may grow of off seeded fluctuations in the plasma to move energy into waves, which are no longer resonant for upconversion. The higher order processes responsible for coupling, which may result in upconversion, occur from the expansion of the electron Lorentz factor . The first term in this expansion gives the four-photon coupling,8 and the next provides the six-photon coupling.12 These resonances are unique in that they result from the full consideration of the electron dynamics without averaging the ponderomotive force. The subtleties between different resonance regimes and the precise scaling of the coupling in these regimes are not discussed here, as the particular expressions can be extensive and are best understood in a more thorough treatment.
The parameters to produce a validation of the resonance outline by Eq. (5) are restrictive. To validate the non-integer resonance, without assuming it a priori requires a full wave model. This full wave model must have a high resolution to resolve the high frequency mode. However, amplification distances resulting from the weak coupling given by Eq. (7) result in significant plasma and pulse lengths, which are difficult to resolve. To maximize the coupling, we use Eq. (7) to guide the parameter selection toward a high intensity a0 and high plasma density n. The more highly resolved model allows us to work in a regime where higher order corrections to the dispersion relation and to the equations of motion are more accurately captured than in a nonlinear series expansions or in an envelope based simulation.
III. NUMERICAL METHODOLOGY
Simulations are initialized to test for the growth of the proposed six-photon scattering when waves are evolved by Eqs. (8)–(10). In a homogeneous plasma of density of n, we initialize a strong “pump” wave packet with and a “seed” wave packet of varying k2. An example of the initial conditions can be seen in Fig. 1. The seed frequency is given a weak amplitude and, to account for different group velocities, is initially leading the pump wave. In the ideal case, the pump wave is completely depleted, and all of the energy is in the low and high frequency waves post interaction. The initial conditions for a, e, and l, are chosen such that the time derivative of each field are initially zero to the lowest non-zero order. Through examination of the spectrum, and in particular, the mode at , we demonstrate that Eq. (5) describes a growing mode, which is only seen when it is satisfied.
The pump wave is initialized with a super-Gaussian envelope, which is led by a much weaker low frequency seed. If pump depletion could be reached, the ideal outcome would be such that the high frequency wave and low frequency wave are amplified until the pump is completely depleted. While at equivalent amplitude, the high frequency wave would contain the majority of the remaining energy due to having a much higher frequency.
The pump wave is initialized with a super-Gaussian envelope, which is led by a much weaker low frequency seed. If pump depletion could be reached, the ideal outcome would be such that the high frequency wave and low frequency wave are amplified until the pump is completely depleted. While at equivalent amplitude, the high frequency wave would contain the majority of the remaining energy due to having a much higher frequency.
IV. RESULTS
The presence of a uniquely resonantly growing mode can be demonstrated in the evolution of Eqs. (8)–(10). For a strong laser pulse, many different scattering channels operate coincidently resulting in the growth of a large number of discrete modes. The intensity of these modes is oscillatory, unless conditions are precisely chosen such that waves are pushed into resonance.
The nature of this resonance can be demonstrated by varying the wavenumber, k2, of a low frequency seed and examining the energy that is scattered into the expected resonance at . Both waves start in a homogenous plasma and are evolved until a minority of the seed overlaps with the pump, as diagramed in the co-propagating coordinate scheme in Fig. 1. The plasma is initialized such that . While the plasma density is such that it is relatively underdense for the pump wave, we note that the large frequency separation between the pump and the seed wave means the plasma density is not much smaller than the quarter critical density of the seed wave. The pump is given an , and the seed is overlapped with a factor of 200 lower intensity. The pulse length, which is set by the slippage condition between the pump and the seed is 450 wavelengths, which propagates for a duration of 4000 laser periods. This propagation length is much shorter than the estimated gain length of laser wavelengths, and reaching the nonlinear regime will not be possible with the current code. This duration is set such that at the end of the simulation, the low frequency seed will have slipped past half of the pump wave, at which point the growth of the high frequency mode has been arrested. The seed frequency is varied around the proposed resonance, with a minimum frequency bounded by two plasmon decay , and the maximum frequency bounded by the coupling of the third harmonic of the low frequency seed to the stokes shifted pump wave, where . Both processes at the bounds are lower order and quickly produce a large electron plasma wave, whose dynamics dominate the desired higher order upshift.
The evolution of the simulation spectrum shows the many nonlinear processes at play, but the resonant six-photon process is the most significant. The spectrum created by the wave simulation for a resonant simulation is shown in Fig. 2. For a resonant k2, the strongest growth is of the six-photon resonance; however, energy may be transferred to many other modes. For example, the second most defined new peak in the spectrum is the third harmonic of the pump wave. However, even though this mode is created by a lower order and thus more strongly coupled, it is not resonant so it cannot continually grow and is saturated at this level. The evolution of the energy in different linear combinations of the pump, seed, and plasma frequency shown in Fig. 2 demonstrates this. Other bands oscillate in energy, while the high frequency six-photon resonance grows to a level higher than that of other nonlinear processes. The growth of the six-photon resonance shows good agreement with the analytical theory, where the analytical estimate is and a fit to the growth in simulations gives a growth rate of . The growth rate peaks when there is good overlap between the pump and the seed, declining at later times as the seed spatially lags behind the pump wave. Many other peaks corresponding to combinations of the pump and seed wave occur in the spectrum, but at weaker levels. Of particular note is the four-photon process, corresponding to , which experiences some scattering, but is not resonant in the colinear geometry without further tuning of the pump and seed wave.9
The evolution of the spectrum of a pump ( ) and seed [ ]. Both waves are initialized within a homogenous plasma and evolved until a majority of the seed has slipped past the pump. Nonlinear distortions of the pump and seed are weak, but many peaks can be seen in the final spectrum compared to the initial spectrum (top). In k space, many spontaneous, but non-growing, modes are instantaneously excited across the spectrum. The particular mode of interest at may grow if Eq. (5) is satisfied allowing it to achieve higher intensities than those at modes resulting from lower order processes. This can be seen through looking at the energy in each frequency band, (bottom), where only the resonant mode is growing, and other nonlinearly driven bands oscillate in energy. The energy for each band is calculated by integrating the spectrum in a window of each peak of width .
The evolution of the spectrum of a pump ( ) and seed [ ]. Both waves are initialized within a homogenous plasma and evolved until a majority of the seed has slipped past the pump. Nonlinear distortions of the pump and seed are weak, but many peaks can be seen in the final spectrum compared to the initial spectrum (top). In k space, many spontaneous, but non-growing, modes are instantaneously excited across the spectrum. The particular mode of interest at may grow if Eq. (5) is satisfied allowing it to achieve higher intensities than those at modes resulting from lower order processes. This can be seen through looking at the energy in each frequency band, (bottom), where only the resonant mode is growing, and other nonlinearly driven bands oscillate in energy. The energy for each band is calculated by integrating the spectrum in a window of each peak of width .
Varying the seed frequency demonstrates the uniqueness of the resonance point. Figure 3 shows the energy in the post interaction spectrum for simulations at varying seed frequency, k2. Only the nonlinear growing modes are presented as the pump and seed wave are strong enough to dominate the weak nonlinear modes. The strongest nonlinearly created mode is the predicted mode. The parameter scan is restricted by other resonances in the system. As the seed frequency approaches the upper value plotted, the third harmonic of the seed approaches the Stokes line of the pump, resulting in significant Raman scattering.
A parameter scan across seed wavenumber k2 results in the varied excitation of many modes. Each scatter point corresponds to the energy in the final spectrum for the dominant nonlinear modes. The expected resonance at is more highly amplified than neighboring wavenumbers. The range of the parameter scan over k2 is heavily restricted by the growth of electron plasma waves at the bounds of k2. Shorter duration simulations highlight a larger difference in the energy of the amplified across a wider range of k2 at the cost of less amplification of the k1 wave.
A parameter scan across seed wavenumber k2 results in the varied excitation of many modes. Each scatter point corresponds to the energy in the final spectrum for the dominant nonlinear modes. The expected resonance at is more highly amplified than neighboring wavenumbers. The range of the parameter scan over k2 is heavily restricted by the growth of electron plasma waves at the bounds of k2. Shorter duration simulations highlight a larger difference in the energy of the amplified across a wider range of k2 at the cost of less amplification of the k1 wave.
Some initial small amount of energy can be transferred into the mode regardless of the seeded k2. However, this energy will oscillate if k2 is not chosen correctly, which occurs as the seed frequency is moved away from the resonance, shown by the red line. As k2 becomes closer to the resonance, the period of oscillation grows, and the energy may coherently grow. The efficiency of this process remains low, with the energy in the high frequency band more than below the initial pump energy. To reach the nonlinear growth regime, the resonance has the condition . For the parameters in the simulations presented, this corresponds to a maximal detuning of . This is stricter than the results from Fig. 3 as the simulations are not propagated for a long duration. Achieving more significant growth, and correspondingly lowering the requirements on the resonance, will require a significantly greater plasma density.
Other resonances may result in growing parasitic modes and become limiting. The desired scattering becomes dominated by Raman scattering if the seed frequency can create an electron plasma wave. Raman scattering may be seeded off of noise; however, there are two other limits also present as the seed frequency changes. If the seed frequency approaches the threshold for two plasmon decay at low k2, or the third harmonic of the seed approaches the Stokes shifted frequency of the pump, Raman scattering may become strong. At these modes, the spectra starts to produce many ωpe shifted bands, which distort the spectrum. When these two resonances occur, the electron plasma wave, and corresponding density fluctuations, quickly grow beyond the limits of the cold hydrodynamic model. We aim to avoid Raman scattering, so this limitation should not be relevant for k2 detuned from these operating points. However, at high density, these bounds may be close to the resonance where . This limits the plasma density, and thus the coupling of the six-photon coupling. Fluctuations in plasma density may push the low frequency wave into resonance and produce a large electron plasma wave. A noisier environment than the shown simulations will also provide a more fecund environment for parasitic Raman scattering. Further manipulation of the plasma and laser pulses will be needed such that Raman scattering does not dominate any desired higher order nonlinearity.
V. SUMMARY AND DISCUSSION
This paper demonstrates through numerical simulation the unique non-integer harmonic resonance driven by six-photon coupling in a plasma. This coupling, when at resonance, driven by the fifth order corrections from the electron Lorentz factor, is large enough that in the post interaction spectrum, the mode is larger than other nonlinear processes. This is true even though Raman scattering and relativistic four-wave mixing are lower order processes. It is isolated in low noise pseudospectral simulations. However, for the pulse durations and propagation distances considered, it is still below the strength of the injected waves. The non-integer multiple of the frequency is required for the six-photon scattered wave to be strong. If the frequency is detuned, growth of the high frequency mode is orders of magnitude lower, and if it is significantly detuned, the injected wave can drive other parasitic processes.
The six-photon process is distinguishable at demonstrated parameters, but reaching efficient upconversion will require further modeling. The pulses in this paper have a picosecond duration and propagate through 2 mm of plasma, producing a high frequency wave of lower intensity than the pump. This is highly inefficient, but, given the significant pump energy, might be detectable in a sensitive experiment. There are challenges that must be addressed before this can be accomplished. While we have considered how waves in the system might seed the growth of plasma waves, our treatment has assumed that electrons are perfectly cold. Kinetic effects may significantly seed the growth of parasitic plasma waves, and it may be difficult to find a parameter regime in which these effects are balanced against significant growth. More detailed study will need to be done in higher growth regimes to examine the balance of these effects. Growth falls short of the requirements for application, where scaling the growth length results in a plasma channel length that is unachievable. This is primarily because a lower plasma density has been used than in the most promising upconversion schemes.12 The use of high densities and high growth has been limited by the development of plasma instabilities, which might be overcome in further simulation. Inhomogeneous plasma conditions, which disrupt the growth of plasma waves, could allow density, and thus coupling to be much higher. Even in a higher coupling regime, pulses will propagate for many more laser cycles. Resolving the shortest wavelengths and highest frequencies would be costly, but can be overcome by working with an envelope model, which would only need to resolve envelope length scales and nonlinear process timescales.
Working around these challenges to model and, eventually, experimentally produce a high energy, short wavelength pulse could be done through the application of a more sophisticated approach. Notably, the six-photon resonance here, which uses a single pump, has a much lower growth rate than those consisting of detuned pump waves. When the beating of the pumps is chosen precisely, the resonance achieves a smaller frequency multiplication, but with much higher growth.12 Parasitic Raman scattering might also be suppressed through varying the density of the plasma. Varying the density would be favorable, as the six-photon resonance considered here results from the coupling of electromagnetic waves, which are less sensitive to the plasma density than Langmuir waves; thus, the resonance is less sensitive than the Raman resonance. Though it is important to note that this will not remain the case if the beat wave between pumps becomes relevant for other six-photon processes. Higher order process requires a high intensity over a large propagation distance, so filamentation may become a problem if low diffraction and, thus, high power is needed. If these conditions are met, simulations performed using particle-in-cell methods could directly show the six-photon process taking into account kinetic effects. Six-photon scattering contains great promise, and iterated scattering, which might be automatically seeded, could produce geometric increases in frequency.
ACKNOWLEDGMENTS
The authors thank Vladimir Malkin for informative discussions. This research was supported by the NSF PHY-2206691, DOE DE-SC0021248, and DOE NNSA DE-NA0004167.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Alec Griffith: Investigation (equal); Methodology (equal); Software (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal). Kenan Qu: Investigation (equal); Supervision (equal); Writing – original draft (equal); Writing – review & editing (equal). Nathaniel J. Fisch: Conceptualization (equal); Funding acquisition (equal); Project administration (equal); Supervision (equal); Writing – review & editing (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.