The resources needed for particle-in-cell simulations of laser wakefield acceleration can be greatly reduced in many cases of interest using an envelope model. However, the inclusion of tunneling ionization in this time-averaged treatment of laser-plasma acceleration is not straightforward, since the statistical features of the electron beams obtained through ionization should ideally be reproduced without resolving the high-frequency laser oscillations. In this context, an extension of an already known envelope ionization procedure is proposed, valid also for laser pulses with higher intensities, which consists in adding the initial longitudinal drift to the newly created electrons within the laser pulse ionizing the medium. The accuracy of the proposed procedure is shown with both linear and circular polarization in a simple benchmark where a nitrogen slab is ionized by a laser pulse and in a more complex benchmark of laser plasma acceleration with ionization injection in the nonlinear regime. With this addition to the envelope ionization algorithm, the main phase space properties of the bunches injected in a plasma wakefield with ionization by a laser (charge, average energy, energy spread, rms sizes, and normalized emittance) can be estimated with accuracy comparable to a nonenvelope simulation with significantly reduced resources, even in cylindrical geometry. Through this extended algorithm, preliminary studies of ionization injection in laser wakefield acceleration can be easily carried out even on a laptop.

The optimization and advanced study of a laser-plasma electron injector are presented based on a truncated ionization injection scheme for high quality beam production. The smilei code is used with laser envelope approximation and a low number of particles per cell to reach computation time performances enabling the production of a large number of accelerator configurations. The developed and tested workflow is a possible approach for the production of a large dataset for laser-plasma accelerator optimization. A selection of functions of merit used to grade generated electron beams is discussed. Among the significant number of configurations, two specific working points are presented in detail. All data generated are left open to the scientific community for further study and optimization.

Exploring novel target schemes for laser wakefield accelerators is essential to address the challenge of increasing repetition rates while ensuring the stability and quality of the produced electron beams. This paper introduces and discusses the prototyping of a two-chamber gas target, integrated into the beamline and operating under continuous gas flow, in the framework of ionization injection. We present the numerical fluid modeling employed to assist the density profile shaping, with a focus on gas mixing and dopant confinement. The importance of localized high-Z gas for ionization injection is demonstrated through particle-in-cell simulations using the simulated gas profiles. We describe the test bench used for prototype evaluation, specifically addressing the plasma electron density and the longitudinal distribution of species relevant to ionization injection. The lifetime of the target at 10 Hz and 60 mJ is measured for different materials, and its effect on the resulting electron beam is assessed using particle-in-cell simulations. Finally, we outline perspectives on high-power operation.

This study reports on a dayside magnetic reconnection event detected by the Magnetospheric Multiscale mission in the presence of a moderate guide field ( times the reconnecting magnetic field on the magnetosphere side) and assumed to be present in the whole reconnection region. The spacecraft traversed the compressed magnetospheric separatrix region, observing cold ions with densities up to 10 and a large magnetosheath density up to 150 . We provide a detailed analysis of current densities, generalized Ohm's law, and energy conversion processes in both the spacecraft and the fluid frames during the separatrix crossing. The normal electric field is directed away from the separatrix due to the cold ion drift on the magnetosphere side and to the magnetosheath ion drift in the presence of a guide field in the exhaust region. In the spacecraft frame, energy transfers from the plasma to the fields due to the convective field associated with the earthward motion of the magnetopause and the ion diamagnetic current associated with the large density gradient. In the fluid frame, energy conversion reverses due to the magnetic field‐aligned current density and electric field produced by the divergence of the electron pressure tensor . Additionally, we give insights into the local changes in electromagnetic, bulk flow, and thermal energies. We show that flow and thermal energy variations of the plasma are mostly driven by the compressible term of the electron pressure strain at the separatrix.

Future, high-fidelity simulations of multi-GeV-class laser Wakefield accelerators (LWFAs) will need to model the propagation of high-intensity laser drivers over meter-scale plasmas with high spatial and temporal resolutions, thus requiring high amounts of computational resources. Various techniques have been devised over the years to reduce the computational cost of such simulations, including the time-averaged ponderomotive approximation, and the use of the Lorentz boosted frame technique. In this paper we discuss the combination of these two computational techniques, highlighting the resulting significant reduction in the computational cost of LWFA simulations and the limitations of this approach. The combination of the two techniques can potentially become essential for the modeling of a multi-TeV, LWFA-based collider.

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 is a challenging task. In this paper, we first revisit ionization rates as given by the Smirnov-Chibisov and Perelomov-Popov-Terent'yev formulas within the paradigm of sequential tunnel ionization. We analyze the limit of validity and possible inconsistencies of this approach. We show that a strongly limiting factor to a precise description of ionization is the competing contribution of different sequential ionization processes. To solve this an algorithm is proposed that allows one to find the dominant nonsequential path of tunnel ionization and significantly improves the precision in simulations. This procedure is implemented in the PIC code SMILEI, and includes the dependence of the ionization rates on the magnetic quantum number of the level. The sensitivity to variations in the ionization model is studied via full simulations of the ionization of an argon target by an incident high-intensity laser pulse. Finally, we analyze generalizations of the Perelomov-Popov-Terent'yev rate developed to describe the barrier suppression ionization in high fields and discuss the necessity and possibility of including these extensions in PIC simulations.