Introduction

Laser-acceleration of ions is currently a topic of high relevance for a wide range of possible applications1,2,3 and associated with many fundamental processes in the laboratory and in astrophysics4,5. In particular, ion acceleration processes in laser-driven scenarios are being studied to realise compact laser-based particle accelerators capable of producing ion beams of high quality. Understanding and controlling the different ion acceleration mechanisms6,7,8,9,10,11, as well as their interplay, is of paramount importance for future developments. Therefore, it is critical to identify scenarios where this control and detailed analysis can be performed, both from a fundamental perspective but also for experimental conditions now being pursued12,13,14,15,16,17.

Inspired by previous work on the formation of shock shells in nano plasmas18,19, we propose a laser irradiation of a macroscopic finite size target to tailor the phase space of the accelerated ions by modifying the target density in an appropriate way. In this paper, we show that the shaping of the ion phase space can be achieved by triggering and combining different acceleration mechanisms. By using a first laser pulse for a pre-heating stage and a second pulse to actually drive the acceleration process19,20, one initial target setup can be used to enter the different acceleration regimes, by changing the delay between the two pulses and the intensity of the first pulse.

A first low-intensity laser pulse with normalised laser vector potential a0 = eA0/mec2 < 1, where e is the electron charge, A0 the laser vector potential, me the electron mass and c the speed of light, irradiates a spherical solid density target (Fig. 1a). The electron energy increases which is followed by a hydrodynamic expansion of the target (Fig. 1b) with the ion sound speed , where Z is the ion charge, kB the Boltzmann constant, Te the electron temperature and mi the ion mass21,22. For laser intensities below 1016 W/cm2, corresponding to , collisional energy absorption processes dominate (inverse bremsstrahlung), while for larger intensities collisionless processes (resonance absorption, Brunel heating) are more important. For simplicity, a uniform target density n0 after expansion is assumed with a radius increase from R* to R by a factor . The expansion time at which the target density has been reduced by a factor x is then given by

Figure 1
figure 1

Simulation setup: (a) The high-density target is irradiated by a laser with large focal area (b) leading to hydrodynamic expansion of the target. (c) A second, stronger laser pulse launches (d) the acceleration process.

We expressed the temperature in terms of the laser potential according to ponderomotive heating, 23 for a small argument expansion . This model is rather simple, e. g. neglecting the density ramps at the target surface or non-uniform target heating and expansion, but it does capture the main questions addressed in this paper.

The actual ion acceleration is then triggered by a second high intensity laser pulse incident on the pre-heated target (Fig. 1c,d). The dominant acceleration process depends on the target density at the time of arrival of the main pulse. Thus, the correct timing of the interaction of the second laser pulse with the target is critical to guarantee that collisionless shock acceleration can prevail over target-normal-sheath-acceleration (TNSA), which is always present at the surface of overdense targets (compare also the 1D and 2D simulations in refs 24 and 25 which show the efficiency of collisionless shock acceleration towards TNSA). We now review the main acceleration mechanisms that can play a role in this configuration.

Typically in the TNSA scheme, μm scale overdense targets are irradiated by a high-intensity laser pulse. The laser cannot penetrate the target beyond densities greater than the relativistic critical density , with the correction due to the relativistic gamma factor of the electrons and the laser frequency ω0. The laser pulse is partially absorbed at the target front surface accelerating electrons to relativistic energies due to ponderomotive heating. These hot electrons propagate across the target and escape at the rear side. The subsequent electrostatic field at the target rear side accelerates the positive ions along the target normal direction to an energy of roughly Ei,TNSA = ZEe9,23. The highest proton energies experimentally achieved with TNSA are 85 MeV cite (F. Wagner et al., PRL 116 205002, 2016)8. TNSA is very sensitive to the particular laser and target configurations and the initial setup26. The ion energy spectrum can typically be described by a thermal-like exponentially decreasing distribution which leads to only a small number of accelerated particles close to the maximum energy cutoff.

Even if the laser cannot propagate through overdense targets, the ponderomotive force is strong enough to push the electrons at the target front surface. As a consequence and since the ions remain quasi-stationary due to their higher inertia, a quasi-static electric field is created at the front surface oriented into the target and propagating at the hole boring velocity with 9,23. Ions reflected from the hole boring potential will acquire a velocity into the target of approximately , which determines the scaling for hole boring acceleration (HBA).

In near-critical density targets, if the hole boring velocity is larger than the ion sound speed, vhb > cs, which is equivalent to the condition n0/ncr < a0 for , the propagation of a shock can be observed ahead of the hole boring compression with a velocity and the adiabatic coefficient κad = 5/3 from ideal 3D gas theory27. A lower limit is obtained from the condition that the density of the downstream (the compressed plasma moving into the target) has to be larger than the density of the upstream (the plasma target at rest with initial density n0), 27. Also in this case, the ions are picked up by the propagating potential and accelerated with approximately twice the shock velocity to , resulting in an energy Ei,sh = mic2i,sh − 1) with 12. If the shock velocity is non-relativistic and , the final ion energy can be approximated as

cf.28,29. Introducing y = n0/γncr ≈ n0/(a0ncr), the ion energy is given as

We note that in general vsh > vhb and thus electrostatic collisionless shock acceleration (CSA) will always lead to higher ion energies than HBA.

If the target density is further decreased to underdense conditions, i.e. , the laser can penetrate the target. However, if the target is thin enough and the laser pulse is long enough that the laser can break through the rear surface, a significant increase in energy is expected from an enhanced TNSA effect, which we also observe in our simulations.

Results

This framework can now be used to interpret and to explore the proposed scenario of ion acceleration. We illustrate this set-up with full scale three dimensional particle-in-cell (PIC) simulations modelling the interaction of a laser pulse with a macroscopic finite size target. The simulations were performed with the PIC code OSIRIS30,31. We model a full realistic experimental setup in three dimensions, which allows us to determine quantitatively all the relevant beam parameters, i.e. (maximum) proton energy, the spectrum and distributions. The initial setup consists of a spherical target with R* = 10 μm being the radius of the sphere, with homogeneous initial density corresponding to a frozen hydrogen pellet with mp/me = 1836. Different target shapes (spherical and cubic) and initial radii have also been tested maintaining the same volume and total number of particles. Additionally, a cylindrical jet with infinite extension along an axis perpendicular to the laser propagation direction has been simulated. For simplicity, the simulation starts with the interaction of the second high intensity laser pulse, which is injected along x1 with normalised intensity in the range a0 = 2–47 with a transverse Gaussian profile and with linear polarisation in x2 direction. We also modelled the full set-up, including the interaction of the first laser pulse with the cold, spherical target, leading to a density gradient in the target density, to the interaction of a second laser pulse with the self-consistently expanded target. We found that the observed lower density at the target surface leads only to a small shift on the reference density for the different regimes (cf. Fig. 2). Since the impact of the inclusion of the self-consistent expansion on the dominant acceleration mechanism is very low, we adopted the simplest model possible, i.e. a spherical target with uniform density characterised by just the uniform density parameter.

Figure 2
figure 2

Maximum proton energies expected from theory as function of the relativistic target density for a0 = 10 due to different acceleration mechanisms: target normal sheath acceleration (green), collisionless shock acceleration (blue), hole boring (brown).

The time axis corresponds to the time delay between the heating pulse with a0 = 0.1 initiating a hydrodynamic expansion and the second laser pulse. The maximum proton energies in the simulations measured inside the target are compared for a cubed (black squares) and a spherical target (red dots).

In our simulations, the parameters of the second laser pulse are the laser wavelength λ0 = 1 μm = 2π c0 with a pulse length , focused on the target front side. The spot size is w0 = 3.5 μm for a0 = 10 and is scaled as a0w0 = const. according to with an approximated focal area of . These are typical parameters of current laser systems32. The cubic simulation box was adapted to the increasing target volume with a maximum box length of L = 855 c0 and 2760 cells per dimension with two particles per cell. The cell size was chosen to guarantee a sampling frequency of at least two points per characteristic wavelength in the overdense plasma, . For a detailed list of the parameters, see Tables 1, 2, 3.

Table 1 Simulations with a0 = 10: Density after hydrodynamic expansion n0/γncr, radius of the sphere R or alternative box length of the cube Lc, simulation box length Lbox and simulation cell size Δx for a laser wavelength λ = 1 μm.
Table 2 Simulations with a0 = 47: Density after hydrodynamic expansion n0/γncr, radius of the sphere R or alternative box length of the cube Lc, simulation box length Lbox and simulation cell size Δx for a laser wavelength λ = 1μm.
Table 3 Simulations with n0/ncr = 6.7, box length Lbox = 72.29 μm, cell size Δx = 0.0277 μm and sphere radius R = 18.15 μm: Normalised laser potential a0 and spot size w0 in SI- and normalised units.

The central results of the parameter scan with constant a0 = 10 (Table 1) performed to illustrate the ion acceleration set-up are presented in Fig. 2. The highest energies in the simulations match with the theoretical prediction. TNSA dominates in the regime of highly overdense targets, which in this setup generates protons with ≈5 MeV maximum energy. The results seem to be independent of the actual target shape (cube or sphere). The maximum proton energy is observed at near-critical density due to CSA and is highest (Ep ≈ 25 MeV) for a target density n0/γncr ≈ 0.5. We observe the same trend as in Fig. 2, with a peak energy for shock acceleration at n0/γncr ≈ 0.5, if the initial density is kept fixed and the laser intensity is varied (parameters given in Table 3). In the supplementary material the proton phase spaces of representative cases are shown in Fig. S1.

Discussion

A detailed analysis of the overdense regime n0/γncr ≥ 0.5 shows a strong density compression propagating through the target. In the pure TNSA-regime n0/γncr > 2, the associated electric field is not strong enough to accelerate the ions inside the target and acceleration occurs mainly at the target rear surface (see Fig. S1a).

In the case of near-critical densities 0.5 ≤ n0/γncr ≤ 2 a strong peak of the electric field is formed at the target front which is co-moving with the region of density compression. It picks up the upstream protons and accelerates them to approximately twice the peak velocity, Eq. (3). This is illustrated in Fig. 3(a) for n0/(γncr) = 1, where two velocities can be clearly distinguished, the hole boring velocity vhb = 0.0522 c and the shock velocity vsh = 0.0696 c > vhb. Due to the thermal expansion of the target, the protons have an additional velocity v0 = 0.02 c, which slightly increases their final velocity to ≈2 vsh + v0. The electrons are heated to approximately 0.8 MeV, thus the shock Mach number is M = vsh/cs ≈ 2.4. The bunch of accelerated particles can be clearly identified in the ion phase space in Fig. S1(b), where collisionless shock acceleration at the target front dominates over TNSA at the target rear surface.

Figure 3
figure 3

Density (a) and magnetic field (b) evolution along x1 at x2 = x3 = 200 c0 for initial target density n0/γncr = 1, laser intensity a0 = 10 and . The grey and black lines show the hole boring velocity, vhb and shock velocity, vsh, respectively. The green dashed lines have a slope v = 0.01 c and the inset shows B3 at 0 = 806.

A collimating azimuthal magnetic field is built-up, see inset of Fig. 3(b) showing B3 (a similar picture is obtained for B2 with B2 and B3 being the components of the azimuthal magnetic field). It is related to the balance of the hot electron beam through the target by a fast cold return current I = recB/2 = 12 mec3/e with a maximum Larmor radius of re = 16 c0 for the hottest electrons and a maximum magnetic field strength B ≈ 1.5 me0/e33,34. The magnetic field guides the electrons in the direction of the back of the target at a constant velocity v = 0.01 c after the laser is off, see Fig. 3(b). We stress that this mechanism to enhance and to guarantee the conditions for shock propagation was not previously identified in the literature. The density compression is maintained by the collimation of the electron beam.

We note that for densities 0.5 < n0/γncr < 1 the compression at the target surface due to the radiation pressure of the laser makes the target opaque to the laser and shock acceleration inside the target is still observed. In the underdense scenario, n0/γncr < 0.5, this compression is not strong enough and the laser can propagate. Proton acceleration inside the target is rather weak in this case.

In addition to the ion energy, it is also important to consider the total number of accelerated protons to the relevant energy ranges. Fully quantitative results can only be obtained in three dimensional simulations. Figure 4 shows, for a0 = 10 and varying target density, the proton spectra and the number of protons with energies above a threshold energy Eth. The clear dominance of shock acceleration can also be recognised here. There are still more than 108 protons at energies E > 15 MeV (see inset in Fig. 4), which would guarantee a sufficient number of particles for many applications35,36,37. See also the comparison of the beam divergence (Fig. S2) which shows the advantage of collisionless shock acceleration towards TNSA. The TNSA case shows an almost uniform distribution of the protons while the CSA case shows a narrow distribution at low angles. Figure 5 shows that the proton energy is saturated during the quasi-steady state of shock propagation. We calculated the energy conversion efficiency, Ekin,total/Elaser, from the laser to protons with energies >4 MeV being accelerated within an angle θ < 10° where θ is defined as tan θ = p/p with p and p the transverse and parallel components of the proton momenta, respectively. The energy conversion is 6–7 orders of magnitude higher for CSA than for pure TNSA.

Figure 4
figure 4

(a) Proton energy spectra for a0 = 10 at t = 0.8 ps and varying target densities given in y = n0/γncr. (b) Number of protons with energy above a threshold energy Eth vs. threshold energy and target density n0/γncr.

Figure 5
figure 5

Maximum proton energy vs. time.

It has to be mentioned that in 3D simulations the maximum energies are significantly smaller and overall TNSA acceleration is much weaker than in a comparable 2D run as was previously clarified in38. This is due to the different radial dependence of the accelerating electric field and the lower efficiency of the heating mechanism in the case of a linearly polarised laser.

We also tested the robustness of our study to varying the laser intensity and pulse length, also with a detailed analysis of the PW regime for and , parameters see Table 2. We found that the discussion presented before holds, showing a similar scaling behaviour with the same dominant acceleration mechanisms as in the case for a0 = 10, see Fig. 6. The maximum proton energy is increased by a factor according to Eq. 3. (Note: We also explored the regime of radiation pressure acceleration (RPA) for highly compressed and thin targets. The target is destroyed and almost all the ions are accelerated, see Fig. S1c). The overall higher proton energy has its maximum Ep = 124 MeV at n0/γncr ≈ 0.5 in the overdense regime due to CSA. Thus, this scheme seems to be robust in order to explore the different acceleration mechanisms and energy scalings at higher laser intensities and is promising for achieving proton energies above 100 MeV in a realistic experimental setup with 10 PW laser systems.

Figure 6
figure 6

Maximum proton energy vs. varying target density n0γncr: Simulation with (blue squares) compared to case a0 = 10 (cf. Fig. 2) showing .

In conclusion, we proposed a scheme that allows to test and to tune the dominant process for ion acceleration. An initially highly opaque spherical target is heated by a low intensity laser pulse leading to a reduction of the target density by a hydrodynamic expansion. By timing a second, intense driver pulse, the acceleration process is initiated. Our fully relativistic 3D particle-in-cell simulations of a realistic experimental setup show a maximum proton energy close to critical density, n0/γncr = 0.5, due to collisionless shock acceleration.

Collisionless shock acceleration also shows a clear dominance when the total number of accelerated protons and the beam divergence are considered. The highly collimated acceleration of particles along the shock propagation direction, due to the self-generated collimating magnetic field, guarantees a high quality of the beam. Our results thus indicate the fully quantitative path to laser-driven acceleration in order to achieve high energy and high quality beams.

Additional Information

How to cite this article: Novo, A. S. et al. Optimizing laser-driven proton acceleration from overdense targets. Sci. Rep. 6, 29402; doi: 10.1038/srep29402 (2016).