2 Model and simulation parameters

To verify our scheme, we use the particle-in-cell (PIC) simulation code EPOCH^[ Reference Arber, Bennett, Brady, Lawrence-Douglas, Ramsay, Sircombe, Gillies, Evans, Schmitz, Bell and Ridgers^31]. The two-dimensional (2D) simulation box is $x\times y=100{\lambda}_0\times 120{\lambda}_0$ , with 100 cells per ${\lambda}_0$ and 36 macroparticles per cell, where ${\lambda}_0=1\, \unicode{x3bc} \mathrm{m}$ is the laser wavelength. The outer diameter and thickness of the hemisphere target are $20{\lambda}_0$ and $3{\lambda}_0$ , respectively. Its material depends on the desired ion beam. For simplicity, we consider a hydrogen plasma target of density ${n}_\mathrm{e}={n}_\mathrm{p}=100{n}_\mathrm{c}$ , where ${n}_\mathrm{c}={m}_\mathrm{e}{\omega}_0^2{\varepsilon}_0/{e}^2$ is the critical plasma density, ${\omega}_0$ is the laser frequency, ${\varepsilon}_0$ is the vacuum permittivity and $-e$ and ${m}_\mathrm{e}$ are the electron charge and rest mass. For the three laser beams irradiating a wire-hemisphere target (TWH) case, a copper (Cu) nanowire array of density $100{n}_\mathrm{c}$ , radius ${r}_0=0.1{\lambda}_0$ and length ${L}_0=3{\lambda}_0$ is attached to the front surface of the hemisphere. All wires are oriented toward the center of the hemisphere and are placed with angular uniformity. The angle between two adjacent wires is $\alpha =10{}^{\circ}$ . Three $\mathrm{p}$ -polarized laser beams incident at angles $\theta =-60{}^{\circ}$ , 0° and 60° relative to the $x$ direction are focused on the target front surface. The spatiotemporal profile of each laser beam is $\exp \left[-{\left(r/{\sigma}_0\right)}^2\right]$ and ${\sin}^2\left(\pi t/2{\tau}_0\right)$ , where ${\sigma}_0=3{\lambda}_0$ is the focal spot radius, ${\tau}_0=10{T}_0$ is the pulse duration and ${T}_0=3.3$ fs is the laser cycle. The normalized laser amplitude is ${a}_0={eE}_0/{m}_\mathrm{e}{\omega}_0c=5.8$ , where ${E}_0$ is the peak electric field strength of the laser and $c$ is the speed of light in vacuum. This corresponds to laser intensity of $4.56\times {10}^{19}$ W/cm², laser power of 6.4 TW and laser energy of 0.16 J. For the cases with a single laser beam (a single laser beam irradiating a wire-hemisphere target (SWH), a single laser beam irradiating a hemisphere target (SH) and a single laser beam irradiating a planar target (SP)), ${a}_0=10$ is used, so that the total laser energy in the five cases is the same.

Figure 2 (a) Evolution of the TNSA proton flux $\mathrm{\mathcal{F}}$ for the five cases in Figure 1. Here, only protons with energies higher than 0.5 MeV are counted. The blue dots and orange squares are for the peak flux ${\mathrm{\mathcal{F}}}_\mathrm{p}$ and total number ${\mathcal{N}}_\mathrm{t}$ , respectively, of the protons. (b) Angular distributions of the protons behind the targets at $t=100{T}_0$ . The blue dots and orange squares are for the averaged angular deviation $\left\langle \Delta \varphi \right\rangle$ (i.e., the average angle between the momentum direction of each proton and its direction with respect to the target center) of all protons, as well as the standard deviation of the proton angular distribution $\sigma$ . For the SP case, the averaged angular deviation, standard deviation and proton number are divided by 2, 4 and 6, respectively. (c) Proton energy spectra at $t=100{T}_0$ . The proton temperatures (obtained from the gradients of the curves) are also given. (d) The maximum laser-to-proton (with energies higher than 0.5 MeV) energy conversion efficiency ${\eta}_\mathrm{p}$ for the five cases.

3 Results and discussion

We record the properties of the protons entering an imaginary (i.e., protons can pass through but not be reflected from it, and nothing happens inside) spherical sample of radius $1{\lambda}_0$ located at the center of the target-back hollow, as indicated by the green circle in Figure 1. The evolution of the proton flux $\mathrm{\mathcal{F}}$ (shadow colored slice profiles) of the five cases is shown in Figure 2(a). One can see that in the TWH case, the peak flux ${\mathrm{\mathcal{F}}}_\mathrm{p}$ and total number (i.e., cumulative flux) ${\mathcal{N}}_\mathrm{t}$ of the protons are $1.55\times {10}^{30}$ cm⁻² s⁻¹ and $1.53\times {10}^{16}$ , respectively. They are an order of magnitude higher than that of the SP case, and about double or triple that of the SH and TH cases, since the wires lead to higher laser energy absorption efficiency^[ Reference Jiang, Zou, Zhao, Hu, Han, Yu, Yu, Yin and Shao^32] and more uniform laser intensity distribution on the hemisphere front surface. This includes the contributions of the wire array, laser beamlets and curved structure. Although the single-pulse intensity in the TWH case is only one-third of that in the SWH case, ${\mathrm{\mathcal{F}}}_\mathrm{p}$ is comparable and ${\mathcal{N}}_\mathrm{t}$ is even higher than in the SWH case. Moreover, the full-width at half-maximum (FWHM) of $\mathrm{\mathcal{F}}$ in the TWH case is twice that of the SWH case, indicating the high $\mathrm{\mathcal{F}}$ can remain for an extended time. This is because in the TWH case the spatial uniformity of the three pulses can counteract the loss of the low single-pulse intensity. Figure 2(b) shows the angular distributions of the protons behind the targets at $t=100{T}_0$ . It is found that for single laser beams the divergence angle of protons is small. With the use of multiple laser beams, the accelerated protons are much more uniformly distributed within $-90{}^{\circ}<\varphi <90{}^{\circ}$ , where $\varphi$ is the angle of the proton momentum with respect to the $x$ direction (with the anticlockwise direction as positive). For the TH case, one can see that there are three small peaks, corresponding to the three laser beams.

Figure 3 (a) Hot-electron energy-density distributions in the region behind the wire-hemisphere structure at $t=100{T}_0$ . The solid black curve shows a typical constant-energy-density contour relatively far away from its back surface. (b) Energy-density distribution of the wire electrons at $t=45{T}_0$ and $55{T}_0$ for the (b1), (b2) TWH and (b3), (b4) SWH cases, respectively. (c) Evolution of the total electron energy density in the hollow behind the targets. The red and yellow dashed curves are for the contribution of the wire electrons in the TWH and SWH cases. (d) Electron-energy spectra at $t=55{T}_0$ . The blue and black dashed curves show the spectra of the wire electrons in the TWH and SWH cases. The inset shows the electron energy spectra in the $x>20{\lambda}_0$ (green and blue curves) and $x<20{\lambda}_0$ (black and red curves) regions of panels (b2) and (b4) for the TWH and SWH cases. (e) Laser-to-electron energy conversion efficiencies. The overlapping inner bars in the TWH and SWH cases are those of the laser-to-wire electrons only, from which one can clearly see the effect of the wires.

Of interest is that the proton angular distribution in the TWH case is much more uniform than that of the TH case. This is because the most intense part of each beam is on the wires. That is, the wires act like a buffer for the Gaussian lasers, resulting in a significantly more uniform angular distribution of the light intensity. Angular uniformity of the protons can be described by two parameters, namely, the averaged angular deviation $\left\langle \Delta \varphi \right\rangle ={\sum}_{i=1}^N\left|\Delta {\varphi}_i\right|/N$ of all protons and the standard deviation of their angular distribution:

(1) $$\begin{align}\sigma =\sqrt{\sum \limits_{i=1}^q{\left({N}_i-N/q\right)}^2/q},\end{align}$$

where $\Delta {\varphi}_i$ is the angle between the momentum direction of each proton and the direction of the proton toward the sample center, $N$ is the total number of protons behind the target and $\varphi \in \left[-90{}^{\circ},90{}^{\circ}\right]$ . We divide this ${180}^{\circ }$ range into $q$ equal parts (here $q=180$ ), and the proton number of each part is ${N}_i$ , so that $N={\sum}_{i=1}^q{N}_i$ . One can see that for the SP case, $\left\langle \Delta \varphi \right\rangle$ and $\sigma$ are about 29° and $1.5\times {10}^{15}$ , respectively, whereas in the TWH case, they are only 5.5° and $6.2\times {10}^{13}$ . That is, in the TWH case a relatively large number of protons can converge toward the sample. Note that the angular uniformity of protons in the TWH case is also obviously better than that of the SWH case. Therefore, the TWH case has significant advantages in producing high-flux and angularly uniform proton beams. The proton spectra at $t=100{T}_0$ and laser-to-proton energy conversion efficiencies ${\eta}_\mathrm{p}$ are shown in Figures 2(c) and 2(d). We can see that the TWH setup is also advantageous for enhancing the proton energy. The maximum proton energy ${\varepsilon}_{\mathrm{p},\max }$ , proton temperature ${T}_\mathrm{p}$ and ${\eta}_\mathrm{p}$ are about 15 MeV, 4.9 MeV and $20.6\%$ (slightly lower than that of the SWH case with high single-pulse intensity), respectively, and for the SP case they are only about 5 MeV, 1.2 MeV and $3.4\%$ , respectively. Note that ${\varepsilon}_{\mathrm{p},\max }$ and ${\eta}_\mathrm{p}$ for the TWH and SWH cases are considerably higher than that of the TH, SH and SP cases, indicating the effectiveness of the wire microstructure. Moreover, ${\varepsilon}_{\mathrm{p},\max }$ and ${\eta}_\mathrm{p}$ for both the TH and SP cases are lower than that of the SH case. This is because in addition to the curvature of the target, proton acceleration is also affected by the laser intensity^[ Reference Daido, Nishiuchi and Pirozhkov^9, Reference Macchi, Borghesi and Passoni^10].

Next we consider the causes for the uniform distribution of the generated high-flux proton beams. Figure 3(a) compares the energy-density distribution of the hot electrons behind the five targets at $t=100{T}_0$ . One can see that in the TWH case, the electrons have much more uniform angular distribution and higher energy density. Figure 3(b) shows the energy-density distributions of electrons from only the wires in the TWH and SWH cases at $t=45{T}_0$ and $55{T}_0$ . They are in the form of small periodic beamlets. This is because the dense attosecond electron bunches are pulled out from the surface of metal wires by the transverse electric field of the laser^[ Reference Jiang, Zou, Zhao, Hu, Han, Yu, Yu, Yin and Shao^32, Reference Yu, Zou, Yu, Yu, Yin, Shao, Zhuo, Zhou and Ruan^33]. They are then accelerated forward by the ponderomotive force of the laser as well as the longitudinal electric field of a transverse magnetic (TM) mode that is excited in the vacuum channels between the wires^[ Reference Yu, Zou, Yu, Yu, Yin, Shao, Zhuo, Zhou and Ruan^33]. When these electrons pass through the hemisphere, the laser pulse is reflected back by the high-density hemisphere shell. Thus, the electrons continue to move forward in almost the original direction over a short distance of a few micrometers. As these bunches leave the shell, some are well confined near the back surface of the shell by the space-charge sheath electric field ${E}_\mathrm{s}$ , while others escape from the sheath field and move into the $x>20{\lambda}_0$ region. The wire electrons are uniformly distributed in space and are responsible for the angularly uniform sheath fields and the subsequent ion acceleration. In contrast, in the SWH, SH and SP cases, the hot electrons are concentrated mainly in the axial region behind focal spot. Figure 3(c) shows the temporal evolution of the total energy density of all electrons behind the target. We see that at all times the energy density in the TWH case is several times higher than that in the cases of TH, SH and SP. Almost all high-energy-density electrons are from the Cu wires and produced between $t=40{T}_0$ and $60{T}_0$ . As shown in Figure 3(d), the highest hot-electron temperature and maximum energy in the TWH case are about 2.6 and 17 MeV, respectively, and most electrons with energy greater than 2 MeV are from the wires. It is important to distinguish the high-energy region of wire electrons in the spectrum. The inset shows the spectra of the electrons in the $x<20{\lambda}_0$ and $x>20{\lambda}_0$ regions of Figures 3(b2) and 3(b4) for the TWH and SWH cases. In both cases, the electrons at $x>20{\lambda}_0$ have a lower temperature and cut-off energy because these electrons lose most of their energies as they escape the sheath field. The electrons trapped in the sheath electric field at $x<20{\lambda}_0$ are more energetic because they can also gain energy when reflected back to the hemisphere. Figure 3(e) shows the laser-to-electron energy conversion efficiency ${\eta}_\mathrm{e}$ for the five cases. We see that ${\eta}_\mathrm{e}$ is as high as $70\%$ in the TWH case, which is much higher than that of the other cases. The energy conversion efficiency from the laser to the wire electrons is up to 61%.

Figure 4 (a) Distributions of the electric field strength behind the five targets at $t=50{T}_0$ . The arrows show the electric field magnitudes and directions. (b) Profiles of the sheath electric field at a distance $0.5{\lambda}_0$ perpendicular to the local target-rear surface. (c) Standard deviation of the sheath electric fields shown in (b).

Figure 5 Proton energy-density distribution in the circular region behind the five targets at $t=100{T}_0$ .

Figure 4(a) is for the sheath electric field ${E}_\mathrm{s}$ behind the targets at $t=50{T}_0$ . We see that the electric fields are directed toward the center region of the hollow. In addition, ${E}_\mathrm{s}$ in the TWH case is stronger than that in the TH case due to the higher sheath-electron temperature, as discussed earlier (recall also that the sheath field of the hot electrons is given by ${E}_\mathrm{s}\propto \sqrt{n_\mathrm{e}{T}_\mathrm{e}}$ ^[ Reference Mora^34]). For completeness, Figure 4(b) shows some typical electric fields at a distance $0.5{\lambda}_0$ away from the local rear surface of the targets. One can see that the field intensity distribution in the TWH case is generally uniform, with only a few relatively small fluctuations. The uniformity of the sheath fields for TNSA of the target protons can be seen from the standard deviations of their strength distributions:

(2) $$\begin{align}\sigma \left({E}_\mathrm{s}\right)=\sqrt{\sum \limits_{i=1}^q{\left({E}_{\mathrm{s},i}-\sum \limits_{i=1}^q{E}_{\mathrm{s},i}/q\right)}^2/q},\end{align}$$

where ${E}_{\mathrm{s},i}$ is the field strength at each simulation grid. Figure 4(c) shows $\sigma \left({E}_\mathrm{s}\right)$ for the five cases. We see that the TWH and TH cases have the lowest $\sigma \left({E}_\mathrm{s}\right)$ and thus are the most uniform. Figure 5 shows the energy-density distributions of the protons in a circular region behind the target at $t=100{T}_0$ . As expected, in the TWH case there is a large spot of high-energy-density protons.

It should be noted that the simulation is 2D, which tends to overestimate the hot-electron temperature, resulting in higher proton energy and flux^[ Reference Liu, Chen, Zheng, Sheng and Liu^35, Reference Héron, Adam and Mora^36]. Moreover, the longitudinal electric field of waveguide TM mode in the vacuum regions (no longer channels in three dimensions) between the wires could differ from that of two dimensions and can thus also affect the proton acceleration^[ Reference Yu, Zou, Yu, Yu, Yin, Shao, Zhuo, Zhou and Ruan^33].

Figure 6 Dependence of the peak flux ${\mathrm{\mathcal{F}}}_\mathrm{p}$ of protons on (a) the wire radius ${r}_0$ , angle $\alpha$ between the two adjacent wires and (b) wire length ${L}_0$ in the TWH case. The values (in units of $1\times {10}^{30}$ cm⁻² s⁻¹) marked in (a) are for ${\mathrm{\mathcal{F}}}_\mathrm{p}$ .

4 Parametric effects and experimental considerations

In the TWH case, the effect of proton acceleration is closely related to the geometric parameters of the wire array. Figure 6 shows the dependence of ${\mathrm{\mathcal{F}}}_\mathrm{p}$ on the radius ${r}_0$ , length ${L}_0$ of the wires and angle $\alpha$ between the two adjacent wires. One can see that the TWH case is robust in producing high-flux proton beams over a relatively wide range of wire parameters, as shown in Figure 6(a). For the same $\alpha$ , ${\mathrm{\mathcal{F}}}_\mathrm{p}$ is higher as ${r}_0$ becomes smaller since less laser energy is reflected. However, when ${r}_0$ is fixed, there is an optimal $\alpha$ to provide enough electrons while reducing light reflection. Figure 6(b) shows the peak flux ${\mathrm{\mathcal{F}}}_\mathrm{p}$ of protons at different ${L}_0$ for $\alpha =10{}^{\circ}$ and ${r}_0=0.1{\lambda}_0$ . We find that the optimal length is ${L}_0=3{\lambda}_0$ to obtain the highest ${\mathrm{\mathcal{F}}}_\mathrm{p}$ . This is because in this case the interval between the two adjacent wires at the opening is about $2.3{\lambda}_0$ , and the distance between the tip of the wires and the back surface of the hemisphere is $6{\lambda}_0$ . The focusing condition of the electron beams is thus satisfied, so that they are well concentrated on the target-back surface, resulting in more intense sheath electric fields^[ Reference Zou, Pukhov, Yi, Zhou, Yu, Yin and Shao^37]. The optimum parameters ${r}_0=0.1{\lambda}_0$ and ${L}_0=3{\lambda}_0$ are also comparable to that in recent experiments^[ Reference Curtis, Calvi, Tinsley, Hollinger, Kaymak, Pukhov, Wang, Rockwood, Wang, Shlyaptsev and Rocca^38, Reference Kong, Zhang, Shou, Xu, Mei, Cao, Pan, Wang, Qi, Lou, Ma, Lan, Wang, Li, Rubovic, Veselsky, Bonasera, Zhao, Geng, Zhao, Fu, Luo, Ma, Yan and Ma^39].

Figure 7(a) shows the effect of the laser incident angle $\theta$ on the proton angular distribution and flux for the TWH case. The angles $\theta$ of the upper and lower wires are symmetrically varied. As expected, as $\theta$ is decreased, the proton distribution becomes more like that of the single-laser case SH. As $\theta$ increases, the proton distribution becomes more uniform. However, as shown in Figure 7(b), $\theta$ has little effect on ${\mathrm{\mathcal{F}}}_\mathrm{p}$ and ${\mathcal{N}}_\mathrm{t}$ . This is because at small incident angles laser interference can enhance the laser field strength^[ Reference Jiang, Shao, Zou, Yu, Hu, Guo, Huang, Zhang, Wu, Zhang, Yu, Yin, Zhuo and Zhou^40], which compensates for the deviation in proton uniformity. Due to experimental factors such as laser spot jitter, the actual incident direction of the laser would have angle deviation $\Delta \theta$ and transverse drift $\Delta s$ from the expected direction. Figures 7(c)–7(f) show the influence of $\Delta \theta$ and $\Delta s$ on the angular distribution and flux of the TNSA protons, respectively. We can see that with increasing $\Delta \theta$ and $\Delta s$ , the angular uniformity of the proton beam and ${\mathrm{\mathcal{F}}}_\mathrm{p}$ decrease slightly, but ${\mathcal{N}}_\mathrm{t}$ remains roughly the same. This indicates that the TWH case is relatively robust in obtaining high-flux angularly uniform proton beams. It is worth noting that the above conclusions are on a small focal spot with 3 μm radius. For a large spot, the influence of $\theta$ , $\Delta \theta$ and $\Delta s$ on the angular distribution and flux of protons is relatively weakened.

For higher power laser pulses, it is necessary to consider the damage of the microstructure by the laser prepulses. To see the effect of the pre-expansion preplasma produced by the prepulse, we performed radiation hydrodynamic simulations using FLASH code^[ Reference Fryxell, Olson, Ricker, Timmes, Zingale, Lamb, Macneice, Rosner, Truran and Tufo^41]. The laser prepulse consists of an amplified spontaneous emission (ASE) pedestal and a picosecond exponentially rising leading edge. The pulse shape and intensity of the prepulse are calculated by the laser peak intensity and contrast. The simulations are performed in 2D cylindrical geometry, using the equation-of-state and multi-group opacity tables from PROPACEOS68^[ Reference MacFarlane, Golovkin, Woodruff and Quant^42]. FLASH uses an adaptive mesh refinement (AMR) scheme and the coarsest/finest mesh size used is 0.1/0.02 μm in both directions. The Courant–Friedrich–Lewy (CFL) time limit is 0.4. The laser deposition is by 3D ray-tracing projection on the cylindrical domain.

Figure 7 (a), (c), (e) Angular distribution and (b), (d), (f) peak flux ${\mathrm{\mathcal{F}}}_\mathrm{p}$ (left) and total number ${\mathcal{N}}_\mathrm{t}$ (right) of the TNSA protons versus the laser incident angle $\theta$ , misalignment $\Delta \theta$ of the laser incident angle and transverse drift $\Delta s$ of the focal spot. The panels on the far left are for $\theta =30{}^{\circ}$ , with misalignment $\Delta \theta =5{}^{\circ}$ and transverse drift $\Delta s=3{\lambda}_0$ .

Figure 8 Results of radiation hydrodynamic simulations. (a) Density distribution of preplasma produced by the prepulse at different laser power $P$ , and its (b) angular and (c) radial profiles along the white lines in (a). The density is in logarithmic color scale and these data are extracted 3 ps before the peak of the main pulse arrived. The black curves in (a) are the contour of the overcritical region. The above density distributions from FLASH simulations are then used in the PIC code EPOCH as input parameters of the wire-hemisphere target. Dependence of (d) the peak flux ${\mathrm{\mathcal{F}}}_\mathrm{p}$ and (e) angular distribution of protons on $P$ both without and with preplasma from the simulation.

We now consider the effect of the laser power $P$ (from 6.4 TW to 1 PW), with contrast ${10}^{10}$ . Figure 8(a) shows the density distributions of preplasma produced by the prepulse under five different values of $P$ in the TWH case. We see that as $P$ increases, the surface of the wire-hemisphere target is ablated to smooth the sharp boundaries, and more preplasma of overcritical density fills the gaps between the wires. When $P\le 200$ TW, the density of the preplasma is usually below the critical plasma density. As $P$ continues to rise, the preplasma in most regions is over critical density, as shown in Figures 8(b) and 8(c). The density distributions from FLASH simulations are then imported into EPOCH code as the initial density distributions of the wire-hemisphere target. Figure 8(d) shows the dependence of ${\mathrm{\mathcal{F}}}_\mathrm{p}$ on $P$ without and with preplasma. We find that the preplasma has little effect on ${\mathrm{\mathcal{F}}}_\mathrm{p}$ as $P$ is less than 200 TW, since the preplasma has low enough density to be classically transparent. However, when $P$ is higher, the overcritical preplasma between the wires will block the channels and prevent laser injection, resulting in a dramatic reduction of ${\mathrm{\mathcal{F}}}_\mathrm{p}$ by approximately 50%. This indicates that for high-contrast lasers with low energy and modest intensity, the effect of preplasma on proton acceleration in the proposed scheme is negligible. For long pulses of high-energy lasers, the contrast needs to be further improved to avoid ablation of the microstructure. The effect of preplasma at different values of $P$ on the angular distribution of protons is shown in Figure 8(e). We see that with increasing $P$ , the angular distribution of the protons is more uniform, since the roughness of the periodic microstructure is diminished by the ablative effect of the stronger laser.

5 Summary

A practical scheme of producing high-flux and angularly uniform proton acceleration using three laser beams interacting with a wire-hemisphere target is proposed and demonstrated by PIC simulations. The peak proton flux and angular uniformity are significantly better than that from the traditional methods. The scheme can be useful for creating warm dense plasma states that are relevant to laboratory astrophysics^[ Reference Remington, Arnett, Drake and Takabe^43], as well as fundamental material property^[ Reference Collins, Da Silva, Celliers, Cold, Foord, Wallace, Ng, Weber, Budil and Cauble^14] and controlled fusion^[ Reference Nuckolls, Wood, Thiessen and Zimmerman^44] research. It should, however, be mentioned that we have invoked ultrashort relativistic laser pulses with duration of approximately 33 fs and energy of approximately 0.16 J per beam. For applications such as proton-beam-driven fast ignition of ICF, the protons must deposit sufficient energy (10–20 kJ) in a hot spot to ignite the fuel^[ Reference Atzeni, Schiavi, Honrubia, Ribeyre, Schurtz, Nicolaï, Olazabal-Loumé, Bellei, Evans and Davies^45]. That is, long-pulse (a few ps) high-energy (hundreds of kJ) lasers are required. With such high-energy picosecond laser pulses, even with high contrast, microstructures such as that used here would be rapidly filled with overcritical plasma^[ Reference Bailly-Grandvaux, Kawahito, McGuffey, Strehlow, Edghill, Wei, Alexander, Haid, Brabetz, Bagnoud, Hollinger, Capeluto, Rocca and Beg^24] that could distort the proton acceleration.