Hostname: page-component-7c8c6479df-hgkh8 Total loading time: 0 Render date: 2024-03-19T01:42:49.931Z Has data issue: false hasContentIssue false

Enhancement of laser-driven ion acceleration in non-periodic nanostructured targets

Published online by Cambridge University Press:  17 January 2020

J. Ferri*
Affiliation:
Department of Physics, Chalmers University of Technology, SE-41296Göteborg, Sweden
I. Thiele*
Affiliation:
Department of Physics, Chalmers University of Technology, SE-41296Göteborg, Sweden
E. Siminos
Affiliation:
Department of Physics, Gothenburg University, SE-41296Göteborg, Sweden
L. Gremillet
Affiliation:
CEA, DAM, DIF, F-91297Arpajon, France
E. Smetanina
Affiliation:
Department of Physics, Gothenburg University, SE-41296Göteborg, Sweden
A. Dmitriev
Affiliation:
Department of Physics, Gothenburg University, SE-41296Göteborg, Sweden
G. Cantono
Affiliation:
Department of Physics, Lund University, SE-22100Lund, Sweden
C.-G. Wahlström
Affiliation:
Department of Physics, Lund University, SE-22100Lund, Sweden
T. Fülöp
Affiliation:
Department of Physics, Chalmers University of Technology, SE-41296Göteborg, Sweden
*
Email addresses for correspondence: julien.ferri@chalmers.se, illia-thiele@web.de
Email addresses for correspondence: julien.ferri@chalmers.se, illia-thiele@web.de
Rights & Permissions [Opens in a new window]

Abstract

Using particle-in-cell simulations, we demonstrate an improvement of the target-normal-sheath acceleration (TNSA) of protons in non-periodically nanostructured targets with micron-scale thickness. Compared to standard flat foils, an increase in the proton cutoff energy by up to a factor of two is observed in foils coated with nanocones or perforated with nanoholes. The latter nano-perforated foils yield the highest enhancement, which we show to be robust over a broad range of foil thicknesses and hole diameters. The improvement of TNSA performance results from more efficient hot-electron generation, caused by a more complex laser–electron interaction geometry and increased effective interaction area and duration. We show that TNSA is optimized for a nanohole distribution of relatively low areal density and that is not required to be periodic, thus relaxing the manufacturing constraints.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2020

1 Introduction

Laser-driven ion acceleration has become a well-established technique to produce compact, high-energy ion beams, owing to the ultra-strong accelerating fields that can be achieved at the surfaces of solid targets (Daido, Nishiuchi & Pirozhkov Reference Daido, Nishiuchi and Pirozhkov2012; Macchi, Borghesi & Passoni Reference Macchi, Borghesi and Passoni2013). Such ion sources show great potential for a number of applications ranging from radiography (Romagnani et al. Reference Romagnani, Fuchs, Borghesi, Antici, Audebert, Ceccherini, Cowan, Grismayer, Kar and Macchi2005) to nuclear photonics (Habs et al. Reference Habs, Thirolf, Gross, Allinger, Bin, Henig, Kiefer, Ma and Schreiber2011) and proton therapy (Bulanov & Khoroshkov Reference Bulanov and Khoroshkov2002). However, even though proton energies close to 100 MeV have been demonstrated in recent experiments using petawatt-class laser facilities (Wagner et al. Reference Wagner, Deppert, Brabetz, Fiala, Kleinschmidt, Poth, Schanz, Tebartz, Zielbauer and Roth2016; Higginson et al. Reference Higginson, Gray, King, Dance, Williamson, Butler, Wilson, Capdessus, Armstrong and Green2018), the few tens of MeV energies that are routinely attained using multi-terawatt-class laser systems are insufficient for many of the foreseen applications, therefore limiting the applicability of laser-driven ion sources. This spurs the development of novel schemes yielding significantly increased proton energies.

The most robust, and extensively investigated, acceleration scheme is the so-called target-normal-sheath acceleration (TNSA) (Snavely et al. Reference Snavely, Key, Hatchett, Cowan, Roth, Phillips, Stoyer, Henry, Sangster and Singh2000; Wilks et al. Reference Wilks, Langdon, Cowan, Roth, Singh, Hatchett, Key, Pennington, MacKinnon and Snavely2001), whereby surface ions are driven outwards by the charge-separation field set up by the laser-accelerated relativistic electrons escaping into vacuum. Because of their large charge-to-mass ratio, the protons that are naturally present due to hydrogen-containing contaminants at the target surfaces respond the fastest to the electric sheath field, and reach the highest velocities. Their final energy spectrum has typically the form of a decreasing exponential with a sharp high-energy cutoff.

Different strategies have been explored in recent years to increase the proton cutoff energies resulting from TNSA. With micrometric foil targets, this requires enhancing the fast electron generation at the target front side. One option is to manipulate the laser temporal profile so as to create a preplasma with an optimal scale length (Kaluza et al. Reference Kaluza, Schreiber, Santala, Tsakiris, Eidmann, Meyer-ter Vehn and Witte2004; Nuter et al. Reference Nuter, Gremillet, Combis, Drouin, Lefebvre, Flacco and Malka2008; Brenner et al. Reference Brenner, Robinson, Markey, Scott, Gray, Rosinski, Deppert, Badziak, Batani and Davies2014), or to induce an optimal electromagnetic interference pattern (Ferri, Siminos & Fülöp Reference Ferri, Siminos and Fülöp2019). Another option is to modify the target properties: reduction of the target thickness (Neely et al. Reference Neely, Foster, Robinson, Lindau, Lundh, Persson, Wahlström and McKenna2006; Ogura et al. Reference Ogura, Nishiuchi, Pirozhkov, Tanimoto, Sagisaka, Esirkepov, Kando, Shizuma, Hayakawa and Kiriyama2012) or transverse size (Buffechoux et al. Reference Buffechoux, Psikal, Nakatsutsumi, Romagnani, Andreev, Zeil, Amin, Antici, Burris-Mog and Compant-La-Fontaine2010) thus results in higher proton energies and numbers. An alternative, which is addressed in the present paper, is to employ nano- (or micro-) structured targets. Enhancement of the proton energy has been discussed for periodic nanohole targets (Nodera et al. Reference Nodera, Kawata, Onuma, Limpouch, Klimo and Kikuchi2008; Psikal et al. Reference Psikal, Grym, Stolcova and Proska2016), periodic nanobrush targets (Yu et al. Reference Yu, Zhou, Jin, Cao, Zhao, Hong, Li and Gu2012), targets with a symmetrical nanocone distribution (Brantov & Bychenkov Reference Brantov and Bychenkov2013) and grating surfaces (Andreev et al. Reference Andreev, Kumar, Platonov and Pukhov2011; Sgattoni et al. Reference Sgattoni, Sinigardi, Fedeli, Pegoraro and Macchi2015; Andreev et al. Reference Andreev, Platonov, Braenzel, Lübcke, Das, Messaoudi, Grunwald, Gray, McGlynn and Schnürer2016; Blanco et al. Reference Blanco, Flores-Arias, Ruiz and Vranic2017). Besides, foils with periodic surface structures have been experimentally shown to yield a twofold increase in proton energy (Margarone et al. Reference Margarone, Klimo, Kim, Prokůpek, Limpouch, Jeong, Mocek, Pšikal, Kim and Proška2012; Ceccotti et al. Reference Ceccotti, Floquet, Sgattoni, Bigongiari, Klimo, Raynaud, Riconda, Heron, Baffigi and Labate2013). However, little attention has been paid so far to the potential of non-periodic structures, even if the addition of nanospheres with irregular diameters on the front side of the target has been studied (Klimo et al. Reference Klimo, Psikal, Limpouch, Proska, Novotny, Ceccotti, Floquet and Kawata2011), as have recently been irregularities in grating targets (Blanco, Flores-Arias & Vranic Reference Blanco, Flores-Arias and Vranic2019) and surfaces coated with foams or nanowires (Zigler et al. Reference Zigler, Eisenman, Botton, Nahum, Schleifer, Baspaly, Pomerantz, Abicht, Branzel and Priebe2013; Fedeli et al. Reference Fedeli, Formenti, Cialfi, Pazzaglia and Passoni2018). Importantly, relaxing the constraint on the structure periodicity would enable simpler and more robust target fabrication methods (Langhammer, Kasemo & Zoric Reference Langhammer, Kasemo and Zoric2007; Zigler et al. Reference Zigler, Palchan, Bruner, Schleifer, Eisenmann, Botton, Henis, Pikuz, Faenov and Gordon2011), as is required to bring laser-driven proton sources closer to applications.

In this paper, we investigate by means of particle-in-cell (PIC) simulations the potential of non-periodically structured targets to enhance TNSA. Two target types are considered, consisting of a flat foil either coated on the front surface with randomly positioned nanocones (‘nanocone targets’) or perforated by nanoholes (‘nanohole targets’). In both cases, the proton cutoff energy is increased by up to a factor of two compared with flat foils. The causes of this improvement are a more complex interaction geometry, combined with locally enhanced charge-separation fields, and a larger effective laser–matter interaction area.

Nanohole targets which, in the case of relatively low areal density, yield the highest enhancement, further lead to efficient electrostatic confinement, hence lengthening the interaction duration and sustaining a high hot-electron density. In § 2, we describe the physical and numerical set-up and in § 3 we compare the results obtained with the structured targets to the flat-foil case and investigate the origin of proton energy enhancement. In § 4, a parametric study of the nanohole targets is presented. Finally, we summarize our results in § 5.

2 Physical and numerical set-ups

We will investigate how the aforementioned two types of structured foils behave with respect to proton acceleration by means of two-dimensional PIC simulations with the smilei code (Derouillat et al. Reference Derouillat, Beck, Pérez, Vinci, Chiaramello, Grassi, Flé, Bouchard, Plotnikov and Aunai2018). The considered structures are visualized in figure 1(ac). The reference target is a flat foil of thickness $d=100~\text{nm}$. It is composed of cold gold atoms, assumed to be 11-times ionized, with an ion number density of $5.85\times 10^{22}~\text{cm}^{-3}$. Simulations with different ionization states of gold atoms (up to 29-times ionized) were also performed but did not result in significant differences in proton spectra. A 10 nm thin proton–electron plasma layer with a number density of $1.74\times 10^{23}~\text{cm}^{-3}$ is added on the foil surfaces to model the hydrogen contaminants.

Figure 1. Schematic representation of the flat foil (a), nanocone (b) and nanohole (c) targets. The targets are irradiated by the laser pulse from the left under an incidence angle of 45°.

The nanocone target, sketched in figure 1(b), is composed of the above-mentioned flat foil coated with a distribution of cones, each having an opening angle of 44° and a variable base size of $b$ (specific values will be set below). The nanohole target, displayed in figure 1(c), consists of the reference foil pierced by holes of the same width and location as the above nanocones. The surfaces of both types of structured targets are coated with a proton layer.

For the initialization of the nanocones’ positions, we choose the following model. Given the position $y_{1}$ of the first nanocone, the average distance $D=2.24b$ – motivated by practical considerations for target manufacturing – and the distance spread parameter $s=D/2$, the position of the $i$th nanocone is set iteratively to $y_{i}=y_{i-1}+D+(u-1/2)s$ as long as $y_{i}<y_{\max }$, where $y_{\max }$ is the maximal position and $u$ is a uniformly distributed random number taking values from 0 to 1. This leads to an average areal cone density of $\unicode[STIX]{x1D70C}=b/D=44\,\%$. The targets are located at $x=20~\unicode[STIX]{x03BC}\text{m}$.

The $p$-polarized laser pulse has a wavelength of $\unicode[STIX]{x1D706}=0.8~\unicode[STIX]{x03BC}\text{m}$ and a maximum intensity of $I_{0}=5\times 10^{19}~\text{W}~\text{cm}^{-2}$, corresponding to a normalized vector potential $a_{0}=4.8$ and a peak electric field of $194~\text{GV}~\text{cm}^{-1}$. It has Gaussian space and time profiles with a focal spot of $5~\unicode[STIX]{x03BC}\text{m}$ FWHM (full width half maximum) and a duration of 38 fs FWHM. It is incident on the target from the lower left-hand side, at 45° from the surface normal. The peak intensity of the laser pulse reaches the foil after 175 fs in the simulation. In the following, we will set this as the time reference $t=0$.

The numerical discretization of the simulations is $\unicode[STIX]{x1D6FF}x=\unicode[STIX]{x1D6FF}y=5~\text{nm}$ and $\unicode[STIX]{x1D6FF}t=11.7~\text{as}$. We use 100 macro-particles per cell and per species in the bulk plasma, while the surface proton–electron layers are represented by 1000 macro-particles per cell and per species.

3 Enhancement of electron heating and ion acceleration

Periodic cone structures have been shown to enhance proton acceleration due to a modification of electron trajectories, hence maximizing laser absorption (Blanco et al. Reference Blanco, Flores-Arias, Ruiz and Vranic2017). It is thus interesting to investigate whether completely relaxing the restriction of periodicity would impact the acceleration process. Here, we consider a non-periodic arrangement of the structures. The cone and hole-base size in this section is fixed to be $b=300~\text{nm}$.

The strength of the rear-side sheath field that underpins TNSA, and which therefore determines the efficiency of the latter, is controlled by the energy density of the laser-generated hot electrons (Daido et al. Reference Daido, Nishiuchi and Pirozhkov2012; Macchi et al. Reference Macchi, Borghesi and Passoni2013). Figure 2(a) plots the energy spectrum of the electrons located in the vacuum region behind the rear side of the target, recorded at $t=175~\text{fs}$. Those electrons mainly account for the generation of the sheath electric field in the early stages of TNSA, when the approximation of the one-dimensional plasma expansion holds. Compared with the flat foil, both structured targets lead to a significantly increased number of relativistic electrons above 7 MeV, with the nanohole target yielding the largest enhancement – by approximately an order of magnitude.

Figure 2. Energy spectra of electrons (a) and rear-side protons (b) from a flat foil (solid black), a nanocone (dashed-dotted yellow) and a nanohole (dashed red) target of $d=100~\text{nm}$ thickness. In (a), only electrons located in the vacuum behind the target backside are considered. Electron spectra are recorded at $t=175~\text{fs}$ and proton spectra at $t=455~\text{fs}$.

Figure 3. Energy spectra of the electrons (a) and rear-protons (b) from a flat foil (solid blue), a $d=100~\text{nm}$ periodic nanohole target (dashed-dotted black) and a $d=100~\text{nm}$ non-periodic nanohole target (dashed red). In (a), only electrons located in the vacuum behind the target backside are considered. Electron spectra are recorded at $t=175~\text{fs}$ and proton spectra at $t=455~\text{fs}$.

Figure 4. (a)–(c) Maximum $\unicode[STIX]{x1D6FE}$ factor (locally averaged over the initial particle distribution) reached by electrons as a function of their initial position (yellow–red) and longitudinal electric field $E_{x}$ (units $m_{e}c\unicode[STIX]{x1D714}_{0}/e$, with $c$ the speed of light, $\unicode[STIX]{x1D714}_{0}$ the laser frequency, $m_{e}$ the electron mass and $e$ the elementary charge) at $t=0$, averaged over a laser period $T_{0}$ (blue–red). (d)–(f) Magnetic field $B_{z}$ (units $m_{e}\unicode[STIX]{x1D714}_{0}/e$) at $t=0$. The initial target densities are indicated in light grey, and correspond to $100~\text{nm}$-thick flat (a,d) nanohole (b,e) and nanocone targets (c,f) and the center of the laser spot on target is localized at $y=33.5~\unicode[STIX]{x03BC}\text{m}$.

A similar behaviour is found for the rear-side proton energy spectra, as plotted in figure 2(b) at $t=455~\text{fs}$. Both nanocone and nanohole targets give rise to much enhanced TNSA: the best performance is observed using nanoholes, with an almost doubled proton cutoff energy compared to that from the flat foil (${\sim}30~\text{MeV}$ protons versus ${\sim}17~\text{MeV}$).

Figure 3 demonstrates that the electron and proton spectra are almost the same for periodic and non-periodic targets, consistent with the conclusion reached by Klimo et al. (Reference Klimo, Psikal, Limpouch, Proska, Novotny, Ceccotti, Floquet and Kawata2011) when investigating irregular spherical structures. It is thus to be expected that the individual structuring units, rather than their periodic arrangement (possibly leading to the excitation of surface plasma waves), are responsible for the enhancement. This is a favourable result from an experimental perspective since it reduces the target manufacturing constraints.

Moreover, figure 3(a,b) presents the electron and rear-proton energy spectra obtained for a flat foil (blue solid line) of reduced thickness ($d=55~\text{nm}$) such that it contains the same total amount of matter as the $d=100~\text{nm}$ nanohole targets, whether periodic or not. This thinner flat foil produces particle energy spectra very similar to the $d=100~\text{nm}$ foil; hence, the enhanced performance of the nanohole targets cannot be ascribed to the direct effect of their reduced volume (or area in two dimensions) on the electron kinetic energy density (which would naturally increase assuming the same amount of laser energy is converted into hot electrons), but rather results from a strongly modified hot-electron dynamics.

Figure 5. Trajectories of a few energetic electrons in the nanohole target (a) and the flat target (b). These electrons are selected randomly amongst those verifying $x>48~\unicode[STIX]{x03BC}\text{m}$ at $t=175~\text{fs}$, and as long as they have reached a threshold energy $E_{\text{th}}$ during the simulation, with $E_{\text{th}}=12~\text{MeV}$ in the nanohole case and $E_{\text{th}}=5~\text{MeV}$ in the flat target case. The initial target density is indicated in light grey. The colour of the trajectories (blue–red) represents the rate of change ($(p_{x}E_{x}+p_{y}E_{y})/\unicode[STIX]{x1D6FE}$) of electron energy (units $m_{e}^{2}c^{2}\unicode[STIX]{x1D714}_{0}/e$). Two trajectories of single electrons (black) are represented for each case in the different insets. The red dots indicate the initial positions of these particles.

In order to gain insight into the electron energization process in the nanostructured target, we record the maximum Lorentz factor $\unicode[STIX]{x1D6FE}$ reached by each electron during the simulation, and plot its locally averaged value as a function of the initial electron position (figure 4ac). In the case of the flat foil (figure 4a), the resulting map shows that, as expected, most of the accelerated electrons originate from a ${\sim}25~\text{nm}$-thick surface layer at the directly irradiated front side of the target (with mean energies ${\sim}6~\text{MeV}$ being reached), and that rear-side electrons undergo negligible acceleration. In the case of the nanohole (figure 4b) and nanocone (figure 4c) targets, some of the highest energy electrons stem from additional regions, namely the nanohole walls and the nanocone sides. The nanostructuring of the target surface therefore increases the effective interaction area, leading to a larger number of hot electrons. The mean energy (${\sim}10~\text{MeV}$) reached by these electrons is also significantly larger than in the flat foil. These effects are supported by the local enhancement of the electrostatic field which appears when using nanostructured targets. The electrostatic field is indeed strongly enhanced at the corners of the nanoholes and at the tips of the nanocones, which correspond to the surfaces from which the most energetic electrons arise. Note that the value of the electrostatic field in these regions becomes of the same order as the laser field (figure 4df), which is also favourable for electron acceleration (Paradkar, Krasheninnikov & Beg Reference Paradkar, Krasheninnikov and Beg2012). Interestingly, when the target is nanostructured, the change in the laser field pattern on the front side of the target indicates a decrease in the laser reflection compared with the flat foil. This is confirmed by the reflection coefficient, which drops from $C_{R_{f}}=0.67$ to $C_{R_{c}}=0.15$ when adding the nanocones. While the flat and nanocone targets remain essentially opaque (with a transmission coefficient $C_{T_{f,c}}<0.02$), the laser light is partially transmitted through the nanoholes ($C_{T_{h}}=0.13$), as is evident from figure 4. Consequently, the absorption coefficient is lower in nanohole ($C_{A_{h}}=0.66$) than in nanocone targets ($C_{A_{c}}=0.83$). As further discussed below, this finding somewhat contradicts the widely shared notion that the absorbed laser energy fraction is the single figure of merit for TNSA.

This general behaviour can be complemented by examining individual electron trajectories. We focus on those high-energy electrons breaking through the target rear side, and therefore contributing to the accelerating sheath field. For this reason, figure 5 plots the trajectories of a sample of the most energetic electrons in both the flat and nanohole targets – with a lower energy cutoff for the selection of 12 MeV in the nanohole target (figure 5a) and of 5 MeV in the flat target case (figure 5b). Only those electrons that are located at a longitudinal position $x\geqslant 48~\unicode[STIX]{x03BC}\text{m}$ at $t=175~\text{fs}$ are selected. The colour of each trajectory is indexed on the rate of change of electron energy in the local electromagnetic fields. It can be seen that, on average, higher values are attained in the nanohole target. While the selected electrons mainly exhibit acceleration in the front- and rear-side vacuum regions, sizable energy transfer is also seen within the cavities, which is consistent with the laser being partially transmitted through the nanoholes (see figure 4e).

Figure 6. Position of the front of the accelerating rear-proton layer at $t=175~\text{fs}$ for a flat foil (solid black), a nanocone (dashed-dotted yellow) and a nanohole (dashed red) targets of $d=100~\text{nm}$ thickness.

The sub-micron dense regions that make up the nanohole targets effectively behave as mass-limited targets (Psikal et al. Reference Psikal, Tikhonchuk, Limpouch, Andreev and Brantov2008; Buffechoux et al. Reference Buffechoux, Psikal, Nakatsutsumi, Romagnani, Andreev, Zeil, Amin, Antici, Burris-Mog and Compant-La-Fontaine2010), leading to efficient electrostatic confinement of the hot electrons. This is evidenced by the single particle trajectories displayed in the insets of figure 5: the nanohole target causes the electron to recirculate in both longitudinal (across the front and rear sides of the target) and transverse (across the nanohole walls) directions. A favourable consequence of this is a longer effective laser–electron interaction time. Also, the laser–electron interaction occurs under various geometrical conditions, and so with increased degrees of freedom. This results in a more complete exploration of the phase space, which ultimately allows the electrons to be accelerated to higher energies. Finally, being prevented from leaving the laser-irradiated region, the hot electrons are able to sustain a strong sheath field over longer times. The confinement of the hot electrons entails a reduced transverse extent of the sheath field, and therefore of the expanding proton cloud. This can be seen in figure 6, which shows the boundary of the proton cloud as resulting from the three target types, the nanohole target giving rise to a more narrow proton distribution. Compared to nanocone targets, the more localized and intense sheath fields induced in nanohole targets translate into lower numbers of accelerated protons but of higher energies, despite lower laser absorption. A drawback of nanohole targets, however, is the significant transverse sheath fields which build up at the sharp transverse vacuum–plasma interfaces, and tend to increase the divergence of the proton beam (measured to be approximately twice higher in nanohole targets than in flat foils).

4 Parametric scan for nanohole targets

Figure 7. Particle density spectra of the rear-side electrons (a,d,g), front-side protons (b,e,h) and rear-side protons (c,f,i) from nanohole targets. The electron spectra are presented at $t=175~\text{fs}$, the front proton spectra at $t=315~\text{fs}$ and the rear proton spectra at $t=455~\text{fs}$. The foil thickness is $d=100~\text{nm}$ in (ac), $d=300~\text{nm}$ in (df), $d=600~\text{nm}$ in (gi). The nanohole diameter is varied in the range $0\leqslant b\leqslant 600~\text{nm}$, as indicated in the legend of each row.

We now perform a parametric scan where we vary the foil thickness $d$ and hole diameter $b$. Henceforth, we consider only the cases where the hole diameter is at least as large as the foil thickness, i.e. $b\geqslant d$.

Figure 7(a) presents the spectra of the rear-side electrons for foil thicknesses from 100 to 600 nm and various hole diameters. One can see that the electrons from the flat foil (grey line) reach lower energies than electrons from perforated foils, whatever the hole diameter. Especially above 5 MeV, perforated targets produce more high-energy electrons. The enhancement is very similar for hole sizes from 100 to 600 nm.

While we are mostly interested in the protons accelerated from the target back side, we also plot for completeness the energy spectra of the protons originating from the front side. Figure 7(b,c) presents the corresponding front (b) and rear (c) proton spectra. In both cases, the structuring significantly enhances the proton cutoff energies. While a 100 nm hole size already enhances the proton cutoff energy by approximately 40 %–50 %, the largest enhancements are reached for 300 to 600 nm diameter holes.

When increasing the gold foil thickness, the trends remain the same (see figure 7di). The presence of the holes increases the number of high-energy electrons and boosts the proton energy by approximately a factor of 2. The thinner the gold foil, the higher the proton energies. However, while this nanohole structuring can be used for all target thicknesses, the relative improvement is slightly more pronounced for thicker foils (by 95 % at $d=300~\text{nm}$ versus 70 % at $d=100~\text{nm}$).

The PIC simulation results suggest that as long as the parameters are in the above-described range leading to enhancement, they can be freely chosen depending on other experimental constraints, e.g. ease of handling and fabrication.

Note that our simulations consider sharp-gradient targets, meaning that they do not address the possible influence of finite preplasmas caused by laser prepulses. The generation of a significant preplasma would modify the picture: the holes might be filled with electrons before the arrival of the main pulse, thus suppressing the benefit of the target structuring. Larger hole diameters may then be preferable under actual experimental conditions.

Figure 8. Rear-side electron and proton energy spectra under conditions similar to those of figure 2, but with different hole densities $\unicode[STIX]{x1D70C}$ according to the legends in (a) and (b).

The areal density of the holes can be expected to play a role in enhancing TNSA. Figure 8 displays the energy spectra of the rear-side electrons and protons obtained from 100 nm thick nanohole targets of nanohole density varying in the range $0\leqslant \unicode[STIX]{x1D70C}\leqslant 67\,\%$. An optimum electron heating and proton acceleration is reached when almost half of the surface area is covered by the holes. However, already for a density of approximately 20 %, a notable increase in the cutoff energy is reached.

The optimum is likely formed by two counteracting effects. On the one hand, when increasing the hole density, the mechanisms described in § 3 further develop: more electrons can get an energy boost due to the increase in the interaction surface and recirculation. On the other hand, raising the hole density reduces the effective target volume, eventually resulting in higher laser transmission, and therefore weaker laser–target coupling. To mitigate the latter effect, one could use thicker foils instead of the considered 100 nm thin foils. However, this solution is limited by the fact that the proton cutoff energy tends to be reduced when the nanohole target gets thicker – in the same manner as for the standard flat-foil target.

5 Conclusion

In summary, by means of two-dimensional PIC simulations, we investigate laser-driven proton acceleration from non-periodic nanohole and nanocone targets. We demonstrate a significant increase in the proton cutoff energy in both types of structured targets compared to flat foils. This better performance is found to originate from several factors: an enlarged effective interaction surface between the laser and the target, a larger number of degrees of freedom in the electron dynamics and a local intensification of the charge-separation fields. In the case of nanohole targets, we identify a large parameter space in terms of hole diameter, foil thickness and hole areal density yielding significant enhancement of the accelerated proton spectra. The specific nanohole geometry yields additional control through the limitation of the hot electrons’ transverse expansion. Interestingly, the enhancement of the ion acceleration does not stem only from the increased laser-to-electron conversion efficiency, since nanohole targets yield higher proton energies than nanocone targets despite a lower laser absorption.

Finally, our results show that the production of structured targets for improved ion acceleration can be relaxed to non-periodic structures with a relatively low areal density of structuring units.

Acknowledgements

The authors acknowledge fruitful discussions with L. Yi and the rest of the PLIONA team. This work was supported by the Knut and Alice Wallenberg Foundation, the European Research Council (ERC-2014-CoG grant 647121) and by the Swedish Research Council, grants no. 2016-05012 and no. 2017-04684. The simulations were performed on resources provided by GENCI (project A0060507594) and at Chalmers Centre for Computational Science and Engineering (C3SE) provided by the Swedish National Infrastructure for Computing (grants SNIC2018-2-13 and SNIC2018-3-297).

References

Andreev, A., Kumar, N., Platonov, K. & Pukhov, A. 2011 Efficient generation of fast ions from surface modulated nanostructure targets irradiated by high intensity short-pulse lasers. Phys. Plasmas 18 (10), 103103.CrossRefGoogle Scholar
Andreev, A., Platonov, K., Braenzel, J., Lübcke, A., Das, S., Messaoudi, H., Grunwald, R., Gray, C., McGlynn, E. & Schnürer, M. 2016 Relativistic laser nano-plasmonics for effective fast particle production. Plasma Phys. Control. Fusion 58 (1), 014038.CrossRefGoogle Scholar
Blanco, M., Flores-Arias, M. T., Ruiz, C. & Vranic, M. 2017 Table-top laser-based proton acceleration in nanostructured targets. New J. Phys. 19, 033004.CrossRefGoogle Scholar
Blanco, M., Flores-Arias, M. & Vranic, M.2019 Suitability and robustness of triangular nanostructured targets for proton acceleration. arXiv:1902.05641.Google Scholar
Brantov, A. & Bychenkov, V. Y. 2013 Laser-triggered proton acceleration from micro-structured thin targets. Contributions Plasma Phys. 53 (10), 731735; arXiv:https://onlinelibrary.wiley.com/doi/pdf/10.1002/ctpp.201310046.CrossRefGoogle Scholar
Brenner, C. M., Robinson, A. P. L., Markey, K., Scott, R. H. H., Gray, R. J., Rosinski, M., Deppert, O., Badziak, J., Batani, D., Davies, J.R. et al. 2014 High energy conversion efficiency in laser-proton acceleration by controlling laser-energy deposition onto thin foil targets. Appl. Phys. Lett. 104 (8), 081123.CrossRefGoogle Scholar
Buffechoux, S., Psikal, J., Nakatsutsumi, M., Romagnani, L., Andreev, A., Zeil, K., Amin, M., Antici, P., Burris-Mog, T., Compant-La-Fontaine, A. et al. 2010 Hot electrons transverse refluxing in ultraintense laser–solid interactions. Phys. Rev. Lett. 105, 015005.CrossRefGoogle ScholarPubMed
Bulanov, S. V. & Khoroshkov, V. S. 2002 Feasibility of using laser ion accelerators in proton therapy. Plasma Phys. Rep. 28, 453456.CrossRefGoogle Scholar
Ceccotti, T., Floquet, V., Sgattoni, A., Bigongiari, A., Klimo, O., Raynaud, M., Riconda, C., Heron, A., Baffigi, F., Labate, L. et al. 2013 Evidence of resonant surface-wave excitation in the relativistic regime through measurements of proton acceleration from grating targets. Phys. Rev. Lett. 111, 185001.CrossRefGoogle ScholarPubMed
Daido, H., Nishiuchi, M. & Pirozhkov, A. S. 2012 Review of laser-driven ion sources and their applications. Rep. Prog. Phys. 75 (5), 056401.CrossRefGoogle ScholarPubMed
Derouillat, J., Beck, A., Pérez, F., Vinci, T., Chiaramello, M., Grassi, A., Flé, M., Bouchard, G., Plotnikov, I., Aunai, N. et al. 2018 Smilei: a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation. Comput. Phys. Commun. 222, 351373.CrossRefGoogle Scholar
Fedeli, L., Formenti, A., Cialfi, L., Pazzaglia, A. & Passoni, M. 2018 Ultra-intense laser interaction with nanostructured near-critical plasmas. Sci. Rep. 8 (1), 3834.CrossRefGoogle ScholarPubMed
Ferri, J., Siminos, E. & Fülöp, T. 2019 Enhanced target normal sheath acceleration using colliding laser pulses. Commun. Phys. 2, 40.CrossRefGoogle Scholar
Habs, D., Thirolf, P. G., Gross, M., Allinger, K., Bin, J., Henig, A., Kiefer, D., Ma, W. & Schreiber, J. 2011 Introducing the fission–fusion reaction process: using a laser-accelerated Th beam to produce neutron-rich nuclei towards the $n=126$ waiting point of the $r$-process. Appl. Phys. B 103 (2), 471484.CrossRefGoogle Scholar
Higginson, A., Gray, R. J., King, M., Dance, R. J., Williamson, S. D. R., Butler, N. M. H., Wilson, R., Capdessus, R., Armstrong, C., Green, J. S. et al. 2018 Near-100 MeV protons via a laser-driven transparency-enhanced hybrid acceleration scheme. Nature Commun. 9 (1), 724.CrossRefGoogle Scholar
Kaluza, M., Schreiber, J., Santala, M. I. K., Tsakiris, G. D., Eidmann, K., Meyer-ter Vehn, J. & Witte, K. J. 2004 Influence of the laser prepulse on proton acceleration in thin-foil experiments. Phys. Rev. Lett. 93, 045003.CrossRefGoogle ScholarPubMed
Klimo, O., Psikal, J., Limpouch, J., Proska, J., Novotny, F., Ceccotti, T., Floquet, V. & Kawata, S. 2011 Short pulse laser interaction with micro-structured targets: simulations of laser absorption and ion acceleration. New J. Phys. 13 (5), 053028.CrossRefGoogle Scholar
Langhammer, C., Kasemo, B. & Zoric, I. 2007 Absorption and scattering of light by Pt, Pd, Ag, and Au nanodisks: absolute cross sections and branching ratios. J. Chem. Phys. 126 (19), 194702.CrossRefGoogle ScholarPubMed
Macchi, A., Borghesi, M. & Passoni, M. 2013 Ion acceleration by superintense laser–plasma interaction. Rev. Mod. Phys. 85, 751.CrossRefGoogle Scholar
Margarone, D., Klimo, O., Kim, I. J., Prokůpek, J., Limpouch, J., Jeong, T. M., Mocek, T., Pšikal, J., Kim, H. T., Proška, J. et al. 2012 Laser-driven proton acceleration enhancement by nanostructured foils. Phys. Rev. Lett. 109, 234801.CrossRefGoogle ScholarPubMed
Neely, D., Foster, P., Robinson, A., Lindau, F., Lundh, O., Persson, A., Wahlström, C.-G. & McKenna, P. 2006 Enhanced proton beams from ultrathin targets driven by high contrast laser pulses. Appl. Phys. Lett. 89 (2), 021502.CrossRefGoogle Scholar
Nodera, Y., Kawata, S., Onuma, N., Limpouch, J., Klimo, O. & Kikuchi, T. 2008 Improvement of energy-conversion efficiency from laser to proton beam in a laser–foil interaction. Phys. Rev. E 78, 046401.Google Scholar
Nuter, R., Gremillet, L., Combis, P., Drouin, M., Lefebvre, E., Flacco, A. & Malka, V. 2008 Influence of a preplasma on electron heating and proton acceleration in ultraintense laser–foil interaction. J. Appl. Phys. 104, 103307.CrossRefGoogle Scholar
Ogura, K., Nishiuchi, M., Pirozhkov, A. S., Tanimoto, T., Sagisaka, A., Esirkepov, T. Z., Kando, M., Shizuma, T., Hayakawa, T., Kiriyama, H. et al. 2012 Proton acceleration to 40 MeV using a high intensity, high contrast optical parametric chirped-pulse amplification/ti: sapphire hybrid laser system. Opt. Lett. 37 (14), 28682870.CrossRefGoogle ScholarPubMed
Paradkar, B. S., Krasheninnikov, S. I. & Beg, F. N. 2012 Mechanism of heating of pre-formed plasma electrons in relativistic laser–matter interaction. Phys. Plasmas 19 (6), 060703.CrossRefGoogle Scholar
Psikal, J., Grym, J., Stolcova, L. & Proska, J. 2016 Hollow target for efficient generation of fast ions by ultrashort laser pulses. Phys. Plasmas 23 (12), 123121.CrossRefGoogle Scholar
Psikal, J., Tikhonchuk, V. T., Limpouch, J., Andreev, A. A. & Brantov, A. V. 2008 Ion acceleration by femtosecond laser pulses in small multispecies targets. Phys. Plasmas 15 (5), 053102.CrossRefGoogle Scholar
Romagnani, L., Fuchs, J., Borghesi, M., Antici, P., Audebert, P., Ceccherini, F., Cowan, T., Grismayer, T., Kar, S., Macchi, A. et al. 2005 Dynamics of electric fields driving the laser acceleration of multi-MeV protons. Phys. Rev. Lett. 95, 195001.CrossRefGoogle ScholarPubMed
Sgattoni, A., Sinigardi, S., Fedeli, L., Pegoraro, F. & Macchi, A. 2015 Laser-driven Rayleigh–Taylor instability: plasmonic effects and three-dimensional structures. Phys. Rev. E 91, 013106.Google ScholarPubMed
Snavely, R. A., Key, M. H., Hatchett, S. P., Cowan, T. E., Roth, M., Phillips, T. W., Stoyer, M. A., Henry, E. A., Sangster, T. C., Singh, M. S. et al. 2000 Intense high-energy proton beams from petawatt-laser irradiation of solids. Phys. Rev. Lett. 85, 29452948.CrossRefGoogle ScholarPubMed
Wagner, F., Deppert, O., Brabetz, C., Fiala, P., Kleinschmidt, A., Poth, P., Schanz, V. A., Tebartz, A., Zielbauer, B., Roth, M. et al. 2016 Maximum proton energy above 85 MeV from the relativistic interaction of laser pulses with micrometer thick $\text{CH}_{2}$ targets. Phys. Rev. Lett. 116, 205002.CrossRefGoogle ScholarPubMed
Wilks, S. C., Langdon, A. B., Cowan, T. E., Roth, M., Singh, M., Hatchett, S., Key, M. H., Pennington, D., MacKinnon, A. & Snavely, R. A. 2001 Energetic proton generation in ultra-intense laser–solid interactions. Phys. Plasmas 8 (2), 542549.CrossRefGoogle Scholar
Yu, J., Zhou, W., Jin, X., Cao, L., Zhao, Z., Hong, W., Li, B. & Gu, Y. 2012 Improvement of proton energy in high-intensity laser-nanobrush target interactions. Laser Part. Beams 30 (2), 307311.CrossRefGoogle Scholar
Zigler, A., Eisenman, S., Botton, M., Nahum, E., Schleifer, E., Baspaly, A., Pomerantz, I., Abicht, F., Branzel, J., Priebe, G. et al. 2013 Enhanced proton acceleration by an ultrashort laser interaction with structured dynamic plasma targets. Phys. Rev. Lett. 110, 215004.CrossRefGoogle ScholarPubMed
Zigler, A., Palchan, T., Bruner, N., Schleifer, E., Eisenmann, S., Botton, M., Henis, Z., Pikuz, S. A., Faenov, A. Y., Gordon, D. et al. 2011 5.5–7.5 MeV proton generation by a moderate-intensity ultrashort-pulse laser interaction with $\text{H}_{2}\text{O}$ nanowire targets. Phys. Rev. Lett. 106, 134801.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Schematic representation of the flat foil (a), nanocone (b) and nanohole (c) targets. The targets are irradiated by the laser pulse from the left under an incidence angle of 45°.

Figure 1

Figure 2. Energy spectra of electrons (a) and rear-side protons (b) from a flat foil (solid black), a nanocone (dashed-dotted yellow) and a nanohole (dashed red) target of $d=100~\text{nm}$ thickness. In (a), only electrons located in the vacuum behind the target backside are considered. Electron spectra are recorded at $t=175~\text{fs}$ and proton spectra at $t=455~\text{fs}$.

Figure 2

Figure 3. Energy spectra of the electrons (a) and rear-protons (b) from a flat foil (solid blue), a $d=100~\text{nm}$ periodic nanohole target (dashed-dotted black) and a $d=100~\text{nm}$ non-periodic nanohole target (dashed red). In (a), only electrons located in the vacuum behind the target backside are considered. Electron spectra are recorded at $t=175~\text{fs}$ and proton spectra at $t=455~\text{fs}$.

Figure 3

Figure 4. (a)–(c) Maximum $\unicode[STIX]{x1D6FE}$ factor (locally averaged over the initial particle distribution) reached by electrons as a function of their initial position (yellow–red) and longitudinal electric field $E_{x}$ (units $m_{e}c\unicode[STIX]{x1D714}_{0}/e$, with $c$ the speed of light, $\unicode[STIX]{x1D714}_{0}$ the laser frequency, $m_{e}$ the electron mass and $e$ the elementary charge) at $t=0$, averaged over a laser period $T_{0}$ (blue–red). (d)–(f) Magnetic field $B_{z}$ (units $m_{e}\unicode[STIX]{x1D714}_{0}/e$) at $t=0$. The initial target densities are indicated in light grey, and correspond to $100~\text{nm}$-thick flat (a,d) nanohole (b,e) and nanocone targets (c,f) and the center of the laser spot on target is localized at $y=33.5~\unicode[STIX]{x03BC}\text{m}$.

Figure 4

Figure 5. Trajectories of a few energetic electrons in the nanohole target (a) and the flat target (b). These electrons are selected randomly amongst those verifying $x>48~\unicode[STIX]{x03BC}\text{m}$ at $t=175~\text{fs}$, and as long as they have reached a threshold energy $E_{\text{th}}$ during the simulation, with $E_{\text{th}}=12~\text{MeV}$ in the nanohole case and $E_{\text{th}}=5~\text{MeV}$ in the flat target case. The initial target density is indicated in light grey. The colour of the trajectories (blue–red) represents the rate of change ($(p_{x}E_{x}+p_{y}E_{y})/\unicode[STIX]{x1D6FE}$) of electron energy (units $m_{e}^{2}c^{2}\unicode[STIX]{x1D714}_{0}/e$). Two trajectories of single electrons (black) are represented for each case in the different insets. The red dots indicate the initial positions of these particles.

Figure 5

Figure 6. Position of the front of the accelerating rear-proton layer at $t=175~\text{fs}$ for a flat foil (solid black), a nanocone (dashed-dotted yellow) and a nanohole (dashed red) targets of $d=100~\text{nm}$ thickness.

Figure 6

Figure 7. Particle density spectra of the rear-side electrons (a,d,g), front-side protons (b,e,h) and rear-side protons (c,f,i) from nanohole targets. The electron spectra are presented at $t=175~\text{fs}$, the front proton spectra at $t=315~\text{fs}$ and the rear proton spectra at $t=455~\text{fs}$. The foil thickness is $d=100~\text{nm}$ in (ac), $d=300~\text{nm}$ in (df), $d=600~\text{nm}$ in (gi). The nanohole diameter is varied in the range $0\leqslant b\leqslant 600~\text{nm}$, as indicated in the legend of each row.

Figure 7

Figure 8. Rear-side electron and proton energy spectra under conditions similar to those of figure 2, but with different hole densities $\unicode[STIX]{x1D70C}$ according to the legends in (a) and (b).