Nonlinear Laser-plasma Interactions Research Articles

First Simulations for the EuAPS Betatron Radiation Source: A Dedicated Radiation Calculation Code

X-ray production through betatron radiation emission from electron bunches is a valuable resource for several research fields. The EuAPS (EuPRAXIA Advanced Photon Sources) project, within the framework of EuPRAXIA, aims to provide 1–10 keV photons (X-rays), developing a compact plasma-based system designed to exploit self-injection processes that occur in the highly nonlinear laser-plasma interaction (LWFA) to drive electron betatron oscillations. Since the emitted radiation spectrum, intensity, angular divergence, and possible coherence strongly depend on the properties of the self-injected beam, accurate preliminary simulations of the process are necessary to evaluate the optimal diagnostic device specifications and to provide an initial estimate of the source’s performance. A dedicated tool for these tasks has been developed; electron trajectories from particle-in-cell (PIC) simulations are currently undergoing numerical analysis through the calculation of retarded fields and spectra for various plasma and laser parameter combinations. The implemented forward approach evaluation of the fields could allow for the integration of the presented scheme into already existing PIC codes. The spectrum calculation is thus performed in detector time, giving a linear complex exponential phase; this feature allows for a semi-analitical Fourier transform evaluation. The code structure and some trajectories analysis results are presented.

Femtosecond laser-induced sub-wavelength plasma inside dielectrics. II. Second-harmonic generation

Second-harmonic emission at a frequency that is twice the laser frequency is an important diagnostic for nonlinear laser–plasma interaction. It is forbidden for centrosymmetric materials such as the bulk of sapphire. The symmetry, however, can be broken by dielectric discontinuities as a result of plasma generation inside a solid dielectric. In the present work, we explore the basic characteristics of experimentally observed second-harmonic emission during focusing a femtosecond Bessel beam inside sapphire. We employ three-dimensional particle-in-cell simulations and the Helmholtz wave equation for theoretical investigations. We analyze how the efficiency of second-harmonic generation and its polarization depend on the plasma parameters. We find that the second-harmonic is generated either due to the coalescence of two-surface electromagnetic waves or nonlinear interaction between the transverse electromagnetic wave and the longitudinal electron plasma wave driven by linear mode conversion. Experimental results agree with the theoretical predictions and confirm the existence of over-critical plasma inside the sapphire that is essential for the resonance of plasma waves or excitation of surface plasmons.

Inverse ray tracing on icosahedral tetrahedron grids for non-linear laser plasma interaction coupled to 3D radiation hydrodynamics

A novel approach to efficiently model 3-D laser plasma interactions at fluid scales is presented. This method, implemented in the IFRIIT propagation code developed at CELIA, relies on inverse ray tracing to compute laser fields at arbitrary locations in a plasma. This enables to describe the fields at high order in space compared to standard forward ray tracing approaches. In addition, inverse ray tracing enables the use of etalon integral methods to reconstruct caustic fields and greatly speeds up calculations of cross-beam energy transfer by decoupling the ray amplitude and ray phase calculations. A comparison of the inverse and forward methods for 3-D calculations of fields in presence or not of cross-beam energy transfer illustrates the significant advantages of the inverse method. Conversely, while the inverse method is well suited to most spherical plasma profiles, it currently cannot treat concave profiles or target holders. The coupling of IFRIIT with the 3-D ASTER radiative hydrodynamics code developed at the Laboratory for Laser Energetics is then presented. ASTER and IFRIIT resolve their respective equations on separate grids which communicate through interpolation. As such, IFRIIT uses a dedicated laser grid adapted to the computations at play, which also allows to use different parallelization methods for both codes: block decomposition for the hydrodynamics versus domain duplication for the laser. Applications to direct-drive implosions for inertial confinement fusion are presented, for which a geodesic icosahedron grid is implemented in IFRIIT. The performances of the ASTER/IFRIIT coupling are demonstrated by conducting simulations of cryogenic implosions performed on the OMEGA laser system, in presence of various sources of 3-D effects; laser port geometry, cross-beam energy transfer, beam imbalance and target mis-alignment. Comparison with neutron data, measured through bang-time, for a cryogenic implosion experiment shows an excellent agreement for the laser-plasma coupling.

High yield polar direct drive fusion neutron sources at the National Ignition Facility

Polar direct drive neutron source experiments were performed at the National Ignition Facility showing substantial improvement in total neutron yield and efficiency of conversion of laser energy to fusion output. Plastic capsules 3–4 mm in diameter were filled with 1.5 mg/cc of deuterium–tritium (DT) fuel and imploded with laser beam pointing and defocus designed to compensate for polar asymmetry introduced by the facility beam entrance angles. Radiation-hydrodynamics simulations were employed to optimize the multi-dimensional laser and target parameter space, within facility and target fabrication constraints. Ensembles of 1D simulations tuned to match the outputs of early shots in the series were used to design subsequent shots in the series. This allowed the later shots to be designed based on empirically motivated sensitivities to laser and target input parameters, while eliminating the need to explicitly model phenomena such as hydrodynamic instabilities and nonlinear laser–plasma interactions. One experiment with a 3.0 mm diameter CH capsule produced 13.6 kJ (4.81 × 1015 DT neutrons) from a laser input below the NIF optics damage threshold at 585 kJ, 328 TW. Two experiments with 4.0 mm capsules produced 31.3 and 33.6 kJ of fusion output (1.11 × 1016 and 1.19 × 1016 DT neutrons) with 1.10 MJ, 390 TW and 1.26 MJ, 425 TW of laser input, respectively.

Spectral structures of backward stimulated Brillouin scattering driven by a picosecond laser

Laser plasma interaction (LPI) is an important content in laser plasma related research, and it is one of the key issues related to the success or failure of inertial confinement fusion ignition, and has received extensive attention. In order to suppress the relevant LPI process as much as possible, the major laboratories around the world have developed a variety of beam smoothing methods through decades of research. However, the current understanding and suppression of LPI are still far from enough, and further in-depth studies are still needed. Generally, the research of LPI is based on nanosecond laser driving, and focuses mainly on the effects of the related LPI process caused by nanosecond lasers. However, the LPI processes, such as stimulated Brillouin scattering (SBS), stimulated Raman scattering (SRS), etc., occur and develop on a time scale of picoseconds.The comprehensive effect can be studied only on a longer scale of nanosecond. For highly nonlinear LPI processes, the comprehensive effect may be difficult to reflect the real physical laws. The emergence of the picosecond laser has made it possible to study the LPI process in more detail and on a more appropriate time scale. The present research tries to gain an understanding of LPI from a more refined perspective. The experimental research of picosecond laser driving LPI is carried out on the Shenguang-Ⅱ upgrade and picosecond laser facilities. First, a nanosecond laser is used to irradiate a target to generate a large-scale plasma, and a few nanoseconds later, the picosecond laser is injected as an interaction beam to drive the LPI scattering such as SBS and SRS. The spectral signal of backscatter light is measured experimentally by using the method of diffuse reflector. From the research results it is found that the backward signals of the band near the laser wavelength contain, in addition to the true backward SBS component, a large number of interference signals introduced by picosecond laser and nanosecond laser. The interference signal introduced by nanosecond laser can be eliminated by using specific measures, but the interference signal introduced by picosecond laser cannot be eliminated experimentally, which will affect the estimation of the true share of the backward SBS. The comprehensive results show that under different experimental conditions, the backward scatter energy of SBS may be less than half that of the total recorded signals. This result is helpful in further understanding and re-recognizing previous relevant experimental data.

Generation of Terawatt-Scale Vortex Pulses Based on Optical Parametric Chirped-Pulse Amplification

In this letter, we report the demonstration of terawatt (TW)-scale femtosecond vortex laser pulses based on noncollinear optical parametric chirped-pulse amplification (OPCPA) centered at 800 nm. To the best of our knowledge, this is the first experimental study of amplification for OPCPA-based vortex pulses to achieve a powerful ultrafast (TW-scale) vortex pulse output. In this experiment, a 0.5-mJ broadband chirped vortex pulse was amplified to 47 mJ with a 15-mm-long LiB3O5 crystal. The double-pass grating pair compressor with an efficiency of 65% maintained the beam vortex, and a pulse of width 30.1 fs was obtained, corresponding to a peak power of 1.02 TW. This laser can be used in studies on nonlinear optics and relativistic laser-plasma interactions, compact plasma-based accelerators, and light sources.

Characterization of suprathermal electrons inside a laser accelerated plasma via highly-resolved K\u237a-emission

Suprathermal electrons are routinely generated in high-intensity laser produced plasmas via instabilities driven by non-linear laser-plasma interaction. Their accurate characterization is crucial for the performance of inertial confinement fusion as well as for performing experiments in laboratory astrophysics and in general high-energy-density physics. Here, we present studies of non-thermal atomic states excited by suprathermal electrons in kJ-ns-laser produced plasmas. Highly spatially and spectrally resolved X-ray emission from the laser-deflected part of the warm dense Cu foil visualized the hot electrons. A multi-scale two-dimensional hydrodynamic simulation including non-linear laser-plasma interactions and hot electron propagation has provided an input for ab initio non-thermal atomic simulations. The analysis revealed a significant delay between the maximum of laser pulse and presence of suprathermal electrons. Agreement between spectroscopic signatures and simulations demonstrates that combination of advanced high-resolution X-ray spectroscopy and non-thermal atomic physics offers a promising method to characterize suprathermal electrons inside the solid density matter.

Real and complex valued geometrical optics inverse ray-tracing for inline field calculations

A 3-D ray based model for computing laser fields in dissipative and amplifying media is presented. The eikonal equation is solved using inverse ray-tracing on a dedicated nonstructured 3-D mesh. Inverse ray-tracing opens the possibility of using Complex Geometrical Optics (CGO), for which we propose a propagation formalism in a finite element mesh. Divergent fields at caustics are corrected using an etalon integral method for fold-type caustics. This method is successfully applied in dissipative media by modifying the ray-ordering and root selection rules, thereby allowing one to reconstruct the field in the entire caustic region. In addition, we demonstrate how caustics in the CGO framework can disappear entirely for sufficiently dissipative media, making the complex ray approach valid in the entire medium. CGO is shown to offer a more precise modeling of laser refraction and absorption in a dissipative medium when compared to Geometrical Optics (GO). In the framework of Inertial Confinement Fusion (ICF), this occurs mostly at intermediate temperatures or at high temperatures close to the critical density. Additionally, GO is invalid at low temperatures if an approximated expression of the permittivity is used. The inverse ray-tracing algorithm for GO and CGO is implemented in the IFRIIT code, in the framework of a dielectric permittivity described in 3-D using a piecewise linear approximation in tetrahedrons. Fields computed using GO and CGO are compared to results from the electromagnetic wave solver Laser Plasma Simulation Environment. Excellent agreement is obtained in 1-D linear and nonlinear permittivity profiles. Good agreement is also obtained for ICF-like Gaussian density profiles in 2-D. Finally, we demonstrate how the model reproduces Gaussian beam diffraction using CGO. The IFRIIT code will be interfaced inline to 3-D radiative hydrodynamic codes to describe the nonlinear laser plasma interaction in ICF and high-energy-density plasmas.

Real and virtual propagation dynamics of angular accelerating white light beams.

Accelerating waves have received significant attention of late, first in the optical domain and later in the form of electron matter waves, and have found numerous applications in non-linear optics, material processing, microscopy, particle manipulation and laser plasma interactions. Here we create angular accelerating light beams with a potentially unlimited acceleration rate. By employing wavelength independent digital holograms for the creation and propagation of white light beams, we are able to study the resulting propagation in real and virtual space. We find that dephasing occurs for real propagation and that this can be compensated for in a virtual propagation scheme when single plane dynamics are important. Our work offers new insights into the propagation dynamics of such beams and provides a versatile tool for further investigations into propagating structured light fields.

Hybrid-PIC modeling of laser-plasma interactions and hot electron generation in gold hohlraum walls

The walls of the hohlraum used in experiments at the national ignition facility are heated by laser beams with intensities ∼1015 W/cm2, a wavelength of ∼1/3 μm, and pulse lengths on the order of a ns, with collisional absorption believed to be the primary heating mechanism. X-rays generated by the hot ablated plasma at the gold walls are then used to implode a target in the hohlraum interior. In addition to the collisional absorption of laser energy at the walls, non-linear laser-plasma interactions (LPI), such as stimulated Raman scattering and two plasmon decay, are believed to generate a population of supra-thermal electrons which, if present in the hohlraum, can have a deleterious effect on target implosion. We describe results of hohlraum modeling using a hybrid particle-in-cell code. To enable this work, new particle-based algorithms for a multiple-ion magneto-hydrodynamic (MHD) treatment, and a particle-based ray-tracing model were developed. The use of such hybrid methods relaxes the requirement to resolve the laser wavelength, and allows for relatively large-scale hohlraum simulations with a reasonable number of cells. But the non-linear effects which are believed to be the cause of hot electron generation can only be captured by fully kinetic simulations with good resolution of the laser wavelength. For this reason, we employ a two-tiered approach to hohlraum modeling. Large-scale simulations of the collisional absorption process can be conducted using the fast quasi-neutral MHD algorithm with fluid particle species. From these simulations, we can observe the time evolution of the hohlraum walls and characterize the density and temperature profiles. From these results, we can transition to smaller-scale highly resolved simulations using traditional kinetic particle-in-cell methods, from which we can fully model all of the non-linear laser-plasma interactions, as well as assess the details of the electron distribution function. We find that vacuum hohlraums should be stable to both two plasmon decay and stimulated Raman scattering instabilities for intensities ≤1015 W/cm2. In gas-filled hohlraums, shocks may be induced in the blowoff gold plasma, which leads to more complex density and temperatures profiles. The resulting effect on LPI stability depends strongly on the details of the profile, and it is possible for the gas-filled hohlraum to become unstable to two plasmon decay at 1015 W/cm2 if the quarter-critical surface reaches temperatures exceeding 1 keV.

Experimental observation of parametric instabilities at laser intensities relevant for shock ignition

We report measurements of parametric instabilities and hot electron generation in a laser intensity regime up to 6 × 1015 W/cm2, typical of the shock ignition approach to inertial fusion. Experiments performed at the PALS laboratory in Prague show that the incident laser energy losses are dominated by Stimulated Brillouin Scattering (SBS) rather than by Stimulated Raman Scattering (SRS) or Two-Plasmon Decay (TPD). Results are compared to hydrodynamics simulations using a code that includes self-consistent calculations of non-linear laser plasma interactions and accounts for the laser intensity statistics contained in the beam speckles. Good agreement is found for the backscattered SRS light, and for temperature and flux of hot electrons. The effect of high-intensity speckles on backscattered SRS is also underlined numerically and experimentally.

Influence of laser induced hot electrons on the threshold for shock ignition of fusion reactions

The effects of Hot Electrons (HEs) generated by the nonlinear Laser-Plasma Interaction (LPI) on the dynamics of Shock Ignition Inertial Confinement Fusion targets are investigated. The coupling between the laser beam, plasma dynamics and hot electron generation and propagation is described with a radiative hydrodynamics code using an inline model based on Paraxial Complex Geometrical Optics [Colaïtis et al., Phys. Rev. E 92, 041101 (2015)]. Two targets are considered: the pure-DT HiPER target and a CH-DT design with baseline spike powers of the order of 200–300 TW. In both cases, accounting for the LPI-generated HEs leads to non-igniting targets when using the baseline spike powers. While HEs are found to increase the ignitor shock pressure, they also preheat the bulk of the imploding shell, notably causing its expansion and contamination of the hotspot with the dense shell material before the time of shock convergence. The associated increase in hotspot mass (i) increases the ignitor shock pressure required to ignite the fusion reactions and (ii) significantly increases the power losses through Bremsstrahlung X-ray radiation, thus rapidly cooling the hotspot. These effects are less prominent for the CH-DT target where the plastic ablator shields the lower energy LPI-HE spectrum. Simulations using higher laser spike powers of 500 TW suggest that the CH-DT capsule marginally ignites, with an ignition window width significantly smaller than without LPI-HEs, and with three quarters of the baseline target yield. The latter effect arises from the relation between the shock launching time and the shell areal density, which becomes relevant in presence of a LPI-HE preheating.

Towards modeling of nonlinear laser-plasma interactions with hydrocodes: the thick-ray approach.

This paper deals with the computation of laser beam intensity in large-scale radiative hydrocodes applied to the modeling of nonlinear laser-plasma interactions (LPIs) in inertial confinement fusion (ICF). The paraxial complex geometrical optics (PCGO) is adapted for light waves in an inhomogeneous medium and modified to include the inverse bremsstrahlung absorption and the ponderomotive force. This thick-ray model is compared to the standard ray-tracing (RT) approach, both in the chic code. The PCGO model leads to different power deposition patterns and better diffraction modeling compared to standard RT codes. The intensity-reconstruction technique used in RT codes to model nonlinear LPI leads to artificial filamentation and fails to reproduce realistic ponderomotive self-focusing distances, intensity amplifications, and density channel depletions, whereas PCGO succeeds. Bundles of Gaussian thick rays can be used to model realistic non-Gaussian ICF beams. The PCGO approach is expected to improve the accuracy of ICF simulations and serve as a basis to implement diverse LPI effects in large-scale hydrocodes.

Higher-Order Nonlocal Effects of a Relativistic Ponderomotive Force in High-Intensity Laser Fields

We have developed a new formula for a relativistic ponderomotive force of transversely localized laser fields based on the noncanonical Lie perturbation method by finding proper coordinates and gauges in the variational principle. The formula involves new terms represented by second and third spatial derivatives of the field amplitude, so that the ponderomotive force depends not only on the local field gradient, but also on the curvature and its variation. The formula is then applicable to a regime in which the conventional formula is hardly applied such that nonlocal and/or global extent of the field profile becomes important. The result can provide a theoretical basis for describing nonlinear laser-plasma interaction including such nonlocal effects, which is examined via particle-in-cell simulation of laser propagation in a plasma with a super Gaussian transverse field profile.

Intense terahertz emission from relativistic circularly polarized laser pulses interaction with overdense plasmas

During the interaction of a relativistic circularly polarized laser pulse with an overdense plasma target, the longitudinal motion of bunches of electrons under the action of light pressure and electrostatic restore force can emit intense terahertz (THz) pulses. This mechanism allows high pump laser intensity and large electron number participating in the emission. Two-dimensional particle-in-cell simulations are carried out to investigate the THz emission. The results suggest that such a source can produce remarkably intense THz pulses with energy of several mJ/sr and power of tens of gigawatts, which could find applications in nonlinear studies and relativistic laser-plasma interaction diagnostics.

Attosecond lighthouses from plasma mirrors

The nonlinear interaction of an intense femtosecond laser pulse with matter can lead to the emission of a train of sub-laser-cycle--attosecond--bursts of short-wavelength radiation1, 2. Much effort has been devoted to producing isolated attosecond pulses, as these are better suited to real-time imaging of fundamental electronic processes3, 4, 5, 6. Successful methods developed so far rely on confining the nonlinear interaction to a single sub-cycle event7, 8, 9. Here, we demonstrate for the first time a simpler and more universal approach to this problem10, applied to nonlinear laser-plasma interactions. By rotating the instantaneous wavefront direction of an intense few-cycle laser field11, 12 as it interacts with a solid-density plasma, we separate the nonlinearly generated attosecond pulse train into multiple beams of isolated attosecond pulses propagating in different and controlled directions away from the plasma surface. This unique method produces a manifold of isolated attosecond pulses, ideally synchronized for initiating and probing ultrafast electron motion in matter.

Experimental studies of the high and low frequency electromagnetic radiation produced from nonlinear laser-plasma interactions

We present the results of experiments on the generation of both high and low frequency electromagnetic radiation from the nonlinear interaction of intense laser beams with underdense plasmas. High frequency radiation is generated as a result of the reflection of a near-infrared laser beam by breaking plasma waves. In the proof-of-principle experiment, a 2 TW, 76 fs Ti:sapphire laser pulse was focused to generate wake waves, and a ~0.1 TW laser pulse was focused to collide with the electron density modulations. Frequencies of the reflected light were 56–110 times higher than that of the initial laser. The estimated number of photons per solid angle was 3×107 photons/sr. Low frequency radiation is emitted by the plasma as well. We found that the polarization of the radiation was similar to that of the driver laser. These properties are in accordance with the generation mechanism by relativistic solitons.

Nonlinear laser focusing using a conical guide and generation of energetic ions

Using conventional methods, a laser pulse can be focused down to around 6-8 microm, but further reduction of the spot size has proven to be difficult. Here it is shown by particle-in-cell simulation that with a hollow cone an intense laser pulse can be reduced to a tiny, highly localized, spot of around 1 microm radius, accompanied by much enhanced light intensity. The pulse shaping and focusing effect is due to a nonlinear laser-plasma interaction on the inner surface of the cone. When a thin foil is attached to the tip of the cone, the cone-focused light pulse compresses and accelerates the ions in its path and can punch through the thin target, creating highly localized energetic ion bunches of high density.

Measurement of preheat due to fast electrons in laser implosions of cryogenic deuterium targets

Preheat due to fast electrons generated by nonlinear laser-plasma interactions can reduce the gain in laser-imploded fusion targets. The first measurements of electron preheat in directly driven, imploding cryogenic-deuterium targets are reported. The preheat level is derived directly from the measured hard-x-ray spectrum. The hard-x-ray detectors were calibrated in situ with a parylene-coated molybdenum solid-sphere target where both the hard x rays and the Mo Kα line were measured. Using electron and radiation transport calculations to relate the measurements to preheat, the fraction of the incident laser energy that preheats the cryogenic deuterium fuel is found to be less than 0.1%, suggesting that preheat will have a negligible impact on target performance. The results are encouraging for the success of planned high-gain direct-drive-ignition experiments on the National Ignition Facility [W. J. Hogan, E. I. Moses, B. E. Warner, M. S. Sorem, and J. M. Soures, Nucl. Fusion 41, 567 (2001)].

Studies of stimulated Raman backscattering and associated nonlinear laser–plasma interactions

Coupled theoretical and computational work is presented aimed at understanding and modelling stimulated Raman backscattering (SRBS) relevant to laser–plasma interactions in large-scale, homogeneous plasmas. With the aid of a new code for simulating and studying the nonlinear coupling in space–time of a large number of modes, and an Eulerian Vlasov–Maxwell code for studying the evolution of large amplitude electron plasma waves, we report results and their interpretations to elucidate the following five observed, nonlinear phenomena associated with SRBS: coupling of SRBS to Langmuir decay instabilities (LDIs); effect of ion-acoustic damping on SRBS; cascading of LDI; stimulated Raman scattering cascades; and stimulated electron acoustic wave scattering (SEAS).