1. Introduction
Since their inception in the early 2000s, laser-driven ion beams have opened novel pathways compared to those powered by traditional accelerators [
1,
2,
3]. Owing to their very compact size (
), ultrafast timescales (
) and very high currents (
) at the source, these beams enable a wide range of applications, spanning flash radiography [
4], radiobiology [
5], nuclear physics [
6], and even ion therapy [
7]. They can also be leveraged to generate bright neutron sources [
8,
9,
10,
11] which already outperform traditional neutron sources in terms of their intensity [
12] and hold promise for material probing [
13] and astrophysics [
14].
Over the past 20 years, research has shown that laser pulses focused at relativistic intensities (
) onto thin solid foils (e.g. a few
m thick) can generate fast (i.e., up to
nowadays [
15]) and highly directional ion beams. This process, named target normal sheath acceleration (TNSA) [
16], relies on the strong (
) electrostatic sheath fields induced by the laser-generated hot electrons [
17] at the target surfaces. Typically, TNSA is most efficient at the nonirradiated target backside [
18,
19,
20], where it produces ion beams with a high particle flux (e.g.,
–
particles exceeding a few MeV) and a relatively small angular divergence (
–30 depending on the energy [
21]) [
22,
23]. Moreover, different kinds of ions can be accelerated simply by changing the target material [
24].
Now that we are at the dawn of the multi-petawatt (PW) laser era [
25], the vast body of theoretical studies [
1,
2,
26] tells us that further progress should be expected via little-explored, yet likely more efficient ion acceleration processes, i.e., radiation pressure acceleration (RPA) [
27,
28,
29], breakout-afterburner (BOA) [
30], and relativistically induced transparency (RIT) [
31,
32,
33]. Numerical simulations predict better scaling of the maximum proton energy with laser intensity when utilizing ultrathin (nanometer-range) targets. In experiments, however, triggering the above mechanisms in a controlled way remains very challenging because the amplified spontaneous emission (ASE) pedestal, extending over picosecond or sometimes nanosecond scales ahead of the peak of the laser pulse, can easily preheat or even fully expand the thin target (less than few micron) [
22], resulting in weak or no acceleration at all.
Plasma mirrors (PMs) are a well-known solution to eliminate the ASE, and hence to improve the temporal laser contrast, defined as the ratio of the peak intensity to the ASE intensity [
34,
35]. PMs are commonly made of polished glass slabs or liquid crystals [
36] through which the low-intensity ASE light is transmitted. The subsequent rise in laser intensity eventually causes the irradiated PM surface to be ionized and, once the local electron density becomes overcritical, to reflect the main pulse. The resulting higher laser contrast can then allow ultrathin foils to survive until the intensity peak, leading to interaction conditions suitable for efficient proton acceleration.
For example, using a PM to achieve a high laser contrast of
, laser pulses of 40 fs duration and 1 J energy have been reported to drive protons up to 10 MeV energies from 20–30 nm thick foils [
37,
38]. Yet such ultrathin targets tend to result in filamentary proton beam profiles, as compared to the smooth profiles obtained with micrometer-range targets [
34,
39,
40]. Concerning energy conversion efficiency, standard TNSA with
thick foils is typically limited to
efficiency for joule-class lasers [
3,
39], which can be boosted up to
if a hybrid TNSA/RPA regime is realized using thinner foils [
41]. In the case of nanometer-range liquid crystal targets driven in the RIT regime by 3 J energy, 30 fs laser pulses, 20 MeV proton cutoff energies have been measured, but with degraded beam profile at low proton energy [
42].
A general observation is that relatively long (a few 100 fs) pulse durations often lead to poor proton beam quality but higher conversion efficiency [
34,
40,
41]. This is explained by the fact that the ASE-driven expansion of the target increases the plasma density scale-length seen by the main laser pulse. This tends to favor the laser absorption into hot electrons but also to disrupt the laser beam profile, entailing, in turn, a degraded imprint in the proton profile. When using a PM, the reduced preplasma due to the higher laser contrast is expected to improve the compactness of the hot-electron distribution and the proton beam quality, yet the associated laser energy loss, typically around
, may detrimentally affect the achievable proton energies. Therefore, it is not clear whether the use of PMs systematically benefits or not proton acceleration.
In this study, we experimentally demonstrate that the combination of an ultrashort (
) and ultraintense (
) laser pulse with a double plasma mirror (DPM) [
34] can significantly improve the maximum energy, conversion efficiency, beam profile, and divergence of the accelerated protons. Numerical particle-in-cell (PIC) simulations support these results, ascribing them to the more compact spatial distribution of hot electrons generated by the ASE-free laser pulse. Our findings have important implications for the development of high-quality proton beams for various applications.
2. Experimental Setup
The experiment was performed at the short-focal area of the Apollon laser facility (Saclay, France), using its F2 secondary beam [
43]. This beam utilizes a Ti:Sapphire laser, delivering pulses with a full-width-at-half-maximum (FWHM) duration of 24 fs, a mean on-target energy of 11 J and a central laser wavelength of 815 nm. Two configurations were tested: direct shots (without a DPM) and shots with a DPM [
34], in which case the on-target laser energy was reduced to 5.7 J, due a measured DPM reflectivity of
reflectivity.
The experimental setup is sketched in
Figure 1. The DPM was placed inside the interaction chamber, between the
off-axis parabola (OAP) and the target. The beam was focused onto the target with a 45incidence angle. A camera recorded the image of the reflected laser beam on a Spectralon scatter plate positioned along the specular reflection axis. The transmitted light was monitored by another Spectralon plate placed along the laser transmission axis. Aluminum foils of
–
thicknesses were used on direct shots, while thinner silicon foils, ranging from 10 nm to 300 nm (with, facing the incident laser, an additional coating of 50 nm of Al for some of them, in order to improve the laser absorption and the proton generation [
44]), were used with the DPM. The proton spectra were measured using stacks of EBT3 Gafchromic radiochromic films (RCFs) interlaced with aluminum filters, placed at 25 mm from the target.
Figure 2 provides information on the focusing of the 140-mm-diameter, near-field laser beam. The focal spot was measured to be elliptical with
FWHM dimensions (see
Figure 2a). It contained
of the on-target laser energy, resulting in a peak intensity of
. This peak intensity is derived directly from the focal spot image shown in
Figure 2a (see caption). These values are very close to those obtained on direct shots during the commissioning of the F2 beam [
43], demonstrating that the DPM did not deteriorate the laser focal spot.
Figure 2b details the laser energy distribution within the laser spot. The graph plots the fraction of the laser energy within an encircled region as a function of the circle radius from the focal spot center.
The temporal intensity profile of the laser beam (measured without a DPM) is shown in
Figure 3a. One can see that before the intensity peak (at time zero), there exists a 750 fs long plateau with a mean intensity of
. This mildly relativistic intensity level is sufficient to destroy Al targets with thickness below
[
43].
The laser irradiated the PMs at 45 incidence angle, with elliptical spots of dimensions
on the PM closest to the focusing OAP and
on the second PM. The peak intensity on the PMs reached
while the prepulse plateau intensity was of
. These values are, respectively, well above and below the ionization threshold (
, considering the
prepulse duration) of the fused silica material of the PM [
45]. The on-target intensity contrast between the prepulse plateau and pulse maximum was then increased by the square of the ratio between the reflection factor of the ionized PM (
) and the antireflection coating factor of each PM (
), i.e., by a factor of
. This would bring the prepulse plateau intensity down to
, causing minimal, if any, disturbance to the foil targets. This prediction is supported by the observed smooth spatial profiles of the protons accelerated from targets as thin as 20 nm (see below). Note that the smooth spatial profiles we observe differ from previously reported results [
34,
39,
40], which could be linked with the fact that the latter used longer laser pulses than in our case. Longer laser pulse leave indeed more time to induce deformation of the target surface during acceleration.
Figure 3b1 shows a typical image of the specularly reflected beam as obtained on the Spectralon plate without a DPM. The beam profile appears to be diffuse and of low intensity, meaning that the laser beam was reflected off an extended preplasma, as a result of the prepulse-induced target expansion. By contrast, the reflected beam profiles recorded on DPM shots (see
Figure 3b2) were much more intense and collimated, suggesting that the laser interacted with a sharp target surface. This behavior is consistent with previous experiments using DPMs [
38,
42].
Figure 3c1 displays the transmitted laser beam profile measured with a DPM in the absence of a target. A very similar image (
Figure 3c2), associated with a laser transmission rate of
, was recorded when shooting at the thinnest (10 nm) Si foils, in which case no proton acceleration was detected. This indicates that, despite the enhanced contrast achieved by the DPM, such ultrathin targets remain vulnerable to the ASE. However, when the target thickness was increased to 50 nm, fast protons were measured up to
while the laser transmission rate plummeted to
(see Spectralon image in
Figure 3c3).
4. Numerical PIC Modeling
To elucidate the mechanism behind the enhanced proton acceleration with a DPM, we have carried 2D PIC simulations with the fully relativistic and electromagnetic
smilei code [
57]. The numerical domain has dimensions of
with a mesh size
, corresponding to 128 cells per micron. The time step is
, where
c is the speed of light. The laser beam, of wavelength
, is injected from the lower left side (
) of the box at 45 incidence angle. It is focused with
p-polarization onto the front side of the target (at
and
). It has a Gaussian envelope with
FWHM spot size and 24 fs FWHM duration. In the no-DPM case, the laser peak intensity is
, corresponding to
. The laser beam also has a 1 ps long prepulse with the same spatial profile and a constant intensity.The intensity of the pre-pulse was adjusted so that the simulated proton energy spectrum matches the experimental one, see
Figure 4a. In the DPM case, we consider only the main 24 fs pulse, with no pre-pulse, but with a reduced peak intensity of
(corresponding to
) to take into account the energy loss associated with the DPM.
In the no-DPM case, the plasma targets consist of neutral Al foils with thickness –. In the DPM case, we use neutral Si foils with –300 nm, followed by a 50 nm thick neutral Al layer. All targets are coated with a 10 nm thick hydrogen contaminant layer on their backside. Due to computational limitations, their initial ion number density is capped at a sub-solid density value, i.e., , where is the critical density at the laser wavelength ( is the vacuum permittivity, the electron mass, and e the elementary charge). Electron and ion species are represented by 256 macroparticles per cell, with a fourth-order shape-function. All simulations describe field and electron impact induced ionization. Boundary conditions for both particles and fields are open along both longitudinal (x) and transverse (y) directions. The integration time of the simulations is after the laser peak.
To illustrate the coupling of the prepulse-free laser (DPM case) with two ultrathin Si foils as a function of their thickness, we display in
Figure 10 the spatial distribution of the transverse electric field (
). We observe near-complete laser transmission across the 10 nm target (
Figure 10a) while the 20 nm target gives rise to partial transmission and reflection (
Figure 10b). This behavior aligns with our experimental findings, i.e., high laser transmission across the 10 nm foil (
Figure 3c2) and no detectable protons for targets thinner than 20 nm (
Figure 5).
The simulated proton energy spectra, for a 200 nm Si + 50 nm Al target (DPM case, red dashed line) and a
Al target (no-DPM case, blue solid line), are plotted alongside the experimental data in
Figure 4a. Both simulation curves, integrated over the entire simulation duration (500 fs after peak laser impact), capture well the exponential shape of the distributions as well as their cutoff energies (
and 24 MeV, respectively).
Moreover, the simulated electron spectra, plotted in
Figure 4b, indicate that the no-DPM configuration generates hot electrons with energies similar to the DPM case, but in much larger numbers. However, these more energetic electron spectrum produced in the no-DPM case does not translate into better ion acceleration, as seen from
Figure 4a. This paradox indicates that the stronger ion acceleration in the DPM case has a different origin, e.g. an enhanced
field due to sharper density gradients at the time of the laser peak, as suggested in the following.
To explain the stronger proton acceleration achieved at high laser contrast, we plot in
Figure 11 the evolution of longitudinal lineouts (across the center of the focal spot) of the particle (electron and proton) densities and electric (
and
) fields, respectively measured in units of
and
. The time origin coincides with the on-target intensity peak. For a 200 nm Si + 50 nm Al target hit by the high-contrast pulse (top row), the accelerating
field seen by the protons (the front of which is located at
) reaches a maximum strength of
just after the laser maximum (
). This field decays over time (down to
at
), but its spatial maximum always coincides with the outer edge of the proton bunch (
Figure 11b,c), where the fastest protons reside (see
Figure 12b). This spatial overlap of the peak electric field and the fastest protons maximizes the acceleration efficiency.
The situation is markedly different in the absence of a DPM, for a 1.5
thick Al target (bottom row). As shown in
Figure 11d, the prepulse-driven expansion of the target backside causes a significant reduction in the
field [
58] just after the laser peak (
). Another consequence of the pre-expanded density profile is that the spatial maximum of
no longer coincides with the front of the proton bunch (see also in the phase space distribution in
Figure 12a).
Another interesting point to notice is the significant
field at
after the main pulse irradiation in the low laser contrast case (
Figure 11e and f). On the contrary, in the high laser contrast case (
Figure 11b), we see almost no
at the location of the protons. Importantly, this field is mostly positive, and hence will bend the protons towards the
direction, i.e., towards the laser propagation axis [
59]. This should contribute to the degraded beam profile seen in
Figure 6a, where the high-energy protons show an elongated shape towards the laser axis. Note that here the
field is not laser-cycle-averaged.
5. Conclusions
In summary, we demonstrate significant improvement in proton acceleration from thin solid targets using a double plasma mirror (DPM) [
34] to enhance the temporal contrast of the ultrashort (24 fs), ultraintense (
) F2 Apollon laser beam. This improvement is observed simultaneously in terms of maximum energy, conversion efficiency, and beam quality. The DPM mitigates the impact of the
long, relativistic-intensity (
contrast) laser pulse preceding the main pulse, thereby enabling solid foils as thin as 20 nm to be used. In contrast, the minimum target thickness exploitable under standard irradiation conditions is 800 nm. With the DPM, we were able to accelerate protons up to 15 MeV from a 20 nm thick foil, and achieved a peak proton energy of 35 MeV from a 250 nm foil, with only 5 J of laser energy. Remarkably, this maximum proton energy is
higher than without the DPM, despite 50% less laser energy. 2D PIC simulations attribute the superior performance of the DPM configuration to a stronger longitudinal electrostatic field and its closer overlap with the expanding proton bunch when the main laser pulse interacts with the sharp-gradient solid foil. These results pave the way for optimizing laser-driven X-ray and neutron sources [
11], opening new perspectives in nuclear material interrogation [
60,
61,
62,
63] and medicine [
7]. Furthermore, we can expect even further improvements in laser-proton acceleration at higher intensities, such as with the upcoming 10 PW-scale F1 beam of the Apollon facility.
Author Contributions
Conceptualization, J.F.; methodology, J.F., F.M. and A.A..; software, W.Y., Y.Y., E.d’H., X.D. and L.G.; investigation, R.L., T.W., J.F., W.Y., I.C., A.B., E.C., C.G., I.P., D.P., D.M. and L.G.; data curation, P.K., F.T., W.Y., J.F., Q.D., E.F., S.P. and R.L.; writing—original draft preparation, W.Y., J.F., R.L., E.F. and L.G.; writing—review and editing, all. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Experimental setup. The DPM was positioned between the off-axis parabola (OAP) and the target. The laser beam hit the target at 45 relative to the target normal. Cameras outside the chamber captured images of the reflected and transmitted laser beams on Spectralon plates. A stack of RCFs (with aluminum filters) was used to measure the proton energy spectrum and spatial profile.
Figure 1.
Experimental setup. The DPM was positioned between the off-axis parabola (OAP) and the target. The laser beam hit the target at 45 relative to the target normal. Cameras outside the chamber captured images of the reflected and transmitted laser beams on Spectralon plates. A stack of RCFs (with aluminum filters) was used to measure the proton energy spectrum and spatial profile.
Figure 2.
(a) Laser focal spot measurement (at low power) on a DPM shot. The intensity is calculated according to the spatial integration of the energy in each pixel, assuming that the full laser energy is contained in the image and taking into account the laser pulse duration. (b) Encircled laser energy distribution.
Figure 2.
(a) Laser focal spot measurement (at low power) on a DPM shot. The intensity is calculated according to the spatial integration of the energy in each pixel, assuming that the full laser energy is contained in the image and taking into account the laser pulse duration. (b) Encircled laser energy distribution.
Figure 3.
(a) Laser contrast measurement during direct (DPM-free) shots. (b) Images of the specularly reflected laser beam collected on a Spectralon plate during a direct shot on a thick Al target and (b1) a DPM shot on a 150 nm (100 nm Si + 50 nm Al) thick target (b2). The shadows present on the images come from objects placed between the target and the Spectralon plate. (c) Images of the transmitted laser beam on another Spectralon plate during DPM shots: (c1) calibration shot (without target), (c2) shot on a 10 nm Si foil, and (c3) shot on a 50 nm Si foil. Colormaps are normalized to the maximum count in the calibrated shot shown in (c1).
Figure 3.
(a) Laser contrast measurement during direct (DPM-free) shots. (b) Images of the specularly reflected laser beam collected on a Spectralon plate during a direct shot on a thick Al target and (b1) a DPM shot on a 150 nm (100 nm Si + 50 nm Al) thick target (b2). The shadows present on the images come from objects placed between the target and the Spectralon plate. (c) Images of the transmitted laser beam on another Spectralon plate during DPM shots: (c1) calibration shot (without target), (c2) shot on a 10 nm Si foil, and (c3) shot on a 50 nm Si foil. Colormaps are normalized to the maximum count in the calibrated shot shown in (c1).
Figure 4.
(a) Measured proton spectra, extracted from RCF data in the DPM case with a 250 nm (200 nm Si + 50 nm Al) thick target (red points, averaged over 3 shots) and in the no-DPM case with a thick Al target (blue points, averaged over 2 shots). Error bars represent the minimum and maximum values from different shots with the same target thickness. Exponential fits are shown as black dotted lines (DPM case) and purple dashed-dotted lines (no-DPM case). 2D PIC simulation results for both cases are plotted as red dashed and blue solid lines, respectively. (b) Simulated electron spectra together with exponential fits in the 4.5–11.5 MeV range. The line styles are the same as in panel (a).
Figure 4.
(a) Measured proton spectra, extracted from RCF data in the DPM case with a 250 nm (200 nm Si + 50 nm Al) thick target (red points, averaged over 3 shots) and in the no-DPM case with a thick Al target (blue points, averaged over 2 shots). Error bars represent the minimum and maximum values from different shots with the same target thickness. Exponential fits are shown as black dotted lines (DPM case) and purple dashed-dotted lines (no-DPM case). 2D PIC simulation results for both cases are plotted as red dashed and blue solid lines, respectively. (b) Simulated electron spectra together with exponential fits in the 4.5–11.5 MeV range. The line styles are the same as in panel (a).
Figure 5.
Target thickness dependence of the (a) proton cutoff energy and (b) laser-to-proton energy conversion efficiency (see text for definition). For both quantities, we distinguish between the maximum (empty symbols) and mean (filled symbols) values for all shots performed with the same target thickness. Direct and DPM shots are represented by squares and dots, respectively. The target thickness for DPM shots refers to the entire Si/Al target. For the thinnest targets, the results are single-shot, hence no mean is reported.
Figure 5.
Target thickness dependence of the (a) proton cutoff energy and (b) laser-to-proton energy conversion efficiency (see text for definition). For both quantities, we distinguish between the maximum (empty symbols) and mean (filled symbols) values for all shots performed with the same target thickness. Direct and DPM shots are represented by squares and dots, respectively. The target thickness for DPM shots refers to the entire Si/Al target. For the thinnest targets, the results are single-shot, hence no mean is reported.
Figure 6.
Spatial dose distribution on the RCF stack. Top row: no-DPM shot on a 800 nm thick Al target. Middle row: DPM shot on a 150 nm (100 nm+50 nm Al) thick target. Bottom row: DPM shot on a 20 nm thick Si target. The first two shots yield similar maximum proton energies, but different spatial beam profiles. The same colormap is used in all images. The arrows locate the laser propagation direction (red) and the target normal (black) direction.
Figure 6.
Spatial dose distribution on the RCF stack. Top row: no-DPM shot on a 800 nm thick Al target. Middle row: DPM shot on a 150 nm (100 nm+50 nm Al) thick target. Bottom row: DPM shot on a 20 nm thick Si target. The first two shots yield similar maximum proton energies, but different spatial beam profiles. The same colormap is used in all images. The arrows locate the laser propagation direction (red) and the target normal (black) direction.
Figure 7.
Measured half-width of the proton angular distribution as a function of proton energy (normalized to the cutoff energy). Filled blue squares: direct shots on 800 nm thick Al foils (symbols represent the average values and error bars the maximum and minimum values over the considered series of shots). Filled red circles: DPM shots on 150 nm (100 nm Si+50 nm Al) thick foils. Empty red circles: DPM shots on 20 nm thick Si foils. The black dashed curve represents the energy-dependent angular distribution observed in many experiments to be characteristic of TNSA protons [
21].
Figure 7.
Measured half-width of the proton angular distribution as a function of proton energy (normalized to the cutoff energy). Filled blue squares: direct shots on 800 nm thick Al foils (symbols represent the average values and error bars the maximum and minimum values over the considered series of shots). Filled red circles: DPM shots on 150 nm (100 nm Si+50 nm Al) thick foils. Empty red circles: DPM shots on 20 nm thick Si foils. The black dashed curve represents the energy-dependent angular distribution observed in many experiments to be characteristic of TNSA protons [
21].
Figure 8.
Spatial dose distribution recorded inside the interaction chamber by a set of dosimeters (details in text). Left panel: results for direct shots on Al foils of , and thicknesses. Right panel: results for DPM shots on Si/Al foils of 150 nm, 200 nm and 250 nm thicknesses. Each circle represents a dosimeter, and its color corresponds to the measured dose. The outer circle in both graphs represents the actual diameter of the target chamber. The grey bar representing the DPM in the right panel is to scale.
Figure 8.
Spatial dose distribution recorded inside the interaction chamber by a set of dosimeters (details in text). Left panel: results for direct shots on Al foils of , and thicknesses. Right panel: results for DPM shots on Si/Al foils of 150 nm, 200 nm and 250 nm thicknesses. Each circle represents a dosimeter, and its color corresponds to the measured dose. The outer circle in both graphs represents the actual diameter of the target chamber. The grey bar representing the DPM in the right panel is to scale.
Figure 9.
X-ray spectra averaged over seven DPM shots on a 200 nm Si foils with a 50 nm Al coating (on the laser-irradiated side) (black curve). The red and blue curves represent flychk atomic physics calculations for Al emission at and , and and , respectively.
Figure 9.
X-ray spectra averaged over seven DPM shots on a 200 nm Si foils with a 50 nm Al coating (on the laser-irradiated side) (black curve). The red and blue curves represent flychk atomic physics calculations for Al emission at and , and and , respectively.
Figure 10.
Spatial distribution of the transverse electric field (in units), in the DPM (no-prepulse) laser configuration with a Si target of (a) 10 nm and (b) 20 nm thickness. The black dotted box shows the target initial position.
Figure 10.
Spatial distribution of the transverse electric field (in units), in the DPM (no-prepulse) laser configuration with a Si target of (a) 10 nm and (b) 20 nm thickness. The black dotted box shows the target initial position.
Figure 11.
Simulations of a DPM shot on a 200 nm Si + 50 nm Al target (top row) and a direct shot on a 1.5 µm thick Al target (bottom row). Longitudinal cuts (across the center of the focal spot at ) of the electron density (red solid curves), proton density (blue dotted curves), electric field (cyan dashed curves) and (averaged over the laser cycle, purple dotted-dashed curves) are plotted at various times after the laser peak has struck the target (). Particle densities are measured in units of while field strengths are shown in units of . The target rear side is initially located at with the DPM and without it.
Figure 11.
Simulations of a DPM shot on a 200 nm Si + 50 nm Al target (top row) and a direct shot on a 1.5 µm thick Al target (bottom row). Longitudinal cuts (across the center of the focal spot at ) of the electron density (red solid curves), proton density (blue dotted curves), electric field (cyan dashed curves) and (averaged over the laser cycle, purple dotted-dashed curves) are plotted at various times after the laser peak has struck the target (). Particle densities are measured in units of while field strengths are shown in units of . The target rear side is initially located at with the DPM and without it.
Figure 12.
Phase space distribution () of protons of (a) the direct shot on a 1.5 µm thick Al target and (b) the DPM shot on a 200 nm Si + 50 nm Al target.
Figure 12.
Phase space distribution () of protons of (a) the direct shot on a 1.5 µm thick Al target and (b) the DPM shot on a 200 nm Si + 50 nm Al target.