1. Introduction
In the last two decades, there has been a significant effort in exploring the generation of high-energy (1–100 MeV) ion beams in plasmas produced by intense ($I>10^{18}$ W cm$^{-2}$) laser–solid interactions (Daido, Nishiuchi & Pirozhkov Reference Daido, Nishiuchi and Pirozhkov2012; Macchi, Borghesi & Passoni Reference Macchi, Borghesi and Passoni2013). It has been shown that accelerating gradients as high as teravolts per metre can be created in the plasma (Wilks et al. Reference Wilks, Langdon, Cowan, Roth, Singh, Hatchett, Key, Pennington, MacKinnon and Snavely2001) and ion beams can be produced with small emittance ($<0.1 {\rm \pi}$ mm mrad (Borghesi et al. Reference Borghesi, Mackinnon, Campbell, Hicks, Kar, Patel, Price, Romagnani, Schiavi and Willi2004; Cowan et al. Reference Cowan2004)) and short bunch duration (${\lesssim }1$ ps (Dromey et al. Reference Dromey2016)). These promising results and the potential to produce high-energy, high-current ion beams in more compact systems than solid-state-based accelerators (Cahill et al. Reference Cahill, Rosenzweig, Dolgashev, Tantawi and Weathersby2018) make the study of laser-driven ion acceleration an active area of research.
Laser-driven ion beams are now routinely used for radiography of high-energy-density plasmas (Borghesi et al. Reference Borghesi2002; Rygg et al. Reference Rygg2008) and hold promise for applications in isochoric heating of materials (Patel et al. Reference Patel, Mackinnon, Key, Cowan, Foord, Allen, Price, Ruhl, Springer and Stephens2003; Tahir et al. Reference Tahir2005), fast ignition of inertial confinement fusion (Roth et al. Reference Roth2001), injectors for conventional accelerators (Antici et al. Reference Antici, Fazi, Lombardi, Migliorati, Palumbo, Audebert and Fuchs2008; Aymar et al. Reference Aymar2020) and tumour therapy (Bulanov et al. Reference Bulanov2008; Kraft et al. Reference Kraft2010; Loeffler & Durante Reference Loeffler and Durante2013; Bulanov et al. Reference Bulanov, Wilkens, Esirkepov, Korn, Kraft, Kraft, Molls and Khoroshkov2014; Linz & Alonso Reference Linz and Alonso2016; Kroll et al. Reference Kroll2022). Important requirements for many of these applications are the ability to produce controllable, quasimonoenergetic (low energy spread) ion beams and at a high-repetition rate. These remain significant challenges for laser-driven ion beams, despite the progress in exploring different acceleration mechanisms, and significant developments in high-repetition rate targets (Kim, Göde & Glenzer Reference Kim, Göde and Glenzer2016; Gauthier et al. Reference Gauthier2017; Göde et al. Reference Göde2017; Obst et al. Reference Obst2017; Curry et al. Reference Curry, Schoenwaelder, Goede, Kim, Rehwald, Treffert, Zeil, Glenzer and Gauthier2020).
The most studied laser-driven ion acceleration mechanism to date – both theoretically and experimentally – is the target normal sheath acceleration (TNSA) (Snavely et al. Reference Snavely2000; Wilks et al. Reference Wilks, Langdon, Cowan, Roth, Singh, Hatchett, Key, Pennington, MacKinnon and Snavely2001; Mora Reference Mora2003). Hot electrons are produced near the front surface of the target during the laser–plasma interaction, typically via $\boldsymbol {J}\times \boldsymbol {B}$ (Kruer & Estabrook Reference Kruer and Estabrook1985; May et al. Reference May, Tonge, Fiuza, Fonseca, Silva, Ren and Mori2011) and Brunel (vacuum) heating (Brunel Reference Brunel1987) mechanisms, which are maximized for linearly polarized lasers (Wilks & Kruer Reference Wilks and Kruer1997; Gibbon Reference Gibbon2005). These hot electrons cross the dense target and escape into the vacuum on the rear side, setting up a strong space-charge sheath field that accelerates the target ions from the back surface in the target-normal direction. Ion beams produced by this mechanism are laminar and possess small emittance, however, their energy spectrum is very broad, being typically characterized by an exponentially decreasing energy distribution (Snavely et al. Reference Snavely2000; Wilks et al. Reference Wilks, Langdon, Cowan, Roth, Singh, Hatchett, Key, Pennington, MacKinnon and Snavely2001; Mora Reference Mora2003).
Alternative ion acceleration schemes have been proposed in order to produce more narrow (quasimonoenergetic) ion spectra, including collisionless shock acceleration (CSA) (Denavit Reference Denavit1992; Silva et al. Reference Silva, Marti, Davies, Fonseca, Ren, Tsung and Mori2004; Fiuza et al. Reference Fiuza, Stockem, Boella, Fonseca, Silva, Haberberger, Tochitsky, Gong, Mori and Joshi2012; Haberberger et al. Reference Haberberger, Tochitsky, Fiuza, Gong, Fonseca, Silva, Mori and Joshi2012) and radiation pressure acceleration (RPA) (Wilks et al. Reference Wilks, Kruer, Tabak and Langdon1992; Esirkepov et al. Reference Esirkepov, Borghesi, Bulanov, Mourou and Tajima2004; Macchi et al. Reference Macchi, Cattani, Liseykina and Cornolti2005; Robinson et al. Reference Robinson, Zepf, Kar, Evans and Bellei2008). Collisionless shock acceleration relies on the reflection of ions off a moving electrostatic shock front produced by the laser–plasma interaction near the front surface, which travels at roughly a constant speed inside the target. A small fraction ($\simeq 1\,\%\unicode{x2013}10\,\%$) of the bulk ions is reflected by the shock producing a narrow energy spread ion beam. Hot electrons are important for driving the ion-acoustic waves that mediate shock formation inside the target, and therefore TNSA will also naturally accompany CSA, which can inadvertently broaden the ion spectrum. Specific shaping of the target density has been proposed as a way to control TNSA and produce high-quality beams from CSA (Fiuza et al. Reference Fiuza, Stockem, Boella, Fonseca, Silva, Haberberger, Tochitsky, Gong, Mori and Joshi2012, Reference Fiuza, Stockem, Boella, Fonseca, Silva, Haberberger, Tochitsky, Mori and Joshi2013). Recently, tuning of the plasma density profile using a second laser was shown to produce narrow energy spread ion beams from CSA (Pak et al. Reference Pak2018). However, achieving precise control of the plasma density profile remains a challenge.
Radiation pressure acceleration relies on the slowly varying, cycle-averaged component of the ponderomotive force exerted by the intense laser pulse on the electrons at the front surface of the solid target. It is this radiation pressure that creates a charge separation between the electrons and the ions that accelerates the latter. It can potentially result in the generation of a quasimonoenergetic and laminar ion beam, provided that the accelerating structure is maintained stable and with uniform velocity, similarly to the case of CSA. In theory, RPA can produce ion beams with very high density because almost all ions in the laser focal region can be accelerated by the space-charge field. However, the experimental characterization of this acceleration scheme and observation of narrow energy spread ion beams have been challenging (Henig et al. Reference Henig2009; Palmer et al. Reference Palmer, Dover, Pogorelsky, Babzien, Dudnikova, Ispiriyan, Polyanskiy, Schreiber, Shkolnikov, Yakimenko and Najmudin2011; Kar et al. Reference Kar2012; Steinke et al. Reference Steinke, Hilz, Schnürer, Priebe, Bränzel, Abicht, Kiefer, Kreuzer, Ostermayr, Schreiber, Andreev, Yu, Pukhov and Sandner2013; Bin et al. Reference Bin2015; Scullion et al. Reference Scullion2017; McIlvenny et al. Reference McIlvenny2021). An important difficulty relies on the requirement of low electron heating for efficient momentum transfer from the laser to the ions, and to avoid other competing ion acceleration mechanisms, such as TNSA and CSA, to develop and dominate. Indeed, recent experiments producing nearly 100 MeV proton beams likely involved the combination of different acceleration schemes and the observed energy spectra were broad (Kim et al. Reference Kim2016; Wagner et al. Reference Wagner2016; Higginson et al. Reference Higginson2018; Shen et al. Reference Shen, Qiao, Pukhov, Kar, Zhu, Borghesi and He2021). Other significant challenges include the mitigation of corrugations at the laser–target interaction surface arising due to instabilities (Palmer et al. Reference Palmer2012; Eliasson Reference Eliasson2015; Sgattoni et al. Reference Sgattoni, Sinigardi, Fedeli, Pegoraro and Macchi2015; Wan et al. Reference Wan, Andriyash, Lu, Mori and Malka2020; Chou et al. Reference Chou, Grassi, Glenzer and Fiuza2022) and finite laser spot effects (Klimo et al. Reference Klimo, Psikal, Limpouch and Tikhonchuk2008; Dollar et al. Reference Dollar2012) and the control of the preplasma level that is naturally formed from the preheating and expansion of the target by a laser prepulse, which poses significant constraints on the laser contrast (Varmazyar, Mirzanejhad & Mohsenpour Reference Varmazyar, Mirzanejhad and Mohsenpour2018).
In all these laser-driven ion acceleration schemes electron heating plays a major role in controlling the dominant acceleration mechanism and the quality of the accelerated ion beams. While it is well established that the laser polarization – linear versus circular – can be important in controlling electron heating via the $\boldsymbol {J}\times \boldsymbol {B}$ mechanism (e.g. May et al. Reference May, Tonge, Fiuza, Fonseca, Silva, Ren and Mori2011), it is not clear how, more generally, the different laser properties affect the interplay between competing processes and instabilities at the front surface of the target in order to ensure a robust control of electron heating and ion acceleration.
Here, we use two-dimensional (2-D) and three-dimensional (3-D) particle-in-cell (PIC) simulations with the fully relativistic electromagnetic code OSIRIS (Fonseca et al. Reference Fonseca, Silva, Tsung, Decyk, Lu, Ren, Mori, Deng, Lee, Katsouleas and Adam2002, Reference Fonseca, Martins, Silva, Tonge, Tsung and Mori2008, Reference Fonseca, Vieira, Fiuza, Davidson, Tsung, Mori and Silva2013) to investigate in detail how the laser and plasma properties determine electron heating and how this will impact the quality of the accelerated ion beams. We focus in particular on RPA, exploring both hole boring (HB) and light sail (LS) regimes. We identify the dominant processes and establish a set of criteria relating the laser and target parameters that enable robust mitigation of electron heating and acceleration of high quality ion beams. This work presents a more detailed analysis of recently published results on the optimization of LS ion acceleration (Chou et al. Reference Chou, Grassi, Glenzer and Fiuza2022) and expands on it by presenting new results on the HB regime. The new understanding and set of conditions provided can have an important impact in the guiding of future experiments and in ensuring a better characterization of different ion acceleration regimes, by isolating the dominant mechanisms.
This paper is organized as follows. The physical regimes considered and the simulation set-up used in this study are described in § 2. In § 3, we show that suppression, or significant mitigation of electron heating is required to produce high-quality ion beams based on HB or LS. When this is not achieved, for thick targets a collisionless shock is formed which, in combination with TNSA, supersedes HB, and, for thin targets LS gives way to rapid decompression and transparency of the target. In § 4 we discuss how the deformation of the target surface due to the growth of instabilities and finite spot size effects controls electron heating. We show that the Rayleigh–Taylor-like instability (RTI) is dominant in determining the onset of strong electron heating and controlling the quality of the accelerated ion beams. Based on this understanding, we define a threshold condition for the duration of the laser pulse that allows mitigation of electron heating and high-quality ion acceleration in both HB and LS regimes. The effect of the laser temporal profile on ion acceleration is studied in § 5, where it is shown that the Gaussian temporal profile leads to an increase of the ion energy spread. In § 6 we discuss the importance of controlling laser prepulse to limit the electron heating due to the formation of preplasma. In § 7 we demonstrate that the new pulse duration conditions indeed minimize the electron heating and the energy spread of the accelerated ion beams over a wide range of laser and target parameters. Furthermore, we show that when the new threshold for pulse duration is combined with the condition for the optimal target thickness for LS, it limits the maximum laser intensity that can be used and the maximum peak energy of the accelerated ion beam. In § 8 we verify these findings with 3-D PIC simulations and demonstrate that based on the established criteria it is possible to generate high-quality ion beams from RPA (in both HB and LS regimes) using realistic laser and target conditions. Finally, in § 9 we present our conclusions and discuss the implications of the work for experimental studies.
2. Radiation pressure acceleration regimes and simulation set-up
We consider an intense laser interacting with a planar target with density $n_0>n_c$ (i.e. an overdense target), where $n_0$ and $n_c=m_e\omega _0^2/4{\rm \pi} \, {\rm e}^2$ are the initial plasma density and the critical density associated with laser propagation in the plasma, with $e$ being the elementary charge, $m_e$ the electron mass, and $\omega _0$ the laser frequency. The intense laser exerts a coherent ponderomotive force on the electrons at the surface of the solid target, creating a charge separation between the pushed electrons and the heavier ions, which in turn accelerates the ions. In practice, the laser radiation pressure acts as a piston pushing the plasma at the front surface inwards as the laser light is reflected from the surface. Using a one-dimensional (1-D) model based on energy- and momentum-flux conservation at the target surface, one finds that the front surface is pushed at the known HB velocity (Wilks et al. Reference Wilks, Kruer, Tabak and Langdon1992)
where $P_L = (1+R)I\cos ^2\theta _0/c$ is the radiation pressure exerted by the laser on the target surface in the normal direction, $Z$ and $A$ are the ion charge and mass numbers, $m_i$ and $m_p$ are the ion and proton masses, $c$ the speed of light, $a_0 \simeq 0.85\sqrt {I [\text {W cm}^{-2}](\lambda _0[\mathrm {\mu }\text {m}])^2/10^{18}}$ the peak normalized vector potential, $\theta _0$ the incidence angle, $\lambda _0$ the wavelength of the laser and $R\leqslant$ 1 the reflection coefficient. We note that $R$ may be a function of $\theta _0$, $\lambda _0$, $a_0$ and target density and composition.
We distinguish between two different RPA regimes: HB (Wilks et al. Reference Wilks, Kruer, Tabak and Langdon1992; Macchi et al. Reference Macchi, Cattani, Liseykina and Cornolti2005) and LS (Esirkepov et al. Reference Esirkepov, Borghesi, Bulanov, Mourou and Tajima2004; Robinson et al. Reference Robinson, Zepf, Kar, Evans and Bellei2008; Macchi & Benedetti Reference Macchi and Benedetti2010). The acceleration regime is determined by the ratio $l_0/(v_{\rm HB} \tau _0)$, where $l_0$ is the target thickness and $\tau _0$ is the laser pulse duration. In the HB regime, $l_0 > v_\text {HB} \tau _0$ and the laser radiation pressure can only push a small fraction of the target. Ions are reflected once off the laser piston acquiring a velocity $v_i \simeq 2v_\text {HB}$, or equivalently a peak energy per nucleon $\epsilon _0 = 2m_pv_\text {HB}^2$. In the LS regime, $l_0 < v_\text {HB} \tau _0$ and the laser can push the whole target repeatedly during the laser pulse duration. In other words, the target is accelerated via multiple HB stages (Macchi & Benedetti Reference Macchi and Benedetti2010; Grech et al. Reference Grech, Skupin, Diaw, Schlegel and Tikhonchuk2011). In this case, the target experiences an acceleration (Macchi & Benedetti Reference Macchi and Benedetti2010)
where $\beta _i=v_i/c$ and $\gamma _i$ are the normalized ion velocity and Lorentz factor. The last equality is the approximation in the non-relativistic limit. For negligible laser electron heating ($R\simeq 1$), an exact solution for the final ion velocity $\beta _{i,0}$ and the corresponding energy per nucleon $\epsilon _0$ exists (Macchi, Veghini & Pegoraro Reference Macchi, Veghini and Pegoraro2009)
where $\xi =c({Zm_en_c}/{Am_pn_0})({a_0^2\tau _0}/{l_0})$.
In our simulations a laser pulse with frequency $\omega _0$ is launched along the $x_1$ direction (which is also the direction of the target normal) from the left boundary and irradiates, unless otherwise stated, an electron–proton plasma (i.e. $m_i=m_p= 1836 m_e$) with initial density $n_0\geqslant$ 40 $n_c$. A minimum density of 40 $n_c$ corresponds to that of high-repetition rate liquid hydrogen targets (Kim, Göde & Glenzer Reference Kim, Göde and Glenzer2016; Gauthier et al. Reference Gauthier2017; Curry et al. Reference Curry, Schoenwaelder, Goede, Kim, Rehwald, Treffert, Zeil, Glenzer and Gauthier2020). The plasma density follows a step-like profile with thickness $l_0$. An initial electron temperature $T_e=100$ eV is used (we have checked that our results are not sensitive to the initial temperature choice in the $10\unicode{x2013}1000$ eV range). The typical size of the simulation box in 2-D (3-D) simulations is 400 (300)$c/\omega _0$ longitudinally, and $250 c/\omega _0$ transversely in $x_2$ (and in $x_3$). The 2-D (3-D) simulations use 16 (8) particles per cell per species and a spatial resolution of 0.2 $(0.5) c/\omega _{pe}$ in each direction, where $\omega _{pe}=\sqrt {4{\rm \pi} \, {\rm e}^2n_0/m_e}$ is the electron plasma frequency. The time step is chosen according to the Courant–Friedrichs–Lewy condition. Open (absorbing) boundary conditions for both particles and fields are used in the longitudinal and transverse directions (except in the cases with a plane wave laser where the transverse boundary conditions are periodic). We have tested different resolutions and numbers of particles per cell to ensure convergence of the results and have used a third-order particle interpolation scheme for improved numerical accuracy. We have also tested different domain sizes to ensure that this domain allows capturing the electron heating and ion acceleration dynamics without the build-up of fields at the boundaries due to the absorption of current from escaping particles. Because in most practical applications the primary interest is in highly directional ion beams, for simulations with a finite laser spot the ion energy spectra are integrated within a $10^\circ$ opening angle from the laser propagation (forward) direction. We have checked that this is consistent with selecting the ions within the area of the focal spot.
To study in detail how the laser–plasma parameters affect electron heating and ion acceleration, we have performed a parameter scan in laser intensity ($a_0=5\unicode{x2013}200$), polarization (P-, S- and circular), incidence angle ($\theta _0 = 0^\circ \unicode{x2013}45^\circ$), full-width half-maximum (FWHM) duration ($\tau _0 = 30\unicode{x2013}2000 \omega _0^{-1}$), focal spot (at $1/e^2$ beam width; $w_0=4\unicode{x2013}50 c/\omega _0$ and plane wave), target composition $1\leqslant A/Z\leqslant 4$ for single-species ions, and CH, density ($n_0 = 40\unicode{x2013}500 n_c$, covering the range from liquid hydrogen to solid-density targets) and thickness ($l_0 = 0.08\unicode{x2013}40 c/\omega _0$).
3. Termination of radiation pressure acceleration due to strong electron heating
In this section, we show that for thick targets (HB regime) when there is significant electron heating at the target surface, HB gives rise to the formation of a collisionless shock that is launched into the target and dominates ion acceleration. We illustrate these results with 2-D simulations where an intense ($a_0 = 12$) laser interacts with an overdense thick target ($n_0 = 42 n_c$ and $l_0 = 12 c/\omega _0$). The laser is either a PP (with the electric field along the $x_2$ direction) or CP plane wave, with a fourth-order super-Gaussian temporal profile with $\tau _0 = 100 \omega _0^{-1}$.
In figure 1(a,b), we show the longitudinal phase spaces of electrons and protons for the PP simulation near the end of the laser–plasma interaction at $t = 150 \omega _0^{-1}$. We can clearly see hot electron bunches produced at a frequency of $2\omega _0$, which is a signature of the $\boldsymbol {J}\times \boldsymbol {B}$ heating mechanism (May et al. Reference May, Tonge, Fiuza, Fonseca, Silva, Ren and Mori2011). As a result, a significant fraction of the laser energy goes into the electron population, weakening HB. Indeed, the measured HB velocity is $v_\text {HB}\simeq 0.026 c$, consistent with a laser absorption into hot electrons of $\simeq 50\,\%$ ($R \simeq 0.5$ in (2.1)). In addition, we observe that a collisionless shock forms and dominates proton acceleration as illustrated in figure 1(a,c). The shock front detaches from the surface and propagates into the target at a nearly constant velocity $v_\mathrm {sh} \simeq 0.035 c$. The proton population at $11 c/\omega _0 \lesssim x_1\lesssim 14 c/\omega _0$ has been accelerated by the shock front to $v_i\simeq 2v_\text {sh}$. It is important to note that even after the laser interaction finishes, the collisionless shock continues to propagate through the target and reflect protons (figure 1c). It is also worth noting that due to the generation of hot electrons, a strong space-charge field develops at the rear surface of the target and leads to TNSA, which is evidenced by the proton phase space at $x_1\gtrsim 19 c/\omega _0$. At later times, the protons accelerated by the shock acquire a large energy spread when they leave the target rear surface and experience TNSA. Tailoring of the rear side density profile is required to control TNSA and enable quasimonoenergetic ion beams from CSA (Fiuza et al. Reference Fiuza, Stockem, Boella, Fonseca, Silva, Haberberger, Tochitsky, Gong, Mori and Joshi2012, Reference Fiuza, Stockem, Boella, Fonseca, Silva, Haberberger, Tochitsky, Mori and Joshi2013). In general, we have found that in configurations in which electrons become relativistic, CSA and TNSA will dominate the ion acceleration mechanisms over HB. This highlights the need to prevent or significantly mitigate electron heating in order to enable HB to be the dominant mechanism and to produce high-quality ion beams.
A commonly employed strategy to mitigate electron heating by both $\boldsymbol {J}\times \boldsymbol {B}$ and Brunel mechanisms in laser–plasma interactions is the use of a CP laser at near normal incidence ($\theta _0 \simeq 0^\circ$; (Macchi et al. Reference Macchi, Cattani, Liseykina and Cornolti2005)). The $\boldsymbol {J}\times \boldsymbol {B}$ mechanism relies on the standing wave created by the incoming and reflected laser field. For linear polarization, the oscillation of the magnetic field at the surface allows for electrons to escape the target and experience the electric field of the laser, being accelerated transversely and then rotated back into the target by the magnetic field. For circular polarization, the magnetic field of the standing wave at the target surface does not decrease to zero – it just rotates – and thus electrons cannot escape the target to be efficiently accelerated (May et al. Reference May, Tonge, Fiuza, Fonseca, Silva, Ren and Mori2011). The Brunel heating mechanism relies on a laser electric field component normal to the surface to directly accelerate the electrons. This is absent for normal incidence, provided that the target surface remains uniform and stable. As we will discuss in more detail in the next section, these conditions can only be maintained for very short interaction times. Figure 1(d–f) shows the results of a simulation with the same laser and plasma parameters as in figure 1(a–c) but using CP. We observe that CP is indeed capable of maintaining reduced electron heating, allowing HB to be the dominant ion acceleration mechanism. We observe that the HB velocity is $v_\text {HB} \simeq 0.031 c$, in good agreement with (2.1) for $R \simeq 1$ and that the accelerated protons have $v_1\simeq 2v_\text {HB}$. We also find that, in contrast to the case of a collisionless shock, the HB velocity abruptly slows down when the laser–plasma interaction ends and proton reflection/acceleration ceases (figure 1f at $t\simeq 165 \omega _0^{-1}$). This is an interesting difference between HB and CSA that impacts the total charge accelerated by each mechanism: while HB tends to reflect a larger fraction of the background ions, it can only do so during a shorter period when compared with CSA.
In addition to the use of a CP laser at near normal incidence ($\theta _0 \simeq 0^\circ$), it has been recently proposed in the context of collisionless shock studies (Grassi et al. Reference Grassi, Grech, Amiranoff, Macchi and Riconda2017) that S-polarization (SP) (with the electric field along the $x_3$ direction) with a large incidence angle can also achieve similar results in terms of mitigating electron heating. This is because Brunel heating is absent for SP and $\boldsymbol {J}\times \boldsymbol {B}$ heating can be significantly reduced since it scales with $\cos \theta _0$. We have explored this possibility over a large range of laser intensities ($a_0 = 5\unicode{x2013}30$) and plasma densities ($n_0 = 40\unicode{x2013}150 n_c$) and confirmed that in 2-D simulations indeed SP with $\theta _0=45^\circ$ can significantly suppress electron heating and lead to HB ion acceleration comparable to the CP case (not shown here). However, the situation changes significantly for more realistic 3-D simulations. In 3-D simulations, the laser–plasma interaction along the direction of laser polarization is effectively PP at $\theta _0 \simeq 0^\circ$ and gives rise to significant electron heating along the laser polarization. In all 3-D SP cases tested, hot electrons acquire relativistic temperature and we observe a transition from HB to CSA similar to the PP case illustrated in figure 1(a–f).
For very thin targets, corresponding to the LS regime, we have observed similar results in terms of the mitigation of electron heating, which was achieved with a CP laser with $\theta _0 \simeq 0^\circ$ and led to the acceleration of ions with a narrow energy spread. We note that for PP and SP, the onset of strong electron heating does not lead to the formation of a collisionless shock. Instead, we observe that the thin target quickly becomes transparent to the laser and the radiation pressure is no longer efficient in accelerating ions. Overall, for both HB and LS regimes, we find that the use of CP is required to significantly mitigate electron heating and optimize RPA. Non-normal incidence angles can be useful for SP, but unfortunately seem to be only effective in 2-D simulations. In the remainder of this paper, we discuss in detail the impact of the laser and plasma parameters on electron heating with CP for both HB and LS regimes in order to understand the set of conditions required for high-quality ion acceleration.
4. Development of surface corrugations
The interaction of an intense laser with an overdense target can significantly modify the shape of the target surface either due to the development of surface instabilities or from finite laser spot size effects. These modulate the surface density profile, which can ultimately trigger strong electron heating (e.g. Klimo et al. Reference Klimo, Psikal, Limpouch and Tikhonchuk2008; Dollar et al. Reference Dollar2012; Paradkar & Krishnagopal Reference Paradkar and Krishnagopal2016), even for a CP laser, and impact the quality and mechanisms of ion acceleration as discussed above. In this section, we discuss the importance of surface corrugations on HB and LS acceleration regimes with a CP laser.
4.1. Growth of surface instabilities
Previous studies have investigated the development of density ripples at the interaction surface and there has been significant discussion on which instabilities are dominant, including the Weibel instability (Sentoku et al. Reference Sentoku, Mima, Kojima and Ruhl2000), RTI (Gamaly Reference Gamaly1993; Pegoraro & Bulanov Reference Pegoraro and Bulanov2007; Palmer et al. Reference Palmer2012; Khudik et al. Reference Khudik, Yi, Siemon and Shvets2014; Eliasson Reference Eliasson2015; Sgattoni et al. Reference Sgattoni, Sinigardi, Fedeli, Pegoraro and Macchi2015), electron–ion coupling instabilities (Wan et al. Reference Wan2016; Wan et al. Reference Wan, Pai, Zhang, Li, Wu, Hua, Lu, Joshi, Mori and Malka2018) or a combination of these (Wan et al. Reference Wan, Andriyash, Lu, Mori and Malka2020). However, the correlation between these instabilities and the onset of electron heating has not been studied systematically for both HB and LS. Here, we begin by illustrating the surface dynamics and electron heating for the interaction of a CP laser with Gaussian longitudinal and transverse intensity profiles with a target at normal incidence using 2-D PIC simulations. In the HB regime, we use a laser with $a_0=27$, $\tau _0=200 \omega _0^{-1}$, $w_0=50 c/\omega _0$ with $n_0 = 40\,n_c$ and $l_0=75 c/\omega _0$; for LS, $a_0=15$, $\tau _0=105 \omega _0^{-1}$, $w_0=50 c/\omega _0$, with $n_0 = 250\,n_c$ and $l_0 =0.085 c/\omega _0$ are used. The results are shown in figure 2(a) for HB and figure 3(a) for LS, respectively.
We observe that indeed, even with CP, there is the onset of strong electron heating after a relatively short interaction time. In the HB regime (figure 2a), the electron and ion phase spaces at $t=150 \omega _0^{-1}$ show that electrons have been heated to relativistic temperatures, which leads to a significant weakening of the laser radiation pressure and to the formation of a collisionless shock and development of strong TNSA field that dominate ion acceleration, similarly to what was observed for the PP case discussed in § 3. A sharp increase of the FWHM energy spread of the accelerated ions ensues, reaching ${\rm \Delta} \epsilon /\epsilon _0 > 100\,\%$ as can be seen in figure 2(c); and by $t = 225 \omega _0^{-1}$ both the electron temperature and ion energy spread have saturated at large values.
In the LS regime (figure 3a), we observe that once the electrons are significantly heated they drive the rapid expansion of the target and broadening of the ion energy spread (figure 3c) due to the associated strong space-charge field and short target thickness. By $t = 115 \omega _0^{-1}$ the target expansion leads to the onset of relativistic transparency as can be seen in figure 3(c) from the evolution of $n_e/(\gamma n_c)$ (the ratio of the electron density, $n_e$, and relativistic critical density at the target front surface, where $\gamma$ is the average Lorentz factor of the electrons). At this point, RPA is terminated and the ion energy distribution ceases to be peaked.
From the temporal evolution of the electron temperature $T_e$ (average kinetic energy of electrons) and FWHM ion energy spread ${\rm \Delta} \epsilon /\epsilon _0$ shown in figures 2(c) and 3(c) for the HB and LS regimes, respectively, we observe that they follow a similar behaviour with the ion beam energy spread increasing sharply following the rapid growth of $T_e$ in both cases. We define this time associated with the onset of strong electron heating, $\tau _\text {heating}$, as the time for which the rate of increase of the electron temperature, ${\rm d} T_e/{\rm d} t$, is maximum.
We have repeated the same simulations in both regimes but using a shorter laser pulse with $\tau _0 = 75 \omega _0^{-1}$ and $\tau _0 = 40 \omega _0^{-1}$, both $< \tau _\text {heating}$, for HB and LS, respectively, to confirm the impact of electron heating on the growth of the ion energy spread and overall target dynamics. The results are shown in figures 2(b,c) and 3(b,c). In these cases, we observe that indeed both $T_e$ and ${\rm \Delta} \epsilon /\epsilon _0$ remain low (figures 2c and 3c). The ion beam energy spread saturates at $t \simeq 2 \tau _0$ with ${\rm \Delta} \epsilon /\epsilon _0 \ll 1$ and remains stable long after the laser–plasma interaction has finished.
We have found that the onset of electron heating in both HB and LS regimes is related to the emergence of large transverse density modulations at the target front surface with a wavelength comparable to that of the laser. For the longer pulse simulations ($\tau _0 > \tau _\text {heating}$), these surface density modulations are visible in the third rows of figures 2(a) and 3(a) at the times where electron heating is also observed. Density modulations with a wavelength comparable to $\lambda _0$ allow the penetration of the laser in the lower density regions and resonant enhancement of its electric field (Eliasson Reference Eliasson2015; Sgattoni et al. Reference Sgattoni, Sinigardi, Fedeli, Pegoraro and Macchi2015), giving rise to effective electron heating, for example via the Brunel mechanism (bottom rows of figures 2a and 3a), with the temperature reached being comparable to that observed in simulations with a linearly polarized laser (not shown here). The spatial distribution of $T_e$ shows indeed that the heating is happening at the walls of these concave valleys and consistent with direct acceleration by the laser electric field. The location of the hot spots of $T_e$ oscillates (from top to bottom of the valleys) in accordance with the phase of the laser electric field. In the simulations with $\tau _0 < \tau _\text {heating}$, the amplitudes of the density modulations at the surface are much smaller during the time of laser interaction (third rows of figures 2b and 3b) leading to a much-reduced electron heating and stable ion acceleration.
In order to study the mechanism responsible for these corrugations, and isolate the effects of surface instabilities, we have performed a parameter scan of 2-D simulations with a long, plane-wave CP laser at normal incidence, with parameters varied in the following ranges: for HB, $5 \leqslant a_0 \leqslant 60, 40 n_c \leqslant n_0 \leqslant 200 n_c$; for LS, $5 \leqslant a_0 \leqslant 200, 40 n_c \leqslant n_0 \leqslant 500 n_c$ and $0.08 c/\omega _0 \leqslant l_0 \leqslant 2 c/\omega _0$. Note that for the LS regime, the initial target thickness $l_0$ is always larger than or equal to the optimal LS target thickness $l_\text {opt}$, defined as $l_\text {opt}= a_0\lambda _0n_c/(\sqrt {2}{\rm \pi} n_0)$, for which the acceleration is maximized by minimizing the total target mass while guaranteeing that the target remains relativistically opaque (Macchi et al. Reference Macchi, Veghini and Pegoraro2009). For both regimes, the target is either composed of single-species ions with $1\leqslant A/Z \leqslant 4$, or CH (plastic). We analyse the growth of ion density modulations by computing the transverse Fourier modes ($k_{x_2}$) of the longitudinal ($x_1$) displacement of the relativistic critical surface ($n \simeq \gamma _0n_c$; where $\gamma _0=\sqrt {1+a_0^2/2}$ is the electron Lorentz factor), as a function of the transverse ($x_2$) position (e.g. figure 4b,e). Previous theoretical studies of surface instabilities considered perturbations of the surface displacement and showed it will grow exponentially due to instability (e.g. Gamaly Reference Gamaly1993; Eliasson Reference Eliasson2015). Alternatively, for LS, one could obtain the Fourier modes of the amplitude of the transverse density profile by integrating over the target longitudinally. We have checked that the obtained growth rates and time scales of the different modes are consistent between both methods.
Figure 4 illustrates the growth of different modes for simulations with the same laser intensity and target parameters of figures 2 and 3. Note that both targets remain opaque to the laser during the time of the analysis and thus the measurements of the growth rates and saturation levels are not affected by the onset of relativistic transparency. The fastest growing modes are observed at $k_{x_2}\gg k_0$, with $k_0=2{\rm \pi} /\lambda _0$ (e.g. $k_{x_2}\simeq 13 \omega _0/c$ in figure 4e at $t\simeq 30 \omega _0^{-1}$), and have been previously described as associated with electron–ion coupling instabilities (Wan et al. Reference Wan2016; Wan et al. Reference Wan, Pai, Zhang, Li, Wu, Hua, Lu, Joshi, Mori and Malka2018, Reference Wan, Andriyash, Lu, Mori and Malka2020). However, these modes saturate at relatively low amplitude and do not lead to significant electron heating. The dominant density modulations are associated with the mode with $k_{x_2}\simeq k_0$ and we observe the onset of strong electron heating during the linear growth and saturation of this mode. This suggests that the onset of strong electron heating is related to laser-driven RTI for which the dominant mode is $k_{x_2}\simeq k_0$ (Eliasson Reference Eliasson2015; Sgattoni et al. Reference Sgattoni, Sinigardi, Fedeli, Pegoraro and Macchi2015).
To further confirm that the RTI mode with $k_{x_2} = k_0$ is the dominant effect on the onset of strong electron heating and degradation of the ion beam quality, we have performed additional simulations with different transverse domains. We observe that for transverse domain sizes $< \lambda _0/2$, where the dominant RTI mode is prohibited, no significant electron heating is observed. In these cases, the fastest growing high-k modes are still captured and thus can still grow, but $T_e$ remains very low (figure 5). For the largest transverse domain size $\gg \lambda _0$, we see that strong electron heating starts near the saturation time of the RTI.
For the laser-driven RTI, the growth rate is $\varGamma _\text {RTI} \propto \sqrt {a_{\rm RPA}k_0}$, where $a_{\rm RPA}=2v_\text {HB}^2/l$ is the acceleration due to radiation pressure ((2.2) with $R\simeq 1$ and $\beta _i\ll 1$). For LS, the target thickness is $l_0$, whereas for HB the effective acceleration layer is characterized by the relativistic electron skin depth $l=\sqrt {a_0}c/\omega _{pe}$ (Gamaly Reference Gamaly1993). Linear fits to the measured growth rates of the $k_0$ mode from the simulations confirm the expected scaling with the following numerical factors (figure 6; see also the example fits in figure 4c,f): $\varGamma _\text {RTI}[\omega _0]\simeq 0.3a_0^{3/4}(n_0[n_c])^{-1/4}( Zm_e/(Am_p))^{1/2}$ for HB and $\varGamma _\text {RTI}[\omega _0]\simeq 0.5a_0(n_0[n_c]l_0[c/\omega _0] Am_p/(Zm_e))^{-1/2}$ for LS.
We find a strong correlation between $\tau _\text {heating}$ and the growth time of the instability, with $\hat {\tau }_{0,\text {RTI}} \equiv \tau _\text {heating} \simeq 5\varGamma ^{-1}_\text {RTI}$ for HB and $\hat {\tau }_{0,\text {RTI}} \equiv \tau _\text {heating} \simeq 3\varGamma ^{-1}_\text {RTI}$ for LS, where $\hat {\tau }_{0,\text {RTI}}$ is defined as the time of the onset of strong electron heating due to the development of the RTI. This is shown in figures 6(b) and 6(d). In order to suppress or significantly mitigate electron heating, the duration of the laser pulse $\tau _0$ should then be smaller than $\hat {\tau }_{0,\text {RTI}}$, which can be written as
It is important to note that although several previous works have studied the development of the RTI, the focus had been on developing strategies to mitigate the penetration of the RTI fingers on the accelerated proton species, which included the use of mixed ion species (Yu et al. Reference Yu, Pukhov, Shvets and Chen2010, Reference Yu, Pukhov, Shvets, Chen, Ratliff, Yi and Khudik2011), advanced laser configurations (Wu et al. Reference Wu, Zheng, Qiao, Zhou, Yan, Yu and He2014; Zhou et al. Reference Zhou, Yan, Mourou, Wheeler, Bin, Schreiber and Tajima2016) and curved targets (Wang, Khudik & Shvets Reference Wang, Khudik and Shvets2021). These are either challenging to implement in practice (Zhou et al. Reference Zhou, Yan, Mourou, Wheeler, Bin, Schreiber and Tajima2016; Wang et al. Reference Wang, Khudik and Shvets2021; and have not yet been proven to be effective experimentally), or still lead to significant electron heating (Yu et al. Reference Yu, Pukhov, Shvets and Chen2010; Wu et al. Reference Wu, Zheng, Qiao, Zhou, Yan, Yu and He2014). As we show here, a quantitative understanding of the detrimental effect that the instability-induced electron heating has on the ion beam quality is critical to produce ion beams with high spectral quality.
We should further note that the results presented here have considered only the regime where ions are non-relativistic, which is appropriate for most current and near-future laser systems. In the relativistic regime, the RTI can still grow as shown in previous numerical studies (Bulanov et al. Reference Bulanov, Echkina, Esirkepov, Inovenkov, Kando, Pegoraro and Korn2010; Sgattoni, Sinigardi & Macchi Reference Sgattoni, Sinigardi and Macchi2014), but how its growth rate changes and, more generally, its impact on electron heating and the spectral quality of the accelerated ions is not well established. One expects that strong electron heating will still be caused by the development of RTI in the relativistic regime – Brunel heating will still be present. The resulting space-charge fields will also lead to an increase of the ion energy spread, however, the rate at which this happens may be more moderate in the relativistic regime when compared with the non-relativistic case. Furthermore, the resulting expansion of the target will also pose limitations on the acceleration due to the onset of relativistic transparency, as in the cases discussed here. A detailed analysis of the relativistic regime is left for future work.
4.2. Finite laser spot size
The transverse variation of the laser intensity due to its spatial profile naturally leads to non-uniform HB velocities across the surface and results in a change of the surface shape over time, also triggering strong electron heating. The onset time of strong electron heating is approximately when the radiation-pressure-driven displacement $d$ of the target surface is comparable to the laser spot size $w_0$, resulting in a significant change in the local incidence angle. This effect has been previously discussed in the case of normal laser incidence ($\theta _0 = 0^\circ$ (Klimo et al. Reference Klimo, Psikal, Limpouch and Tikhonchuk2008; Wan et al. Reference Wan, Andriyash, Lu, Mori and Malka2020)). In fact, this argument can be generalized to other laser incidence angles if we consider the displacement $d$ to be along the axis of the laser resulting in a time for electron heating that is independent of $\theta _0$. We define this time as $\hat {\tau }_{0,\text {FS}}$ (‘FS’ stands for finite spot), and can estimate it as $d(\hat {\tau }_{0,\text {FS}})= w_0$, where for simplicity we assume that the laser intensity is approximately constant in time. For HB, we have $d(t)=v_\text {HB}t$. For LS, we consider first that the target experiences a constant acceleration in its rest frame – i.e. $a_{\rm RPA}\simeq 2v_\text {HB}^2/l_0$ ((2.2) with $R\simeq 1$ and $\beta _i\ll 1$). This then gives $d(t)=(c\sqrt {c^2+a_{\rm RPA}^2t^2}-c^2)/a_{\rm RPA}\simeq (cl_0/2v_\text {HB}^2)(\sqrt {c^2+4v_\text {HB}^4t^2/l_0^2}-c)$ (for all cases considered, we have found that the error introduced by taking the non-relativistic ion velocity limit is $\lesssim 5$ %). Equating $d(\hat {\tau }_{0,\text {FS}})=w_0$ yields
We note that the estimate above considers near diffraction-limited laser focusing where the transverse intensity profile is smooth. If the Strehl ratio is low and speckle-like intensity distributions are present at focus, these will lead to modulations of the surface and laser incidence angle on scales comparable to the speckle size. In that case, in the threshold condition given in (4.2) we should replace $w_0$ by the size of the laser intensity speckles, which can be a significant limitation for small-scale (${\sim }\lambda _0$) speckles. We further note that with oblique incidence ($\theta _0 > 0$), the ion beam direction will be modified. This can happen due to partial absorption of the laser field, which will impart transverse (along the target surface) momentum to the ions (Macchi et al. Reference Macchi, Grassi, Amiranoff and Riconda2019) and also as the surface is modified to be nearly normal to the laser pulse causing the ion beam direction to be primarily along the laser propagation direction.
4.3. Dominance of the RTI
We now compare the constraints on $\tau _0$ imposed by the RTI via (4.1) and the finite spot size effect via (4.2). Rayleigh–Taylor-like instability dominates when $\hat {\tau }_{0,\text {RTI}} < \hat {\tau }_{0,\text {FS}}$, or equivalently $w_0/\lambda _0\gtrsim 2[a_0 n_c/n_0]^{1/4}$ for HB, and $w_0/\lambda _0 \gtrsim 3$ for LS. These conditions are met for targets that are relativistically opaque ($a_0 < n_0/n_c$) and for typically used spot sizes $w_0\gtrsim$ (2–3) $\lambda _0$ (assuming near diffraction-limited laser intensity profiles at focus). Thus, we expect that for conditions of interest for laser-driven ion acceleration in overdense targets the surface corrugations by RTI discussed here impose the main limitation on the pulse duration, (4.1), for high-quality ion beam acceleration.
5. Laser temporal profile
In this section, we evaluate how different laser temporal profiles impact the quality of the ion beam. This can be particularly relevant for the HB regime, as different ion populations will experience different HB velocities due to $v_\text {HB}(t)\propto a_0(t)$. We have performed simulations in both HB and LS regimes, for a plane wave laser with $\tau _0<\hat {\tau }_{0,\text {RTI}}$, where the laser temporal profile is either Gaussian or a flat-top (approximated by a fourth-order super-Gaussian). The total laser energy and $a_0$ is kept the same in both cases. The range of $a_0$, $n_0$ and $l_0$ explored was similar to § 4, and we have varied the pulse duration in the range $0.2 \hat {\tau }_{0,\text {RTI}}\lesssim \tau _0 < \hat {\tau }_{0,\text {RTI}}$.
In general, we find that the accelerated proton beams have a relatively narrow energy spread with both profiles, but the energy spread in the Gaussian case is typically larger, by a factor up to $\simeq 2$. Figure 7 illustrates the typical differences between both profiles for HB and LS. One important feature that we observe in the HB regime is the development of an extended low-energy population (<3 MeV in figure 7a). This is due to the contribution from the low HB velocity phase at the edges of the laser pulse, and can be understood as follows. In a time interval $t$ to $t+{\rm d}t$, a population of ions ${\rm d}N\propto n_iv_\text {HB}(t)\,{\rm d}t$ is accelerated to a velocity $2v_\text {HB}(t)\propto a_0(t)$ (for $R\simeq 1$). Therefore, the resulting energy spectrum will be
for a Gaussian pulse, where $\tilde {\epsilon }_0 = 2 m_i v_\text {HB,0}^2$ is the peak energy of the ions (with $v_\text {HB,0}$ the HB velocity associated with the peak intensity). This matches well the low-energy component of the proton spectrum in figure 7(a), which is responsible for the additional energy spread with respect to the flat-top case.
6. Laser prepulse
In this section, we explore the impact that the development of a preplasma induced by a laser prepulse can have in triggering early electron heating and affecting the spectral quality of ion beams accelerated via RPA. We consider a preplasma with an exponential density profile at the front surface of the target with scale length $l_g$. For the laser and target parameters considered in the short-pulse cases in figure 2(b) and 3(b), with $\tau _0 = 75 \omega _0^{-1}$ and $l_0 = 75 c/\omega _0$ for HB and $\tau _0 = 40 \omega _0^{-1}$ and $l_0 = 0.085 c/\omega _0$ for LS, we have performed additional simulations with different levels of $l_g$, which was varied in the range $0.01 \leqslant l_g/l_0 \leqslant 0.4$ (the total target mass is conserved).
Figure 8 shows that the preplasma can significantly impact electron heating and the ion energy spread for levels of $l_g \gtrsim 0.1\,l_0$. In the case of HB, the electron temperature and ion energy spread for a preplasma scale length $l_g < 0.1 l_0$ ($l_g < 1.2\ \mathrm {\mu }$m for $\lambda _0 = 1\ \mathrm {\mu }$m) are similar to the case of no preplasma, but above this level strong electron heating is triggered, leading to a very fast degradation of ion beam quality even before the pulse arrives at the main target. For LS, for a preplasma with $l_g \leq 0.1 l_0$ the results are also very similar to the no preplasma case. However, for $l_g \gtrsim 0.2 l_0$ ($l_g \gtrsim 2.7$ nm for $\lambda _0 = 1\ \mathrm {\mu }$m) both the electron temperature and ion energy spread are observed to increase to nearly twice the values of the no preplasma case. These results confirm the need to carefully control the level of preplasma to produce ion beams with high spectral quality from RPA and will help inform the laser prepulse contrast requirements for future experimental studies.
7. Optimal regime of radiation pressure acceleration
Our findings make clear the importance of limiting the pulse duration to control electron heating and obtain quasimonoenergetic ion beams from RPA in both HB and LS regimes. In figure 9 we demonstrate that (4.1) is robust over a wide range of laser and target parameters even when realistic Gaussian transverse and temporal pulse profiles are considered. We observe that, indeed, the derived threshold condition marks the transition from low to high electron heating (figure 9a,c) and consequently from low to high energy spread of the accelerated ion beam (figure 9b,d). By repeating some of the simulations in the high-quality regimes using a small, but finite laser incidence angle, we have also confirmed that in general for an incidence angle $\lesssim 10^\circ$, as typically used experimentally, electron heating is still maintained at a low level and the quality of the ion beam remains similar to the case with normal incidence.
For the LS regime, combining the threshold condition for the pulse duration (see (4.1)) and the optimal target thickness condition $l_0=l_\text {opt}$ results in a new condition on the maximum laser intensity with important implications for the optimization of ion beam energy spread in LS acceleration (Chou et al. Reference Chou, Grassi, Glenzer and Fiuza2022). In particular, we obtain that the laser $a_0$ should be limited by
and consequently the maximum energy per nucleon $\epsilon _0$ of the narrow energy spread peak is
In contrast to previous works, and to the common practice of pushing for higher $a_0$ to generate stable LS and higher energy ion beams (e.g. Qiao et al. Reference Qiao, Zepf, Borghesi and Geissler2009), (7.1) and (7.2) indicate the existence of a maximum $\hat {a}_0$ and $\hat {\epsilon }_0$ for high-quality LS ion beams, and show that these are determined primarily by the laser wavelength and pulse duration. In particular, it is worth highlighting that for this optimal regime of acceleration the energy per nucleon of the ion spectral peak does not depend on the target density, composition and laser energy (transverse spot size). These predictions for the LS regime have been recently validated with 3-D simulations in Chou et al. (Reference Chou, Grassi, Glenzer and Fiuza2022).
8. The 3-D simulation results
To further explore RPA in more realistic 3-D configurations and validate our model for the optimal laser duration for high-quality ion beams, we have performed 3-D simulations in both HB and LS regimes.
In the HB regime, a CP laser with $a_0=12$ ($I \simeq 2\times 10^{20} \text {W cm}^{-2}$ for $\lambda _0 = 1\ \mathrm {\mu }$m), $w_0 = 2\ \mathrm {\mu }$m is incident with $\theta _0 = 0^\circ$ on a planar electron–proton target with $n_0 = 40 n_c$ and thickness of $3.5\ \mathrm {\mu }$m, corresponding to the conditions of typical liquid hydrogen jet targets (Kim, Göde & Glenzer Reference Kim, Göde and Glenzer2016; Gauthier et al. Reference Gauthier2017; Göde et al. Reference Göde2017; Obst et al. Reference Obst2017). A fourth-order super-Gaussian temporal profile is used. For these parameters, (4.1) indicates $\tau _0 \lesssim 160$ fs for high-quality proton beams to be produced. We thus run three simulations with pulse durations of $\tau _0=80$, 160 and $265$ fs to test this criterion. In figure 10(a), we show the energy spectra of the proton beams exiting the target from the rear surface for these three cases. When the pulse duration is smaller than the predicted threshold, a quasimonoenergetic proton beam is generated, peaking at roughly $1.8$ MeV, which corresponds to $v_i=2v_\text {HB}\simeq 0.031 c$ and is consistent with (2.1) for $R\simeq 1$. When the pulse duration is comparable to the threshold, we observe that the ion spectral peak at similar energy is still visible, but less prominent. For the longer pulse duration ($\tau _0=265$ fs) we find substantial electron heating ($T_e\simeq 3$ MeV), and a strong TNSA field is generated at the rear surface, which broadens the proton energy spread. No clear spectral peak is observed at the same energy. These results confirm the validity of the derived limit on pulse duration for HB based on the development of surface corrugations and associated electron heating.
We note that the generation of narrow energy spread HB ion beams with similar peak energies around 1 MeV have been reported in Palmer et al. (Reference Palmer, Dover, Pogorelsky, Babzien, Dudnikova, Ispiriyan, Polyanskiy, Schreiber, Shkolnikov, Yakimenko and Najmudin2011), where the laser pulse duration appears to meet the criterion of our model. Specifically, they considered a CP CO$_2$ laser with $\lambda _0 = 10\ \mathrm {\mu }$m, $a_0\simeq 0.7$, $w_0 = 70\ \mathrm {\mu }$m and a hydrogen gas jet target with $n_0\lesssim 10 n_c$, $l_0\sim 800\ \mathrm {\mu }$m. The pulse duration used $\tau _0 \simeq 6$ ps was less than $\hat {\tau }_{0,\text {RTI}}\simeq 9$ ps as required by (4.1). The measured ion beam energies were also shown to be consistent with efficient HB.
For the LS regime, additional 3-D simulations have also been performed to illustrate the change in the ion beam quality for different pulse durations. We simulate a CP laser impinging on a planar electron–proton target with $n_0 = 100 n_c$ and $l_0 = 350$ nm at normal incidence. The laser pulse is temporally Gaussian and the transverse profile is either Gaussian with $w_0 \simeq 7\ \mathrm {\mu }$m or plane wave. These parameters are similar to those used in Qiao et al. (Reference Qiao, Zepf, Borghesi and Geissler2009), where the generation of a GeV proton beam was obtained in 2-D simulations and it was argued that extreme laser intensities were needed to obtain stable high-quality ion beams. Two cases were illustrated: unstable acceleration with $a_0 = 117$ and a stable regime with $a_0 = 214$. Under these conditions, (4.1) predicts that high-quality proton beams require $\tau _0 \lesssim 20$ fs for $a_0 = 117$ and $\tau _0\lesssim 10$ fs for $a_0 = 214$, both of which are much smaller than the duration $\tau _0\simeq 64$ fs used in the cited work.
Figure 10(b) shows the comparison of the proton spectra at the same time as figure 4 in Qiao et al. (Reference Qiao, Zepf, Borghesi and Geissler2009), which corresponds to $t\simeq 2\tau _0$ after the laser interaction. For the case with $a_0=117$, when a pulse duration $\tau _0 = 20$ fs is used, we observe the generation of a high-quality proton beam with $250$ MeV, ${\rm \Delta} \epsilon /\epsilon _0 \simeq 13$ %. For the same parameters but using $\tau _0 = 64$ fs (not shown) we observe the generation of a beam with $460$ MeV, ${\rm \Delta} \epsilon /\epsilon _0 \simeq 32$ %, similar to their green curve. This indicates that it is the onset of electron heating associated with the surface instability in the longer pulse that leads to the increase of the energy spread and to a significant reduction of the coupling efficiency.
For the highest intensity case ($a_0=214$), we observe that similarly the proton beam energy spread is improved from ${\rm \Delta} \epsilon /\epsilon _0 = 20$ % with $\tau _0=64$ fs to ${\rm \Delta} \epsilon /\epsilon _0 = 12$ % with $\tau _0=10$ fs. We further observe that for a transverse Gaussian laser profile the importance of the short pulse duration is even more dramatic. Defining the laser-to-proton energy conversion efficiency $\eta$ as the fraction of the laser energy being carried by the beam within a $10^\circ$ opening angle, for $\tau _0 = 10$ fs we obtain high-quality proton beams with ${\sim }600$ MeV, ${\rm \Delta} \epsilon /\epsilon _0 = 30$ % and $\eta = 27$ %, whereas with $\tau _0 = 64$ fs, we have ${\sim }100$ MeV, ${\rm \Delta} \epsilon /\epsilon _0 = 45$ % and $\eta = 0.5$ %. The proton beam energy obtained with the short pulse ($\tau _0 = 10$ fs) is comparable to the prediction from 1-D theory of 700 MeV (see (2.3a,b)). We note that this is slightly higher than the prediction of (7.2) because the target thickness considered in Qiao et al. (Reference Qiao, Zepf, Borghesi and Geissler2009) was $l_0\lesssim l_\text {opt}$. More importantly, these results confirm that the process that controls the stability and spectral quality of the LS accelerated ion beams is the electron heating via the development of RTI at the target surface and that choosing the appropriate pulse duration and laser intensity is critical for the acceleration of high-quality ion beams.
Figure 10(b) also includes the results of a 3-D PIC simulation for which the laser and target parameters satisfy the optimal LS regime (see (7.1)). We have chosen a set of parameters relevant for near-future short-pulse laser facilities, where a Gaussian laser pulse profile with $\tau _0 = 15$ fs, $w_0 = 5\ \mathrm {\mu }$m and $a_0 = \hat {a}_0 \simeq 122$ according to (7.1) is used. The target has $n_0 = 250 n_c$ and $l_0 = l_\text {opt}$. Under this optimal regime, we indeed observe stable acceleration of the protons via LS leading to the generation of a narrow energy spread proton beam with peak energy $\epsilon _0 \simeq 310$ MeV in very good agreement with the prediction of (7.2), ${\rm \Delta} \epsilon /\epsilon _0 \simeq$ 25 % and $\eta \simeq 2\,\%$. The total laser-to-proton energy conversion efficiency into $4{\rm \pi}$ is $\sim 20\,\%$.
9. Conclusions
In summary, we have presented a detailed study of electron heating and ion acceleration in the interaction of intense laser pulses with overdense plasmas, with a particular focus on radiation pressure acceleration. We have shown that electron heating controls the interplay between different ion acceleration mechanisms. For circularly polarized lasers, the onset of strong electron heating is dominated by the development of corrugations of the target surface on the scale of the laser wavelength due to the laser-driven RTI. In the HB regime, electron heating leads to the development of a collisionless shock and the transition to CSA and TNSA. In the LS regime, electron heating causes fast target expansion and onset of relativistic transparency leading to a significant increase of the ion energy spread. We have shown that to reduce or suppress electron heating and obtain high-quality (low energy spread) ion beams by RPA it is critical to use laser pulses much shorter than the saturation time of the RTI. Using 2-D and 3-D PIC simulations, we have demonstrated that when such short pulses are used high-quality proton beams can be produced with maximum energy comparable to the optimal 1-D RPA theory, even for lasers with Gaussian transverse and temporal intensity profiles.
Interestingly, we note that a few of the previous experimental studies that reported low energy spread ion beams from either the HB (Palmer et al. Reference Palmer, Dover, Pogorelsky, Babzien, Dudnikova, Ispiriyan, Polyanskiy, Schreiber, Shkolnikov, Yakimenko and Najmudin2011) or LS (Henig et al. Reference Henig2009; Steinke et al. Reference Steinke, Hilz, Schnürer, Priebe, Bränzel, Abicht, Kiefer, Kreuzer, Ostermayr, Schreiber, Andreev, Yu, Pukhov and Sandner2013) acceleration appear to meet the laser duration criteria developed in this work. This is encouraging and the understanding provided by our work can help guide future experimental developments in this area. For example, for the parameters of high-power, high-contrast state-of-the-art and near future laser systems with $\tau _0 \simeq 15$ fs, such as the ELI-NP (Doria et al. Reference Doria, Cernaianu, Ghenuche, Stutman, Tanaka, Ticos and Ur2020), Apollon 10 PW (Papadopoulos et al. Reference Papadopoulos2016) and EP-OPAL (Bromage et al. Reference Bromage2019) facilities, our 3-D PIC simulations demonstrate the possibility to produce ${\sim }300$ MeV proton beams with ${\sim }25\,\%$ energy spread and high coupling efficiency.
Acknowledgements
The authors thank the OSIRIS Consortium, consisting of UCLA and IST (Portugal) for the use of the OSIRIS 4.0 framework and the visXD framework. Simulations were performed at Cori (NERSC) and Theta (ALCF) through ERCAP and ALCC computational grants.
Editor L.O. Silva thanks the referees for their advice in evaluating this article.
Declaration of interest
The authors report no conflict of interest.
Funding
This work was supported by the United States Department of Energy SLAC contract no. DEAC02-76SF00515, by the United States DOE Early Career Research Program under FWP 100331 and by the DOE FES under FWP 100182.