1. Introduction
High intensity laser-driven ion acceleration has been studied for more than two decades and nowadays short proton bunches of up to tens of megaelectronvolts with $10^{12}$
particle/sr/shot and divergences of the order of several degrees are measured (Snavely et al. Reference Snavely2000; Steinke et al. Reference Steinke2020). These properties are obtained when the acceleration mechanism at play is Target Normal Sheath Acceleration (TNSA) (Passoni, Bertagna & Zani Reference Passoni, Bertagna and Zani2010) where an ultra-high-intensity laser interacts with a solid target. It creates a hot electron population that leaves the target generating fields capable of extracting ions from the surface of the target to accelerate them. Nevertheless using solid targets requires precise alignment at each shot and leads to the generation of debris and the coating of surrounding high-intensity laser optics rendering it delicate to use on the new high repetition rate ultra-high-intensity laser facilities such as Apollon (Papadopoulos et al. Reference Papadopoulos2016), ELI (Margarone et al. Reference Margarone2018; Radier et al. Reference Radier2022) or VEGA (Huault et al. Reference Huault2017).
To mitigate this issue, self-regenerating targets are under study: comprising of water jets (George et al. Reference George2019; Puyuelo-Valdes et al. Reference Puyuelo-Valdes2022) or droplets (Karsch et al. Reference Karsch, Düsterer, Schwoerer, Ewald, Habs, Hegelich, Pretzler, Pukhov, Witte and Sauerbrey2003; Ter-Avetisyan et al. Reference Ter-Avetisyan, Schnürer, Nickles, Smirnov, Sandner, Andreev, Platonov, Psikal and Tikhonchuk2008), cryogenic ribbons (Margarone et al. Reference Margarone2016) or liquid crystal films (Thoss et al. Reference Thoss, Richardson, Korn, Faubel, Stiel, Vogt and Elsaesser2003; Poole et al. Reference Poole, Andereck, Schumacher, Daskalova, Feister, George, Willis, Akli and Chowdhury2014). Near-critical density targets such as high-density gas jets (Palmer et al. Reference Palmer2011; Haberberger et al. Reference Haberberger, Tochitsky, Fiuza, Gong, Fonseca, Silva, Mori and Joshi2012; Henares et al. Reference Henares2019), potentially cryogenically cooled (Obst et al. Reference Obst2017; Rehwald et al. Reference Rehwald2023) are also under investigation. Using such targets enables other acceleration mechanisms to participate in the particle energisation. Among the regimes unlocked in this configuration are Magnetic Vortex Acceleration (Bulanov et al. Reference Bulanov, Dylov, Esirkepov, Kamenets and Sokolov2005; Nakamura et al. Reference Nakamura, Bulanov, Esirkepov and Kando2010; Helle et al. Reference Helle, Gordon, Kaganovich, Chen, Palastro and Ting2016) or Collisionless Shock Acceleration (CSA) (Silva et al. Reference Silva, Marti, Davies, Fonseca, Ren, Tsung and Mori2004; Fiuza et al. Reference Fiuza, Stockem, Boella, Fonseca, Silva, Haberberger, Tochitsky, Mori and Joshi2013; Wan et al. Reference Wan, Andriyash, Hua, Pai, Lu, Mori, Joshi and Malka2019). In the latter, the laser pulse generates an electrostatic shock leading to the formation of an interface between two plasma slabs of significantly different densities. This shock propagates colinearly with the laser axis in the plasma at the speed $v_{sh}= M C_s$
where $M$
is its Mach number and $C_s$
the ion acoustic velocity. Above a certain critical Mach number, this shock is able to accelerate ions up to twice its speed. In these idealised simulations, the target is a thin near-critical density slab with very sharp density gradients. This acceleration process can also be studied using pre-expanded foils (Gauthier et al. Reference Gauthier2014; Pak et al. Reference Pak2018; Tochitsky et al. Reference Tochitsky, Pak, Fiuza, Haberberger, Lemos, Link, Froula and Joshi2020). Coulomb explosion has also been reported as the responsible mechanism for ion acceleration in underdense targets (Krushelnick et al. Reference Krushelnick1999; Sarkisov et al. Reference Sarkisov, Bychenkov, Novikov, Tikhonchuk, Maksimchuk, Chen, Wagner, Mourou and Umstadter1999) as well as in the near-critical regime (Wei et al. Reference Wei2004). Here, some electrons, heated by the laser pulse, escape the target volume and leave a space charge force behind them responsible for ion energisation. Using near-critical density gas targets, the hole boring acceleration mechanism can also participate in the acceleration of ions. Here, the radiation pressure associated with the laser pulse is intense enough to push an overdense layer inwards towards the target and accelerate the particles within (Wilks et al. Reference Wilks, Kruer, Tabak and Langdon1992; Pukhov & Meyer-ter Vehn Reference Pukhov and Meyer-ter Vehn1997).
The main drawback associated with high density gas jets available nowadays comes from the need to obtain a thin layer (${\approx }10$
s of $\mathrm {\mu }$
m) of near-critical density (${\approx }10^{21}$
e cm$^{-3}$
for 1 $\mathrm {\mu }$
m laser wavelength) at the centre of the target that is positioned in vacuum. Practically, current nozzle designs do not enable the production of a target displaying these ideal characteristics. In fact relatively long (100s of $\mathrm {\mu }$
m), low density (${\approx }10^{19}$
e cm$^{-3}$
) regions are obtained in which the laser has to propagate before reaching the centre of the target and driving the shock (Henares et al. Reference Henares2019) in a region that can be thicker than the optimal value for ion acceleration. The density profile obtained with such targets is presented in figure 1(a), where the solid blue line corresponds to experimental data and the dashed red line represents the result of a Computational Fluid Dynamics (CFD) simulation. It is worth noting that there is no adjusting parameter to ensure the very good agreement observed between these two profiles. Moreover, changing the backing pressure will only result in increasing the density overall and will not modify the shape of the profile.
Figure 1. Density profile 500 $\mathrm {\mu }$
m from the nozzle fed with 100 bar backing pressure of N$_2$
exit retrieved with an interferometry set up (solid blue line) or with CFD simulations (red dashed line) (a). Figure extracted from Henares et al. (Reference Henares2019). Schematic side-view of the laser arrangement in the longitudinal (b) and transverse tailoring (c) studies, respectively.
When interacting with a high intensity laser pulse, initially the nanosecond Amplified Spontaneous Emission (ASE) part of the pulse propagates in the target and modifies the initial plasma density profile through optical tailoring. Then, during the propagation in the underdense region of this modified target, the laser undergoes nonlinear effects such as filamentation, loses its integrity and part of its energy, and hence drives a shock that is not as intense nor necessarily in the forward direction of the laser pulse. Consequently, the properties of the accelerated particles are modified, as was shown by Puyuelo-Valdes et al. (Reference Puyuelo-Valdes2019).
Precisely tuning the ASE of a laser chain is challenging, hence optically tailoring the gas target with additional, dedicated prepulses can help improve the control over the obtained density profile and optimise the interaction. This has been experimentally studied with CO$_2$
lasers ($\lambda _L = 10\ \mathrm {\mu }$
m) by Tresca et al. (Reference Tresca, Dover, Cook, Maharjan, Polyanskiy, Najmudin, Shkolnikov and Pogorelsky2015) and Dover et al. (Reference Dover, Cook, Tresca, Ettlinger, Maharjan, Polyanskiy, Shkolnikov, Pogorelsky and Najmudin2016) with a copropagating laser pulse to create the sharpest gradient, hence driving the most efficient shock to accelerate ions. Promising results were also obtained earlier by Haberberger et al. (Reference Haberberger, Tochitsky, Fiuza, Gong, Fonseca, Silva, Mori and Joshi2012) with a train of pulses sent to a near-critical density target with the same laser wavelength. More recently another scheme was discussed in the 1 $\mathrm {\mu }$
m laser wavelength regime with tailored laser pulses propagating transversely to the driving laser pulse (Bonvalet et al. Reference Bonvalet2021; Marquès et al. Reference Marquès2021). In this scheme the two laser pulses will drive blast waves in front and behind the gas target (with respect to the main laser propagation axis) that will heat and compress the initial density profile to form a sharp gradient, thin plasma slab. The present collaboration led recent experimental campaigns on 1 $\mathrm {\mu }$
m laser facilities where the copropagating scheme, that will be called longitudinal, and the transverse technique were studied. Both schemes are schematically illustrated in figures 1(b) and 1(c), respectively. A publication focused on a more detailed study on the results obtained in the transverse tailoring configuration was published very recently (Marquès et al. Reference Marquès2024).
In this manuscript, the experimental set-up associated with both campaigns will be described and the main results will be presented. Finally, comparisons will be made between the two schemes and conclusions will be drawn.
2. Experimental set-up
2.1. Longitudinal tailoring
This study was conducted at GSI with the Petawatt High-Energy Laser for heavy Ion eXperiments (PHELIX) (Bagnoud et al. Reference Bagnoud2010). The prepulse and the main pulse were issued from the same front end, hence had the same characteristics, as follows: the wavelength was 1054 nm with a pulse duration of 500 fs and a contrast of $10^{11}$
80 ps before and 10$^{10}$
100 ps after the main pulse. It was focused in the target chamber by an $f/1.7$
off-axis parabola down to a focal spot size of 5.3 $\mathrm {\mu }$
m full width at half maximum (FWHM) with ${\approx }40\,\%$
of the energy contained in the central spot. The root-mean-square (r.m.s.) transverse pointing fluctuation was of the order of 15 % of the FWHM focal spot size in both directions. For ion acceleration studies, the delay between the two pulses was scanned between 10 ps and 1 ns while the energy in the prepulse was varied from 0.1 to 50 J. This results in intensities delivered by the prepulse as high as $6.9 \times 10^{19}$
W cm$^{-2}$
. The energy in the main pulse was varied from 50 to 100 J corresponding to intensities up to $1.5 \times 10^{20}$
W cm$^{-2}$
. The r.m.s. energy fluctuation was 10 %. These two spatially overlapping beams are represented in red and orange in figure 2.
Figure 2. Schematic top-view of the experimental set-up implemented for the longitudinal tailoring study. In red/orange the prepulse and main beam are displayed (respectively) while the probe beam path is in green and the three magnetic spectrometers (MS) deployed are indicated in black. The 90$^\circ$
MS is located below the probe beam path. The schematic is not to scale.
The pulses were focused on a high-density gas jet target previously described by Henares et al. (Reference Henares2019) made of a supersonic de Laval nozzle showing an exit diameter of 240 $\mathrm {\mu }$
m, fed by a pressure of up to 1000 bar of hydrogen. The target was precisely aligned thanks to side-view and bottom-view optical diagnostics not depicted in figure 2. The interaction could be probed thanks to a pick off of the beams subsequently frequency doubled to 526 nm and propagating through the target. This 500 fs probe beam could be delayed thanks to a motorised delay stage positioned in the chamber. This beam is depicted in green in figure 2. A shadowgraphy and Nomarski interferometry set-ups were developed to diagnose the generated plasma at different interaction times. These diagnostics were protected with a series of filters including strong bandpass filters (FB530-10 from Thorlabs).
Three Magnetic Spectrometer (MS, represented in black in figure 2) coupled to image plate (IP) were used to detect the accelerated protons. They were positioned at 0$^\circ$
, 38$^\circ$
and 90$^\circ$
from the target with respect to the laser axis. They were shielded with lead to increase the signal-to-noise ratio and equipped with pinholes to improve the energy resolution (at distances of 51.5, 165 and 7.5 cm from the interaction point corresponding to 0.15, 0.03 and 1.4 $\mathrm {\mu }$
sr acceptance solid angle, respectively). The energy resolution of these MS was 0.5 % at 1 MeV and 1.2 % at 6 MeV at 0 and 90$^\circ$
and 2 % at 1 MeV and 6 % at 6 MeV at 38$^\circ$
. The BAS-MS IPs from the Fuji Photo-Film Co. Ltd were used and analysed with a FUJIFILM FLA-7000 reader. The IP response function to protons has been taken from Bonnet et al. (Reference Bonnet, Comet, Denis-Petit, Gobet, Hannachi, Tarisien, Versteegen and Aleonard2013).
2.2. Transverse tailoring
The study on transverse tailoring discussed in the present manuscript took place at the LULI laboratory (Ecole Polytechnique, France) on the PICO2000 laser. In this experiment the tailoring pulses propagated transversely compared with the driving pulse. The latter had a wavelength of 1054 nm, a pulse duration of 1 ps and a contrast of $10^9$
80 ps before and 100 ps after the main pulse. It was focused by an off-axis parabola to a focal spot size of 10 $\mathrm {\mu }$
m FWHM (displayed in orange in figure 3). The energy in this pulse was varied from 40 to 60 J resulting in an intensity as high as $2 \times 10^{19}$
W cm$^{-2}$
. The tailoring pulses were also at a wavelength of 1054 nm but presented a duration of 0.6 ns and were focused thanks to a lens down to 42 $\mathrm {\mu }$
m FHWM. The energy in each arm varied from 4 to 6 J. The focal spots positions for all beams varied by $\pm \sim 20\ \mathrm {\mu }$
m over the experimental run. These two tailoring beams are colinear and represented in red in figure 3.
Figure 3. Schematic top-view of the experimental set-up implemented for the transverse tailoring study. In orange the main beam path is displayed while the nanosecond tailoring beam path is in red prior to the target, the probe beam is in green and the two MS deployed are indicated in black. The schematic is not to scale.
The target used is similar to the one presented in the previous section. The probe beam came from a seeded Nd:YAG laser source which delivers 7 ns pulses at 532 nm (Quanta-Ray LAB-150). This beam is depicted in green in figure 3. The optical diagnostics used in this experiment consisted of a shadowgraphy set-up with a CCD camera coupled to a gated optical intensifier (GOI) recording two-dimensional spatially resolved snap shots and a streak camera to monitor the time evolution along the picosecond axis.
Similarly to the previous set-up, two MS with IPs were used to detect the accelerated protons (represented in black in figure 3). They were positioned at 0$^\circ$
and 30$^\circ$
from the target with respect to the laser axis (at distances of 55 and 45 cm from the interaction point corresponding to 0.65 and 0.97 $\mathrm {\mu }$
sr acceptance solid angle, respectively). The energy resolution of these MS was 0.2 % at 1 MeV. The same acquisition and analysis procedures have been applied as described in the previous section.
3. Experimental results
3.1. Longitudinal tailoring
3.1.1. Plasma target tailoring study
Initially, studies were conducted to investigate the formation of a plasma channel at low laser energy and low hydrogen backing pressures. Thanks to the interferometry diagnostic, it was revealed that a lower density plasma channel is formed, as can be seen in figure 4(a). This image results from firing a 5 J laser pulse corresponding to an intensity of $5.8 \times 10^{18}$
W cm$^{-2}$
on 100 bar hydrogen backing pressure target probed 312 ps after the interaction. On this image the laser propagates from left to right and its focusing geometry is depicted by the red cone. The nozzle is placed on top, not shown on the image. The fringes’ curvature appears to change twice around the laser axis (vertically on the image), which is a signature of a lower density region. Analysing this image to obtain the plasma density is challenging but focusing on a region where the fringes are relatively symmetric around the laser axis, $Z = -300\ \mathrm {\mu }$
m and $Y = -200\ \mathrm {\mu }$
m depicted as region A in figure 4(a), allows us to estimate the value of the density. This is what is depicted in figure 4(b) where the solid blue line is the density profile obtained along the white dashed line in figure 4(a) and the red dashed line is the results of a CFD simulation without a tailoring pulse. In this area (edges of the plasma channel), the analysis showed an increase of the density by a factor $\sim$
2.5 compared with the expected density from the unperturbed gas jet obtained from CFD simulations. Assuming that this density has been pushed by the laser pulse away from the laser axis, where the plasma channel is formed (region B in figure 4a), this means that at this delay the electronic density in the plasma channel is lower than without plasma tailoring.
Figure 4. Raw interferogram of the target firing at 100 bar hydrogen backing pressure 312 ps after the interaction with a 5 J laser pulse corresponding to an intensity of $5.8 \times 10^{18}$
W cm$^{-2}$
. The nozzle is on top of the image (not shown) and the laser propagates from left to right (a). The 0 graduations on the longitudinal (Z) and height (Y) axis mark where the laser is focused. The focusing laser is indicated by the red cone. The letters A and B correspond to regions that will be referred to in the core of the text. The white dashed line is where the density profile is reconstructed and shown in (b) with a solid blue line. The dashed red line in (b) corresponds to CFD simulation results without a tailoring pulse.
For laser–plasma interaction at higher backing pressures, radiation will be emitted at twice the initial frequency of the laser pulse in the regions where it interacts with a high plasma density and sharp gradient. This is the bright emission region that can be seen in figure 5. The figure shows images recorded on the shadowgraphy diagnostic with the nozzle on top and the focusing laser represented by the red cone. The vertical line represents the centre of the nozzle where the density is the highest. For both these shots the backing pressure was fixed at 800 bar of hydrogen and the probing time was 312 ps after the interaction with the main pulse. Figure 5(a) is obtained without prepulse and with 5 J in the main pulse ($5.6 \times 10^{18}$
W cm$^{-2}$
) while figure 5(b) is with a 1 J prepulse ($1.5 \times 10^{18}$
W cm$^{-2}$
) sent 300 ps before the 4 J main pulse ($4.5 \times 10^{18}$
W cm$^{-2}$
). The dark region below the trapezoidal nozzle shape corresponds to a region significantly larger than the volume where the density is above the critical density for the 532 nm probe pulse wavelength because of its refraction in the density gradients. This is described in detail in Marquès et al. (Reference Marquès2024). It is clear that with the prepulse, the frequency-doubled radiation is emitted after the centre of the gas jet (${\approx }100\ \mathrm {\mu }$
m) while without it is emitted ${\approx }100\ \mathrm {\mu }$
m prior to the centre. Measuring the $2\omega$
radiation being emitted from after the centre of the nozzle indicated that the laser propagated farther inside the high-density part of the target and delivered its energy in this region, which was in good agreement with the initial goal.
Figure 5. Shadowgraphy images obtained after the interaction of a laser pulse with a gas jet fed by a 800 bar hydrogen backing pressure. The probe beam is sent 312 ps after the interaction with the main pulse. The focusing laser pulse is depicted with a red cone while the vertical line represents the centre of the nozzle. Here (a) is obtained without a prepulse and with a main pulse intensity of a $5.8 \times 10^{18}$
W cm$^{-2}$
while (b) is with a 1 J prepulse ($1.5 \times 10^{18}$
W cm$^{-2}$
) sent 300 ps before the main pulse for which the intensity was $4.5 \times 10^{18}$
W cm$^{-2}$
.
3.1.2. Particle acceleration
Since the previous section showed the formation of a plasma channel that enabled the laser pulse to deposit its energy farther inside the gas target, the next step was to look at the effect of this optical tailoring on the particle acceleration thanks to the three MS installed.
To do so the laser energy of the main pulse was increased up to 100 J (corresponding to an intensity of $1.5 \times 10^{20}$
W cm$^{-2}$
) and the hydrogen backing pressure to 1000 bar. The energy contained in the prepulse was varied between 0.1 and 1 J (corresponding to intensities of $1.5 \times 10^{17}$
W cm$^{-2}$
and $1.5 \times 10^{18}$
W cm$^{-2}$
, respectively) and the delay between the prepulse and the main pulse was tuned between 10 and 1000 ps. A reference shot without prepulse was also recorded for comparison. The results are presented in figure 6 where figure 6(a–c) correspond to the data obtained with a 0.1 J prepulse on the different MS and figure 6(d–f) to the series obtained with a 1 J prepulse. The rows correspond to the MS at 0$^\circ$
(figure 6a,d), 38$^\circ$
(figure 6b,e) and 90$^\circ$
(figure 6d,f), respectively. The different solid lines correspond to the distinct prepulse delays used while the dashed lines represent the detection limits.
Figure 6. Series of spectra obtained during the prepulse energy and delay scans. Panels (a–c) corresponds to 0.1 J energy in the prepulse while the second one (d)–(f) corresponds to 1 J energy (intensities of $1.5 \times 10^{17}$
W cm$^{-2}$
and $1.5 \times 10^{18}$
W cm$^{-2}$
, respectively). The three rows correspond to the 0$^\circ$
(a,d), 38$^\circ$
(b,e) and 90$^\circ$
(c,f) MS, respectively. The different spectra in the same graphic, shown in different colours, correspond to distinct delays between the prepulse and the arrival of the main pulse. The black solid line indicated the shot without the prepulse. The detection limit is calculated with the noise specific to each shot and indicated by dashed lines.
At 0$^\circ$
both the number of accelerated particles and the maximum energy obtained reduce when adding a prepulse in the delay range explored similarly with a 0.1 or 1 J prepulse (see figures 6a and 6d). It seems to follow a trend where the longer the delay, the less protons are accelerated at lower energies. The effect is stronger at higher prepulse energy.
At 38$^\circ$
structures appear at relatively low delay (10 to 100 ps) with energy peaks of 6 MeV but tend to vanish at larger delays for both prepulse energy. There are two clear structures at a delay of 10 ps with a 0.1 J prepulse at ${\approx }3.5$
and ${\approx }6$
MeV persisting at 100 ps delay (figure 6b). While for the 1 J case a structure (at $\approx$
2.5 MeV) can only be seen for a 100 ps delay (figure 6e).
At 90$^\circ$
, no clear structures can be observed in the range between 1 and 6 MeV. The proton signal detected at 1 J prepulse energy and 10 ps delay is very perturbed and requires more work to be presented. Moreover, the energy limit is set at 6 MeV on these spectra because the spectrometer was positioned much closer to the interaction point ($\sim$
7 cm) than to the other two (>50 cm), leading to a very large x ray spot on the IP during some shots. The result was a very high level of background noise, similar in intensity to the proton signal, which renders the estimation of the proton number challenging above 6 MeV. Nevertheless, in the configuration where a 0.1 J prepulse was fired 100 ps before the main pulse, we were able to observe a structure around $\approx$
7.5 MeV clearly separated from the background, as shown in the raw IP presented in figure 7.
Figure 7. Raw image of an IP positioned in the 90$^\circ$
MS after a 0.1 J prepulse was fired 100 ps before the main pulse at 1000 bar backing pressure. A clear structure can be observed $\approx 7.5$
MeV.
The observed structures at different angles cannot be attributed to TNSA or Coulomb explosion. The latter is most likely responsible for the acceleration of the low energy, continuous proton distribution on axis and at the other angles as the shape and maximum energy observed are compatible with the results reported in Krushelnick et al. (Reference Krushelnick1999) and Sarkisov et al. (Reference Sarkisov, Bychenkov, Novikov, Tikhonchuk, Maksimchuk, Chen, Wagner, Mourou and Umstadter1999). In fact, the maximum energy is the relativistic ponderomotive energy $U_p = m_e c^2(\gamma -1)$
where $\gamma = (1+a_0^2/2)^{(1/2)}$
is the Lorentz factor for the electron transverse motion. Given the $a_0$
(normalised laser intensity) of the main pulse reached during that experiment, it corresponds to $\approx$
3 MeV which is in good agreement with the measured energy. The observed structures could be the result of shocks driven by the filamenting laser pulse. Among the ion acceleration mechanisms accessible during the interaction of a laser and underdense or near-critical density plasmas, CSA is one of the few that can produce structured spectra. If the prepulse does not generate a plasma channel with a low enough density, because the energy is too small or the delay too short for example, the main laser pulse could still filament and drive shocks in various directions leading to the observations of structures in the spectra at distinct angles and for different delays. It is worth noting that no clear structures could be observed without a prepulse. This could be explained because a clear shock cannot easily propagate in the current experimental conditions without tailoring.
3.2. Transverse tailoring
3.2.1. Plasma target tailoring study
A reference shot was taken without tailoring pulses at 800 bar backing pressure with a main pulse intensity of $2 \times 10^{19}$
W cm$^{-2}$
with the GOI triggered 200 ps after the arrival of the main pulse, exposed for 100 ps and the shadowgraphy image associated with this shot is shown in figure 8(a). As in figure 5, the nozzle is on top of the figure and the focusing laser is depicted by the red cone. Similarly to the previous section, the dark region is significantly larger than the actual transverse size of the volume where the density is above the critical one. Figure 8(b) was obtained with two tailoring pulses arriving 2.4 ns before the main pulse with a 200 bar backing pressure and the GOI triggered 100 ps before the arrival of the main pulse and exposed for 100 ps. In that case, blast waves propagate towards the centre of the target and create a dense plasma slab in the centre where the driving laser is focused. The associated streak camera images are shown in figures 8(c) and 8(d) where the main laser pulse arrival time is indicated by a red arrow.
Figure 8. (a) Shadowgraphy data obtained after the interaction of a $1.2 \times 10^{19}$
W cm$^{-2}$
laser pulse with a gas jet fed by a 800 bar hydrogen backing pressure. (b) Similar data obtained at 200 bar backing pressure with two $3 \times 10^{14}$
W cm$^{-2}$
tailoring pulses sent 2.4 ns before the main pulse. The diagnostic was triggered 200 ps or 100 ps before the arrival of the main pulse, respectively, and in both cases exposed for 100 ps. The focusing laser pulse is depicted with a red cone while the vertical line represents the centre of the nozzle. (c,d) Corresponding laser shots on the streak camera where the main laser pulse arrival time is indicated by a red arrow.
Using the streak camera diagnostic, without the prepulses, it can be noticed that the main pulse creates almost instantaneously a 1 mm-wide plasma (figure 8c). With the tailoring beams (figure 8d), the propagation of the two blast waves can be diagnosed temporally until they collide and form the transient target. One can note that, on this shot, the main pulse arrives after the thinnest target size was obtained.
3.2.2. Particle acceleration
The spectra obtained on the magnetic spectrometers at 0$^\circ$
and 30$^\circ$
for the two shots described in the previous section are represented in figures 9(a) and 9(b), respectively. A continuous distribution is obtained with a maximum energy lower than 1.5 MeV at both angles without tailoring (black solid lines). With the tailoring pulses, structures are observed in the distributions around 1–1.5 MeV. These structures are a signature of a different acceleration mechanism compared with the case without tailoring. When the target is not optically modified, the spectrum obtained resembled the ones obtained from the Coulomb explosion mechanism with a characteristic continuous distribution at various angles. When tailoring the target, quasimonoenergetic structures appear that could be the result of CSA due to the presence of sharper gradients enabled by the tailoring of the target. It can also be noted that, overall, fewer protons are produced at 30$^\circ$
rendering the acceleration process more directional as expected from CSA and more protons are accelerated to relatively high energy with plasma tailoring.
Figure 9. Proton spectra obtained at 0$^\circ$
(a) and 30$^\circ$
(b) from the laser axis with (red solid line) and without (black solid line) the tailoring pulses as described in figure 8. The detection limit is calculated with the noise specific to each shot and indicated by dashed lines.
4. Discussions
Several challenges were identified during the experimental campaigns where high-density gas targets were used. First, due to the density profiles obtained with the nozzles available nowadays, optical tailoring strategies must be implemented experimentally to reduce the low-density wings and have thin, near-critical density targets. This optical tailoring is challenging to precisely diagnose in three dimensions as optical probing results in two-dimensional reconstruction that can be misleading. Furthermore, the beam time constraints and repetition rate on such laser facilities also limits the parameter phase space that can be explored in these multibeam experiments (energy and duration of the pulses, Rayleigh lengths, delays, wavelength, $\ldots$
). One can also note that spatially extended particle diagnostics could improve the obtained data as the angles at which the particles are accelerated are difficult to predict and could vary shot-to-shot. Target tailoring, as well as the shot-to-shot variation of the laser pulse focal spot, which are reasonably small in the presented work, could lead to a modification of the symmetry of the target, hence ions could be emitted at angles that are not covered by the particle detectors with the limited acceptance angle used here.
Nevertheless, a clear effect of the optical tailoring of the high-density gas jet targets on the particle spectra obtained has been demonstrated in both configurations. In the longitudinal one, the first observations seem to indicate that a plasma channel can be formed, enabling the laser pulse to propagate and deposit its energy farther towards the centre of the target. However, within the performed parameters scan, this does not result in the acceleration of higher energy protons in the forward direction, as could be expected. Still, clear structures could be observed up to energies of 6 MeV at 38$^\circ$
and 90$^\circ$
only when the target was optically tailored. This could be a signature of the laser pulse propagating to a higher density but still undergoing filamentation and generating shocks at various angles. During this experiment, the explored phase space was limited especially in terms of prepulse energy and delay between the pulses and future experiments are required to extend this phase space. It could also be noted that the profile at the rear of the target was not tailored in this scheme as opposed to the transverse one, which is not ideal for CSA to occur according to the particle-in-cell simulations presented by Fiuza et al. (Reference Fiuza, Stockem, Boella, Fonseca, Silva, Haberberger, Tochitsky, Mori and Joshi2013).
The transverse scheme is more complicated to set up as it requires the synchronisation and spatial alignment of several beams from distinct sources but showed promising results as clear structures were obtained in the forward direction. It is worth noting that this scheme facilitates working at much lower backing pressure which is in favour of using it at higher repetition rate. Nevertheless, the explored phase space was again limited.
To conclude, even though a clear effect associated with the plasma tailoring was measured, leading to the observation of structures in the spectrum at several megaelectronvolts of energy, each of the presented spectra were obtained from a single shot. To mitigate the potentially biased results and to enable the exploration of broader phase spaces, parametric studies would need to be performed on high repetition rate, high intensity laser facilities. This would enable the exploration of common phase spaces between the two configurations, hence making a direct comparison easier but also unlocking a clear understanding of the mechanisms at play to develop the laser plasma acceleration of ions on near-critical gas targets.
Acknowledgements
Part of the results presented here are based on the experiment P-21-00004, which was performed at the PHELIX facility at the GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt (Germany). The authors would like to give special thanks to the PHELIX operation team as well as the staff of the LULI facility.
Editor V. Malka thanks the referees for their advice in evaluating this article.
Declaration of interest
The authors report no conflict of interest.
Funding
This work is supported by IN2P3-CNRS (nos. 2016-2024 ALP-IONS projects); Agence Nationale de la Recherche (ANR-17-CE30- 0026-Pinnacle); the Region Nouvelle Aquitaine (AAPR2022-2021-17029010 VALERIA); LASERLAB-EUROPE (grant agreement no. 871124, European Union's Horizon 2020 research and innovation programme) and the authors acknowledge the financial support of the IdEx University of Bordeaux / Grand Research Program ‘GPR LIGHT’.
Data availability
The data that support the findings of this study are available from the corresponding authors, A.M. or M.T., upon reasonable request.
Open access
