Ion acceleration in underdense plasmas by ultra-short laser pulses

We report on the ion acceleration mechanisms that occur during the interaction of an intense and ultrashort laser pulse ( 10^{18} \text{W}\,\text{cm}^{-2}\,\mu \text{m}^{2}$?> I &lgr; 2 > 10 18 W cm − 2 &mgr; m 2 ?> ) with an underdense helium plasma produced from an ionized gas jet target. In this unexplored regime, where the laser pulse duration is comparable to the inverse of the electron plasma frequency &ohgr; pe ?> , reproducible non-thermal ion bunches have been measured in the radial direction. The two He ion charge states present energy distributions with cutoff energies between 150 and 200 keV, and a striking energy gap around 50 keV appearing consistently for all the shots in a given density range. Fully electromagnetic particle-in-cell simulations explain the experimental behaviors. The acceleration results from a combination of target normal sheath acceleration and Coulomb explosion of a filament formed around the laser pulse propagation axis.

Most of these works were performed by focusing the laser beam over a thin solid target, which allows us to get very sharp density ramps and high accelerating fields [2][3][4][5]10]. An alternative is to use gaseous targets, which enable work at a high repetition rate, to control the accelerated species and the plasma profile in a simple manner. This is the promising approach adopted in recent studies using very energetic infrared CO 2 lasers [17,18], for which suitable gas targets can be fabricated. However, the production of gas jets with overcritical density for near-infrared laser systems and sharp edges constitutes a technological challenge that has just started to be taken up [19], a fact that has limited the development of the gas target approach up to now.
Reported acceleration experiments using gas targets and near-infrared laser systems were performed in the underdense regime ( < n n e c , with n e and n c , respectively, the electron and critical plasma densities), and with laser pulse durations much larger than the inverse electron plasma frequency [1,[20][21][22]. In this regime, the fastest ions gain their radial momenta from Coulomb explosions, when the laser field is still present and bores a channel by expelling plasma electrons. The accelerated ions continue to drift undisturbed due to their inertia when the quasineutrality is restored as the electrons come back to screen the electrostatic field after the pulse is gone. In the first experiments, a 50 J-0.9 ps laser pulse focused over a gas jet with = × − n 5 10 cm e 19 3 allowed the production of a few MeV ions with exit angles around°100 to the axis of propagation and Maxwellian distribution [1]. The contrast between the obtained Maxwellian spectrum and the one expected from a Coulomb explosion (flat with a peak at the energy cutoff) was attributed to the superposition of Coulomb explosion-like spectra for the different values of laser intensity occurring along the pulse propagation in the gas jet [1,20]. Experiments performed with larger laser energies (180 J) and higher densities ( > − n 10 cm e 20 3 ) lead to ion energies of the order of 10 MeV, with plateau-like spectra [21]).
Longitudinal acceleration was achieved by raising the laser energy (340 J) [22][23][24][25]. In this case, a radial component with a cutoff of 10 MeV, arguably coming from the Coulomb explosion of the ion channel, coexisted with a beam accelerated longitudinally with a cutoff of 40 MeV. The longitudinal acceleration is caused by a target normal sheath acceleration (TNSA)-like mechanism, possibly enhanced by the azimuthal magnetic field existing at the plasma boundary [23][24][25].
In this paper, radial acceleration from underdense He plasma ( ∼ − n 10 cm 3 ) with short pulses is demonstrated for the first time. Observations with single-shot diagnostics of spectral patterns close to quasi-monoenergetic helium ion bunches (cutoff ∼ E 200 keV 0 ) are presented and discussed. Owing to the ultrashort laser pulse duration, comparable to the inverse electron plasma frequency, the direct laser-induced Coulomb explosion of a uniform positively-charged core from laser channeling is likely to be an incomplete picture for depicting the momentum transfer on the ion response timescale. The ions gain their radial momenta essentially after the pulse passage [26]. The short pulse excites a strongly nonlinear wakefield, which will radially accelerate ions and produce an annular radial ion density profile (see equation (8) in [27]). At high densities and high laser intensities, electrostatic effects result in the compression of ions situated close to the axis forming a filament, which later Coulomb explodes [28]. For a longer time, when the fast inner ions outrun the slow outer ions, a collisionless shock can form in the phase space, with significant modification of the final energy distribution [28,29].
In the experiment, an ultrashort (duration τ = 35 fs) linearly polarized pulse from the Laboratoire d'Optique Appliquée Ti:sapphire Salle Jaune laser [30] carrying = E 800 mJ ) was between 7 and 50 (the plasma and Rayleigh length values were comparable).
The interaction was spatially adjusted above the nozzle by means of plasma emission imaging ( μ 10 m spatial resolution). The top-view and side-view images allow the positioning of the beam, respectively, along a diameter of the jet and at a chosen height from the nozzle (see figure 1). Ion emission measurements were carried out at°0 and ∼°80 from the laser propagation axis through a 100 μm diameter pinhole (solid acceptance × − 1.2 10 sr 8 ) with a Thomson parabola spectrometer coupled to 40 mm imaging microchannel plate chevron assembly (two microchannel plates (MCPs) and a phosphor screen). This very narrow angle of collection efficiently filters out ions accelerated at different directions from the normal laser axis. A linearized thermo-electric-cooled 16-bit CCD camera ( × 1024 1024 px) enables singleshot recordings of scintillating traces and thus shot-to-shot statistical analysis of the ion spectra. The electric gains of the MCP assembly were set to preserve the full dynamics of 16-bit pictures from saturation and keep the signal-to-noise ratio better than 50:1. With a 2:1 magnification, the system gives a 500 μm pinhole image of the jet on the first MCP plane (10 pixels on the CCD image), which allows us to regard the source as point-like and the ion emission as laminar. The detection threshold is ∼48 keV for + He 2 and ∼10 keV for + He . No ion was detected in the forward direction (see figure 1), i.e. no + He 2 ion with energy above the detection threshold (48 keV) leaves the jet at°0 . Both ionic species, on the other hand, were detected with the spectrometer at°80 . A raw picture from a single laser shot is presented in figure 2 and shows striking gaps in the spectral traces and the presence of both species + He and  3 (0.5% to 2% of n c ). Moreover, changing the laser polarization direction does not affect the spectra. The number of ions reaching the MCP is relatively low (roughly of the order of a few hundreds) because of the very small angular acceptance of the detector.
To gain an insight into the mechanism of ion acceleration, we carry out a numerical study of the interaction process via particle-in-cell (PIC) simulations, performed using the fully electromagnetic PIC code Calder-Circ [32]. In this code, Maxwell equations are solved by projecting over a basis of Fourier modes. Few modes are enough to accurately describe the laser evolution and the wakefield, resulting in simulations needing much less computational resources than fully 3D simulations. The modelling of the acceleration process requires the inclusion in the simulation box of the whole gas jet, and the continuation of the simulations far beyond the wakefield decay time, precluding the use of too many numerically expensive Cartesian 3D simulations. The simulations performed here include two Fourier modes.
The  3(b)) nicely reproduces the main features of the experimental spectrum ( figure 3(a)). The spectra present low and high energy branches. The first one has a Maxwellian-like distribution with a cutoff around 40 keV. The high energy branch presents a broad peak, which goes from 50 to 150 keV for an exit angle of°78 . These values are close to the low and high energy cutoffs of the upper branch in the experiment. Simulation also shows that ions leaving the jet at°0 have energies lower than 30 keV. This value lies below the detection threshold for + He 2 (48 keV), but above the one for + He (10 keV). However, their conversion to + He , which would allow their detection by the Thomson parabola at°0 , is precluded because these ions do not transverse a region of neutral or partially ionized gas. The simulated spectrum does not visibly change by changing the laser polarization direction, in agreement with the experimental result.
The two energy branches found in the spectrum correspond to two distinct ion populations accelerated by different mechanisms and in different regions of the gas jet. The low energy branch is composed of ions accelerated in the central region of the gas jet by the radial field existing in the wakefield. The high energy branch, on the other hand, corresponds to ions accelerated in the plasma sheath that is created in the falling density ramp of the plasma. As the longitudinal plasma density gradient is not sharp enough, the radial field is larger than the longitudinal one, and the main acceleration is radial. Ions are pushed towards the axis by the quasi-electrostatic field in the sheath, and Coulomb explosion of the formed filament further accelerates the ions.
The acceleration of the low energy ion population results from the fact that the average radial field in a non-linear wakefield does not vanish. Close to the axis (typically for μ < r 1 m) the time averaged radial field is negative, especially in the presence of electrons injected and accelerated in the wakefield [28]. The Thomson parabola at°0 was not able to detect these electrons because its magnetic field is too large. However, a test made using a double magnet scheme shows the presence of relativistic electrons longitudinally accelerated in the wakefield, although the setup was not well suited to recovering the electron spectrum.
Beyond the central region, the average field over the ions in the pulse trail is positive. The result of the central focusing field and that surrounded by a region of the defocusing field is the formation of a channel of low density with a central filament. The channel radially expands roughly ballistically up to the formation of a sharp edge on the outer shell of the channel, which cannot be screened, and the emergence of an electrostatic shock [28,29], which will redistribute the energy between ions. Simulations made with a very narrow gas jet and meshes with higher radial resolution show that to accurately describe the effect of the shock on the spectrum requires the use of smaller radial sizes. However, these simulations show that the shock has little effect on the cutoff energy itself, which will remain around 30 keV. The filament, on the other hand, will pinch ballistically up to a sub-micrometric radius, and then Coulomb explode. The cutoff energy after the explosion is close to 40 keV.
The high energy branch of the spectrum corresponds to ions accelerated in the falling density ramp of the gas jet, where a TNSA-like sheath is formed. The dilute electron population is not able to screen the ion charge, and a radially focusing and longitudinally accelerating (for positive ions) electric field will arise (figures 4(a) and (c)) and stand for a few ps. The presence  of a stable magnetic structure, visible in the simulation, could also contribute to the sustainability of the sheath by limiting the electron transport [23][24][25]. The radial field focuses ions coming from μ < r 10 m towards the axis, forming a filament ( figure 5(b)). When the filament reaches its minimum radius, its peak density becomes 20 times larger than the unperturbed original density, whereas for the filament formed in the plasma core, the maximum density is only two times larger than the unperturbed one. The weak screening will result in a stronger Coulomb explosion (figures 5(c) and 4(b)). The result of this explosion is a sudden radial acceleration with an energy cutoff of ∼150 keV. Due to the longitudinal field ( figure 5(b)), ions will leave the target with an angle smaller than°90 . The longitudinal acceleration also results in a finite minimal energy for ions collected at a given angle, and therefore the gap between the spectra of the two ion populations. Ions above 50 keV are emitted with angles between ∼°55 and ∼°85 , symmetrically distributed around the laser propagation axis. The solid angle of the resulting emission is around π sr. Because of this very broad emission, the characterization of the accelerated ion population in terms of beam properties (emittance, length, etc) is not very meaningful.
For densities higher than n 2% c and up to the highest density explored ( n 3.5% c ), the gap disappears and the spectrum becomes Maxwellian-like, a trend also found in simulations. The reason for this is the rise in the energy cutoff of ions accelerated in the plasma core, which becomes higher than 50 keV, thus filling the gap. The spectrum of the two ion populations are therefore overlapped. These results will be detailed in future work.
In conclusion, novel aspects of laser ion acceleration have been revealed in a regime where the laser pulse duration is close to ω − p 1 . An ion population radially accelerated in the falling density ramp of the gas jet was identified. The energy gap between this population and a less energetic one accelerated in the wakefield is a signature of the presence of a longitudinal field in the falling density ramp. This study also suggests that to get more energetic ions in the forward direction, one has to enhance the longitudinal electric field, for example by using targets with much steeper density gradients.