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 3.73. 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.

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 ( 2 + 3 ) ω and ( 2 3 ) ω. 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, n 0 / n c = ω p e 2 ω 2, for laser frequency ω and plasma frequency ω p e 2 = 4 π e 2 n 0 m e 1, or normalized wave amplitude, e | A | / m c 2, for laser vector potential A. When both parameters are not 1, 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 ω p e < ω 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.

We review and reformulate how the coupling of electromagnetic waves may be achieved through relativistic corrections to electron dynamics in strong electromagnetic waves. The model in which the coupling occurs is that of an electron fluid, which is driven purely by electromagnetic forces, neglecting electron temperature and pressure along with any ion dynamics. Thus, the fluid is strictly described by Maxwell's equations and the coupled electron quantities of density and velocity. Working in the Coulomb gauge, the dynamics may be described by non-dimensionalized vector potential a = e A x / m e c 2, electrostatic potential ϕ = e Φ / m e c 2, and electron action s = S / m e c 2, resulting in three coupled equations as written by Malkin and Fisch,9 
(1)
(2)
(3)
These equations result from combing Maxwell's equations and the relativistic dynamics for electron particle motion, they neglect all kinetic effects under the assumption that the relativistic motion dominates the thermal motion.9 In the case we shall consider, we assume propagation purely in the z direction, with the vector potential a only having a transverse, x, component. In this case, the conservation of canonical momentum makes the electron perpendicular momentum p = mca, and the action may be related to the more easily understood quantity of parallel momentum through the relation p z = m c 2 z s. To complete the correspondence between Eqs. (1)–(3) and cold electron hydrodynamics, electron density can be straightforwardly related to the electrostatic potential through Gauss's law, ω p e 2 + c 2 z 2 ϕ = ω p e 2 n e / n 0.

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 γ = 1 + a 2 + ( c z s ) 2. 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 six-photon scattering process may be characterized by a simple set of resonance conditions and the asymptotic behavior of the coupling. Perpendicularly polarized waves have dispersion relation, written as ω j 2 = k j 2 + n, which is normalized for k 0 = 1 and c = 1. These waves may couple to one another through the relativistic nonlinearity. The fifth order nonlinearity supports several different six-photon scattering processes, but the simplest case consists of one where four input photons at frequency ω0 and wavenumber k0 are turned into two photons, one at ω 1 , k 1 and one at ω 2 , k 2, meeting the following resonance conditions:
(4)
These conditions may be solved resulting in the following resonant frequencies:
(5)
The growth of the envelopes of these modes in the system can be represented by a set of dynamical equations in amplitudes, aj, of
(6)
for a dimensionless factor V depending on the particularities of the resonance, which is 1 for the range of parameters considered in this paper. This resonance is unique in that it can occur for co-linear wavevectors in a plasma without renormalization of the laser frequency, which is not possible for lower order relativistic scattering.9 Furthermore, in the limit, this produces a unique non-integer resonance, in contrast to typical integer harmonics. The growth rate of the six-photon coupling, Γ, normalized to the laser frequency, ω0, scales as follows:
(7)
The resulting gain length, L, may be estimated by the expression L = V π ( ω 0 ω p e ) 6 λ 0 a 0 4. The form of the leading coefficient and the corrections in higher orders of plasma density may be seen in the exposition by Malkin and Fisch.12 From Eq. (7), it can be seen that the growth of any six-photon scattering is highly sensitive to both wave amplitude and density.

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.

The numerical evolution of the relativistic scattering process is chosen to maximize the observability of the six-photon scattering process. To capture a distinguishable signal, the pseudospectral package Dedalus32 is used in order to minimize numerical noise in the Fourier spectrum, in contrast to high noise particle in cell results. The nonlinear pulse propagation is non-dimensionalized and transformed into the co-propagating frame η = t , ξ = z t,
(8)
(9)
(10)
where c = 1, and density is normalized such that n = ω p e 2 / k 0 2. Given that only the derivatives of ϕ and s couple to the other equation, for convenience, we reduce the equations with the substitutions z ϕ = e and c z s = l. To ensure numerical stability, a viscosity term is added to curtail the rapid nonphysical growth of high frequency modes driven by the complete absence of physical dissipation.

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 k 0 = 1 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 k 1 = 4 k 0 k 2, we demonstrate that Eq. (5) describes a growing mode, which is only seen when it is satisfied.

FIG. 1.

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.

FIG. 1.

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.

Close modal

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 4 k 0 k 2. 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 k 0 / k p e = 7.5. 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 a 0 = 0.5, 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 8 × 10 5 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 ω 2 = 2 ω p e, and the maximum frequency bounded by the coupling of the third harmonic of the low frequency seed to the stokes shifted pump wave, where 3 ω 2 = ω 0 ω p e. 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 1.6 × 10 7 ω 0 and a fit to the growth in simulations gives a growth rate of 1.8 × 10 7 ω 0. 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 k = 2 k 0 k 2, which experiences some scattering, but is not resonant in the colinear geometry without further tuning of the pump and seed wave.9 

FIG. 2.

The evolution of the spectrum of a pump ( k 0 = 1 , a 0 = 0.5) and seed [ k 2 ( 2 3 ) k 0 , I 2 = I 0 / 200]. 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 4 k 0 k 2 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 k p e / 4.

FIG. 2.

The evolution of the spectrum of a pump ( k 0 = 1 , a 0 = 0.5) and seed [ k 2 ( 2 3 ) k 0 , I 2 = I 0 / 200]. 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 4 k 0 k 2 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 k p e / 4.

Close modal

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 4 k 0 k 2 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.

FIG. 3.

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 k 2 = 4 k 0 [ ( 2 1 + n + 3 ) 2 n ] 1 / 2 k 0 = 0.252 k 0 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.

FIG. 3.

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 k 2 = 4 k 0 [ ( 2 1 + n + 3 ) 2 n ] 1 / 2 k 0 = 0.252 k 0 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.

Close modal

Some initial small amount of energy can be transferred into the 4 k 0 k 2 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 10 5 below the initial pump energy. To reach the nonlinear growth regime, the resonance has the condition δ ω 2 / ω 2 1 2 n 2 ω 0 4 a 0 4. For the parameters in the simulations presented, this corresponds to a maximal detuning of 10 5. 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 ω = 4 ω 0 3 k 0. 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.

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 10 5 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.

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.

The authors have no conflicts to disclose.

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).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
M. R.
Edwards
and
J. M.
Mikhailova
, “
The x-ray emission effectiveness of plasma mirrors: Reexamining power-law scaling for relativistic high-order harmonic generation
,”
Sci. Rep.
10
,
5154
(
2020
).
2.
P.
Emma
,
R.
Akre
,
J.
Arthur
,
R.
Bionta
,
C.
Bostedt
,
J.
Bozek
,
A.
Brachmann
,
P.
Bucksbaum
,
R.
Coffee
,
F.-J.
Decker
,
Y.
Ding
,
D.
Dowell
,
S.
Edstrom
,
A.
Fisher
,
J.
Frisch
,
S.
Gilevich
,
J.
Hastings
,
G.
Hays
,
P.
Hering
,
Z.
Huang
,
R.
Iverson
,
H.
Loos
,
M.
Messerschmidt
,
A.
Miahnahri
,
S.
Moeller
,
H.-D.
Nuhn
,
G.
Pile
,
D.
Ratner
,
J.
Rzepiela
,
D.
Schultz
,
T.
Smith
,
P.
Stefan
,
H.
Tompkins
,
J.
Turner
,
J.
Welch
,
W.
White
,
J.
Wu
,
G.
Yocky
, and
J.
Galayda
, “
First lasing and operation of an ångstrom-wavelength free-electron laser
,”
Nat. Photonics
4
,
641
647
(
2010
).
3.
C.
Pellegrini
,
A.
Marinelli
, and
S.
Reiche
, “
The physics of x-ray free-electron lasers
,”
Rev. Mod. Phys.
88
,
015006
(
2016
).
4.
N. D.
Powers
,
I.
Ghebregziabher
,
G.
Golovin
,
C.
Liu
,
S.
Chen
,
S.
Banerjee
,
J.
Zhang
, and
D. P.
Umstadter
, “
Quasi-monoenergetic and tunable X-rays from a laser-driven Compton light source
,”
Nat. Photonics
8
,
28
31
(
2014
).
5.
O.
Hort
,
M.
Albrecht
,
V. E.
Nefedova
,
O.
Finke
,
D. D.
Mai
,
S.
Reyné
,
F.
Giambruno
,
F.
Frassetto
,
L.
Poletto
,
J.
Andreasson
,
J.
Gautier
,
S.
Sebban
, and
J.
Nejdl
, “
High-flux source of coherent XUV pulses for user applications
,”
Opt. Express
27
,
8871
8883
(
2019
).
6.
U.
Chaulagain
,
M.
Lamač
,
M.
Raclavský
,
K. P.
Khakurel
,
K. H.
Rao
,
K.
Ta-Phuoc
,
S. V.
Bulanov
, and
J.
Nejdl
, “
ELI Gammatron beamline: A dawn of ultrafast hard x-ray science
,”
Photonics
9
,
853
(
2022
).
7.
C. A.
Haynam
,
P. J.
Wegner
,
J. M.
Auerbach
,
M. W.
Bowers
,
S. N.
Dixit
,
G. V.
Erbert
,
G. M.
Heestand
,
M. A.
Henesian
,
M. R.
Hermann
,
K. S.
Jancaitis
,
K. R.
Manes
,
C. D.
Marshall
,
N. C.
Mehta
,
J.
Menapace
,
E.
Moses
,
J. R.
Murray
,
M. C.
Nostrand
,
C. D.
Orth
,
R.
Patterson
,
R. A.
Sacks
,
M. J.
Shaw
,
M.
Spaeth
,
S. B.
Sutton
,
W. H.
Williams
,
C. C.
Widmayer
,
R. K.
White
,
S. T.
Yang
, and
B. M. V.
Wonterghem
, “
National Ignition Facility laser performance status
,”
Appl. Opt.
46
,
3276
3303
(
2007
).
8.
V. M.
Malkin
and
N. J.
Fisch
, “
Towards megajoule x-ray lasers via relativistic four-photon cascade in plasma
,”
Phys. Rev. E
101
,
023211
(
2020
).
9.
V. M.
Malkin
and
N. J.
Fisch
, “
Resonant four-photon scattering of collinear laser pulses in plasma
,”
Phys. Rev. E
102
,
063207
(
2020
).
10.
V. M.
Malkin
and
N. J.
Fisch
, “
Super-resonant four-photon collinear laser frequency multiplication in plasma
,”
Phys. Rev. E
105
,
045207
(
2022
).
11.
A.
Griffith
,
K.
Qu
, and
N. J.
Fisch
, “
Modulation-slippage trade-off in resonant four-wave upconversion
,”
Phys. Plasmas
28
,
052112
(
2021
).
12.
V. M.
Malkin
and
N. J.
Fisch
, “
Six-photon resonant scattering of collinear laser pulses in plasma
,”
Phys. Rev. E
108
,
045208
(
2023
).
13.
V. M.
Malkin
,
G.
Shvets
, and
N. J.
Fisch
, “
Fast compression of laser beams to highly overcritical powers
,”
Phys. Rev. Lett.
82
,
4448
4451
(
1999
).
14.
C.-H.
Pai
,
M.-W.
Lin
,
L.-C.
Ha
,
S.-T.
Huang
,
Y.-C.
Tsou
,
H.-H.
Chu
,
J.-Y.
Lin
,
J.
Wang
, and
S.-Y.
Chen
, “
Backward Raman amplification in a plasma waveguide
,”
Phys. Rev. Lett.
101
,
065005
(
2008
).
15.
Y.
Ping
,
R. K.
Kirkwood
,
T.-L.
Wang
,
D. S.
Clark
,
S. C.
Wilks
,
N.
Meezan
,
R. L.
Berger
,
J.
Wurtele
,
N. J.
Fisch
,
V. M.
Malkin
,
E. J.
Valeo
,
S. F.
Martins
, and
C.
Joshi
, “
Development of a nanosecond-laser-pumped Raman amplifier for short laser pulses in plasma
,”
Phys. Plasmas
16
,
123113
(
2009
).
16.
R. M. G. M.
Trines
,
F.
Fiúza
,
R.
Bingham
,
R. A.
Fonseca
,
L. O.
Silva
,
R. A.
Cairns
, and
P. A.
Norreys
, “
Simulations of efficient Raman amplification into the multipetawatt regime
,”
Nat. Phys.
7
,
87
92
(
2011
).
17.
M. R.
Edwards
,
K.
Qu
,
J. M.
Mikhailova
, and
N. J.
Fisch
, “
Beam cleaning of an incoherent laser via plasma Raman amplification
,”
Phys. Plasmas
24
,
103110
(
2017
).
18.
K.
Qu
,
I.
Barth
, and
N. J.
Fisch
, “
Plasma wave seed for Raman amplifiers
,”
Phys. Rev. Lett.
118
,
164801
(
2017
).
19.
R.
Miriam Cheriyan
,
N.
Varghese
,
R. S.
Sooraj
,
K. H.
Rao
, and
N.
Smijesh
, “
A comprehensive review on amplification of laser pulses via stimulated Raman scattering and stimulated Brillouin scattering in plasmas
,”
Plasma
5
,
499
539
(
2022
).
20.
A. A.
Andreev
,
C.
Riconda
,
V. T.
Tikhonchuk
, and
S.
Weber
, “
Short light pulse amplification and compression by stimulated Brillouin scattering in plasmas in the strong coupling regime
,”
Phys. Plasmas
13
,
053110
(
2006
).
21.
S.
Weber
,
C.
Riconda
,
L.
Lancia
,
J.-R.
Marquès
,
G. A.
Mourou
, and
J.
Fuchs
, “
Amplification of ultrashort laser pulses by Brillouin backscattering in plasmas
,”
Phys. Rev. Lett.
111
,
055004
(
2013
).
22.
M. R.
Edwards
,
J. M.
Mikhailova
, and
N. J.
Fisch
, “
X-ray amplification by stimulated Brillouin scattering
,”
Phys. Rev. E
96
,
023209
(
2017
).
23.
M. R.
Edwards
,
Y.
Shi
,
J. M.
Mikhailova
, and
N. J.
Fisch
, “
Laser amplification in strongly magnetized plasma
,”
Phys. Rev. Lett.
123
,
025001
(
2019
).
24.
R. K.
Kirkwood
,
D. P.
Turnbull
,
T.
Chapman
,
S. C.
Wilks
,
M. D.
Rosen
,
R. A.
London
,
L. A.
Pickworth
,
A.
Colaitis
,
W. H.
Dunlop
,
P.
Poole
,
J. D.
Moody
,
D. J.
Strozzi
,
P. A.
Michel
,
L.
Divol
,
O. L.
Landen
,
B. J.
MacGowan
,
B. M.
Van Wonterghem
,
K. B.
Fournier
, and
B. E.
Blue
, “
A plasma amplifier to combine multiple beams at NIF
,”
Phys. Plasmas
25
,
056701
(
2018
).
25.
R. K.
Kirkwood
,
D. P.
Turnbull
,
T.
Chapman
,
S. C.
Wilks
,
M. D.
Rosen
,
R. A.
London
,
L. A.
Pickworth
,
W. H.
Dunlop
,
J. D.
Moody
,
D. J.
Strozzi
,
P. A.
Michel
,
L.
Divol
,
O. L.
Landen
,
B. J.
MacGowan
,
B. M.
Van Wonterghem
,
K. B.
Fournier
, and
B. E.
Blue
, “
Plasma-based beam combiner for very high fluence and energy
,”
Nat. Phys.
14
,
80
84
(
2018
).
26.
R. K.
Kirkwood
,
P. L.
Poole
,
D. H.
Kalantar
,
T. D.
Chapman
,
S. C.
Wilks
,
M. R.
Edwards
,
D. P.
Turnbull
,
P.
Michel
,
L.
Divol
,
N. J.
Fisch
,
P.
Norreys
,
W.
Rozmus
,
J.
Bude
,
B. E.
Blue
,
K. B.
Fournier
,
B. M.
Van Wonterghem
, and
A.
MacKinnon
, “
Production of high fluence laser beams using ion wave plasma optics
,”
Appl. Phys. Lett.
120
,
200501
(
2022
).
27.
P.
Michel
,
L.
Divol
,
D.
Turnbull
, and
J. D.
Moody
, “
Dynamic control of the polarization of intense laser beams via optical wave mixing in plasmas
,”
Phys. Rev. Lett.
113
,
205001
(
2014
).
28.
G.
Lehmann
and
K. H.
Spatschek
, “
Plasma-based polarizer and waveplate at large laser intensity
,”
Phys. Rev. E
97
,
063201
(
2018
).
29.
M. R.
Edwards
,
V. R.
Munirov
,
A.
Singh
,
N. M.
Fasano
,
E.
Kur
,
N.
Lemos
,
J. M.
Mikhailova
,
J. S.
Wurtele
, and
P.
Michel
, “
Holographic plasma lenses
,”
Phys. Rev. Lett.
128
,
065003
(
2022
).
30.
M. R.
Edwards
and
P.
Michel
, “
Plasma transmission gratings for compression of high-intensity laser pulses
,”
Phys. Rev. Appl.
18
,
024026
(
2022
).
31.
J. M.
Rax
and
N. J.
Fisch
, “
Third-harmonic generation with ultrahigh-intensity laser pulses
,”
Phys. Rev. Lett.
69
,
772
775
(
1992
).
32.
K. J.
Burns
,
G. M.
Vasil
,
J. S.
Oishi
,
D.
Lecoanet
, and
B. P.
Brown
, “
Dedalus: A flexible framework for numerical simulations with spectral methods
,”
Phys. Rev. Res.
2
,
023068
(
2020
).