Search

We discuss the inclusion of the process of multiple ionization of atoms in high-intensity electromagnetic fields into particle-in-cell (PIC) codes applied to the simulation of laser-plasma interactions. The approach is based on the application of the Perelomov-Popov-Terent'yev formulas for ionization rates within the paradigm of sequential tunnel ionization. We analyze difficulties and possible inconsistencies of the approach, which may cause problems in the implementation of PIC codes. To test the sensitivity to the ionization model variations, we perform simulations of an argon target ionization by an incident high-intensity laser pulse. For this, we implement an original procedure in PIC code SMILEI, which allows us to account for the dependence on the magnetic quantum number in ionization rates consistently, and also include barrier suppression at high fields., Comment: 14 pages, 3 figures, 3 tables

Autoresonant phase-locking of the plasma wakefield to the beat frequency of two driving lasers offers advantages over conventional wakefield acceleration methods, since it requires less demanding laser parameters and is robust to variations in the target plasma density. Here, we investigate the kinetic and nonlinear processes that come into play during autoresonant plasma beat-wave acceleration of electrons, their impact on the field amplitude of the accelerating structure, and on acceleration efficiency. Particle-in-Cell simulations show that the process depends on the plasma density in a non-trivial way but can be reliably modeled under specific conditions. Beside recovering previous fluid results in the deeply underdense plasma limit, we demonstrate that robust field excitation can be achieved within a fully kinetic self-consistent modeling. By adjusting the laser properties, we can amplify the electric field to the desired level, up to wave-breaking, and efficiently accelerate particles; we provide suggestions for optimized laser and plasma parameters. This versatile and efficient acceleration scheme, producing electrons from tens to hundreds of MeV energies, holds promise for a wide range of applications in research industry and medicine., Comment: 14 pages, 13 figures

Electromagnetic showers developing from the collision of an ultra-intense laser pulse with a beam of high-energy electrons or photons are investigated under conditions relevant to future experiments on multi-petawatt laser facilities. A semi-analytical model is derived that predicts the shower multiplicity, i.e. the number of pairs produced per incident seed particle (electron or gamma photon). The model is benchmarked against particle-in-cell simulations and shown to be accurate over a wide range of seed particle energies (100 MeV - 40 GeV), laser relativistic field strengths ($10 < a_0 < 1000$), and quantum parameter $\chi_0$ (ranging from 1 to 40). It is shown that, for experiments expected in the next decade, only the first generations of pairs contribute to the shower while multiplicities larger than unity are predicted. Guidelines for forthcoming experiments are discussed considering laser facilities such as Apollon and ELI Beamlines. The difference between electron- and photon seeding and the influence of the laser pulse duration are investigated., Comment: 23 pages, 9 figures

It was suggested [A. R. Bell & J. G. Kirk, PRL 101, 200403 (2008)] that an avalanche of electron-positron pairs can be triggered in the laboratory by a standing wave generated by intense laser fields. Here, we present a general solution to the long-standing problem of the avalanche growth rate calculation. We provide a simple formula that accounts for the damping of the growth rate due to pair and photon migration from the region of prolific generation. We apply our model to a variety of 3D field configurations including focused laser beams and show that i) the particle yield for the full range of intensity able to generate an avalanche can be predicted, ii) a critical intensity threshold due to migration is identified, iii) the effect of migration is negligible at a higher intensity and the local growth rate dominates. Excellent agreement with Monte Carlo and self-consistent PIC simulations is shown. The growth rate calculation allows us to predict when abundant pair production will induce a back-reaction on the generating field due to plasma collective effects and screening. Our model can be applied to study the generation of electron-positron pair avalanches in realistic fields to plan future experiments at ultra-high-intensity laser facilities., Comment: 32 pages, 12 figures, 3 tables

We propose a new mechanism for generating coherent subcycle optical pulse by directing a relativistic electron beam (REB) into a plasma with a density up-ramp. The subcycle pulse is coherently emitted by bubble-sheath electrons in REB-induced superluminal plasma wake. Using three-dimensional particle-in-cell and far-field time-domain radiation simulations as well as analytical modeling, we show that an isolated subcycle optical shock can be produced at the Cherenkov angle. This radiation has ultra-short attosecond-scale duration and high intensity and exhibits excellent directionality with ultra-low angular divergence and stable carrier envelope phase. Its central frequency can be easily tuned over a wide range, from the far-infrared to the ultra-violet, by adjusting the plasma and driver-beam density., Comment: 7 pages, 3 figures

We present a novel approach named as pumping ionizing gating (PIG) for the generation of isolated attosecond pulses (IAPs). In this regime, a short laser is used to ionize a pre-existing gas grating, creating a fast-extending plasma grating(FEPG) having an ionization front propagating with the velocity of light. A low-intensity long counterpropagating pump pulse is then reflected by a very narrow region of the ionization front, only where the Bragg conditions for resonant reflection is satisfied. Consequently, the pump reflection is confined within a sub-cycle region called PIG, and forms a wide-band coherent IAP in combination with the frequency up-conversion effect due to the plasma gradient. This approach results in a new scheme to generate IAPs fromlong picosecond pump pulses. Three-dimensional (3D) simulations show that a 1.6-ps, 1-{\mu}m pump pulse can be used to generate a 330 as laser pulse with a peak intensity approximately 33 times that of the pump and a conversion efficiency of around 0.1%.These results highlight the potential of the PIG method for generating IAPs with high conversion efficiency and peak intensity., Comment: It provides a new way to generate isolated attosecond pulse(IAP) by a picosecond pump, which has a protential to boost the IAP energy to joule level

A new electron acceleration mechanism is identified that develops when a relativistically intense laser irradiates the wedge of an over-dense plasma. This induces a diffracted electromagnetic wave with a significant longitudinal electric field that accelerates electrons from the plasma over long distances to relativistic energies. Well collimated, highly-charged (nC) electron bunches with energies up to 100's MeV are obtained using a laser beam with $I \lambda_0^2 =3.5\times 10^{19}\,{\rm W \mu m^2/cm^2}$. Multi-dimensional particle-in-cell simulations, supported by a simple analytical model, confirm the efficiency and robustness of the proposed acceleration scheme.

Ultra-short high-power lasers can deliver extreme light intensities ($\ge 10^{20}$ W/cm$^2$ and $\leq 30 f$s) and drive large amplitude Surface Plasma Wave (SPW) at over-dense plasma surface. The resulting current of energetic electron has great interest for applications, potentially scaling with the laser amplitude, provided the laser-plasma transfer to the accelerated particles mediated by SPW is still efficient at ultra-high intensity. By mean of Particle-in-Cell simulations, we identify the best condition for SPW excitation and show a strong correlation between the optimum Surface Plasma Wave excitation angle and the laser's angle of incidence that optimize the electron acceleration along the plasma surface. We also discuss how plasma density and plasma surface shape can be adjusted in order to push to higher laser intensity the limit of Surface Plasma Wave excitation. Our results open the way to new experiments on forthcoming multi-petawatt laser systems.

The forthcoming generation of multi-petawatt lasers opens the way to abundant pair production by the nonlinear Breit-Wheeler process, i.e., the decay of a photon into an electron-positron pair inside an intense laser field. In this paper we explore the optimal conditions for Breit-Wheeler pair production in the head-on collision of a laser pulse with gamma photons. The role of the laser peak intensity versus the focal spot size and shape is examined keeping a constant laser energy to match experimental constraints. A simple model for the soft-shower case, where most pairs originate from the decay of the initial gamma photons, is derived. This approach provides us with a semi-analytical model for more complex situations involving either Gaussian or Laguerre-Gauss (LG) laser beams. We then explore the influence of the order of the LG beams on pair creation. Finally we obtain the result that, above a given threshold, a larger spot size (or a higher order in the case of LG laser beams) is more favorable than a higher peak intensity. Our results match very well with three-dimensional particle-in-cell simulations and can be used to guide upcoming experimental campaigns., Comment: 34 pages, 11 figures, submitted to New J. Phys

Type-III-burst radio signals can be mimicked in the laboratory via laser-plasma interaction. Instead of an electron beam generating Langmuir waves (LW) in the interplanetary medium, the LWs are created by a laser interacting with a millimeter-sized plasma through the stimulated Raman instability. In both cases, the LWs feed the Langmuir decay instability which scatters them in several directions. The resulting LWs may couple to form electromagnetic emission at twice the plasma frequency, which has been detected in the interplanetary medium, and recently in a laboratory laser experiment [Marqu\`es et al. Phys. Rev. Lett. 124, 135001 (2020)]. This article presents the first numerical analysis of this laser configuration using particle-in-cell simulations, providing details on the wave spectra that are too difficult to measure in experiments. The role of some parameters is addressed, with a focus on laser intensity, in order to illustrate the behavior of the electromagnetic emission's angular distribution and polarization.

We propose to use ultra-high intensity laser pulses with wavefront rotation (WFR) to produce short, ultra-intense surface plasma waves (SPW) on grating targets for electron acceleration. Combining a smart grating design with optimal WFR conditions identified through simple analytical modeling and particle-in-cell simulation allows to decrease the SPW duration (down to few optical cycles) and increase its peak amplitude. In the relativistic regime, for $I\lambda_0^2=3.4 \times 10^{19}{\rm W/cm^2\mu m^2}$, such SPW are found to accelerate high-charge (few 10's of pC), high-energy (up to 70 MeV) and ultra-short (few fs) electron bunches.

A new experimental platform based on laser-plasma interaction is proposed to explore the fundamental processes of wave coupling at the origin of interplanetary radio emissions. It is applied to the study of electromagnetic (EM) emission at twice the plasma frequency ($2\omega_p$) observed during solar bursts and thought to result from the coalescence of two Langmuir waves (LWs). In the interplanetary medium, the first LW is excited by electron beams, while the second is generated by electrostatic decay of Langmuir waves. In the present experiment, instead of an electron beam, an energetic laser propagating through a plasma excites the primary LW, with characteristics close to those at near-Earth orbit. The EM radiation at $2\omega_p$ is observed at different angles. Its intensity, spectral evolution and polarization confirm the LW-coalescence scenario., Comment: Accepted for publication in Phys. Rev. Lett

Laser-generated plasma gratings are dynamic optical elements for the manipulation of coherent light at high intensities, beyond the damage threshold of solid-stated based materials. Their formation, evolution and final collapse require a detailed understanding. In this paper, we present a model to explain the nonlinear dynamics of high amplitude plasma gratings in the spatially periodic ponderomotive potential generated by two identical counter-propagating lasers. Both, fluid and kinetic aspects of the grating dynamics are analyzed. It is shown that the adiabatic electron compression plays a crucial role as the electron pressure may reflect the ions from the grating and induce the grating to break in an X-type manner. A single parameter is found to determine the behaviour of the grating and distinguish three fundamentally different regimes for the ion dynamics: completely reflecting, partially reflecting/partially passing, and crossing. Criteria for saturation and life-time of the grating as well as the effect of finite ion temperature are presented.

Laser plasma amplification of sub-picosecond pulses above the Joule level is demonstrated, a major milestone for this scheme to become a solution for the next-generation of ultra-high intensity lasers. By exploring over 6 orders of magnitude the influence of the incident seed intensity on Brillouin laser amplification, we reveal the importance of a minimum intensity to ensure an early onset of the self-similar regime, and a large energy transfer with a very high efficiency, up to 20%. Evidence of energy losses of the seed by spontaneous backward Raman is found at high amplification. The first three-dimensional envelope simulations of the sub-picosecond amplification were performed, supplemented by one-dimensional PIC simulations. Comparisons with the experimental results demonstrate the capability of quantitative predictions on the transferred energy. The global behavior of the amplification process, is reproduced, including the evolution of the spatial profile of the amplified seed.

Soon available multi petawatt ultra-high-intensity (UHI) lasers will allow us to probe high-amplitude electromagnetic fields interacting with either ultra-relativistic electron beams or hot plasmas in the so-called moderately quantum regime. The correct modelling of the back-reaction of high-energy photon emission on the radiating electron dynamics, a.k.a. radiation reaction, in this regime is a key point for UHI physics. This will lead to both validation of theoretical predictions on the photon spectrum emitted during the laser-particle interaction and to the generation of high energy photon sources. In this paper we analyse in detail such emission using recently developed models to account for radiation reaction. We show how the predictions on the spectrum can be linked to a reduced description of the electron distribution function in terms of the first energy moments. The temporal evolution of the spectrum is discussed, as well as the parameters for which quantum effects induce hardening of the spectrum.

Radiation-reaction in the interaction of ultra-relativistic electrons with a strong external electromagnetic field is investigated using a kinetic approach in the weakly quantum regime ($\chi \lesssim 1$, with $\chi$ the electron quantum parameter). Three complementary descriptions are considered, their domain of applicability discussed and their predictions on average properties of an electron population compared. The first description relies on the radiation reaction force in the Landau and Lifschitz (LL) form. The second relies on the linear Boltzmann equation for the electron and photon distribution functions. It is valid for any $\chi \lesssim 1$, and usually implemented numerically using a Monte-Carlo (MC) procedure. The third description relies on a Fokker-Planck (FP) expansion and is rigorously derived for any ultra-relativistic, otherwise arbitrary configuration. Our study shows that the evolution of the average energy of an electron population is described with good accuracy in many physical situations by the leading term of the LL equation with the so-called quantum correction, even for large values of the $\chi$. The leading term of the LL friction force (with quantum correction) is actually recovered naturally by taking the FP limit. The FP description is necessary to correctly describe the evolution of the energy variance (second order moment) of the distribution function, while the full linear Boltzmann (MC) description allows to describe the evolution of higher order moments whose contribution can become important when $\chi \rightarrow 1$. This analysis allows further insight on the effect of particle straggling in the deformation of the particle distribution function. A general criterion for the limit of validity of each description is proposed, as well as a numerical scheme for inclusion of the FP description in Particle-In-Cell codes.

SMILEI is a collaborative, open-source, object-oriented (C++) particle-in-cell code. To benefit from the latest advances in high-performance computing (HPC), SMILEI is co-developed by both physicists and HPC experts. The code's structures, capabilities, parallelization strategy and performances are discussed. Additional modules (e.g. to treat ionization or collisions), benchmarks and physics highlights are also presented. Multi-purpose and evolutive, SMILEI is applied today to a wide range of physics studies, from relativistic laser-plasma interaction to astrophysical plasmas., Comment: submitted to Computer Physics Communications

The Weibel instability driven by two symmetric counter-streaming relativistic electron plasmas, also referred to as current-filamentation instability, is studied in a constant and uniform external magnetic field aligned with the plasma flows. Both the linear and non linear stages of the instability are investigated using analytical modeling and Particle-In-Cell (PIC) simulations. While previous studies have already described the stabilizing effect of the magnetic field, we show here that the saturation stage is only weakly affected. The different mechanisms responsible for the saturation are discussed in detail in the relativistic cold fluid framework considering a single unstable mode. The application of an external field leads to a slighlt increase of the saturation level for large wavelengths, while it does not affect the small wavelengths. Multi-mode and temperature effects are then investigated. While at large temperature the saturation level is independent of the external magnetic field, at small but finite temperature the competition between different modes in the presence of an external magnetic field leads to a saturation level lower with respect to the unmagnetized case., Comment: Submitted to Phys. Rev. E

Different electron acceleration regimes in the evanescent field of a surface plasma wave are studied by considering the interaction of a test electron with the high-frequency electromagnetic field of a surface wave. The non-relativistic and relativistic limits are investigated. Simple scalings are found demonstrating the possibility to achieve an efficient conversion of the surface wave field energy into electron kinetic energy. This mechanism of electron acceleration can provide a high-frequency pulsed source of relativistic electrons with a well defined energy. In the relativistic limit, the most energetic electrons are obtained in the so-called electromagnetic regime for surface waves. In this regime the particles are accelerated to velocities larger than the wave phase velocity, mainly in the direction parallel to the plasma-vacuum interface.

The interaction of laser pulses with thin grating targets, having a periodic groove at the irradiated surface, has been experimentally investigated. Ultrahigh contrast ($\sim 10^{12}$) pulses allowed to demonstrate an enhanced laser-target coupling for the first time in the relativistic regime of ultra-high intensity $>10^{19} \mbox{W/cm}^{2}$. A maximum increase by a factor of 2.5 of the cut-off energy of protons produced by Target Normal Sheath Acceleration has been observed with respect to plane targets, around the incidence angle expected for resonant excitation of surface waves. A significant enhancement is also observed for small angles of incidence, out of resonance., Comment: 5 pages, 5 figures, 2nd version implements final corrections

The generation of astrophysically relevant jets, from magnetically collimated, laser-produced plasmas, is investigated through three-dimensional, magneto-hydrodynamic simulations. We show that for laser intensities I ~ 10^12 - 10^14 W/cm^2, a magnetic field in excess of ~ 0.1 MG, can collimate the plasma plume into a prolate cavity bounded by a shock envelope with a standing conical shock at its tip, which re-collimates the flow into a super magneto-sonic jet beam. This mechanism is equivalent to astrophysical models of hydrodynamic inertial collimation, where an isotropic wind is focused into a jet by a confining circumstellar torus-like envelope. The results suggest an alternative mechanism for a large-scale magnetic field to produce jets from wide-angle winds. (abridged version), Comment: Accepted for publication in Physical Review Letters. 5 pages, 4 figures