Spin-polarized 3 He shock waves from a solid-gas composite target at high laser intensities

We investigate collisionless shock acceleration of spin-polarized 3 He for laser pulses with normalized vector potentials in the range a 0 = 100 − 200. The setup utilized in the 2D-particle-in-cell simulations consists of a solid carbon foil that is placed in front of the main Helium target. The foil is heated by the laser pulse and shields the Helium from the highly oscillating fields. In turn, a shock wave with more homogeneous fields is induced, leading to highly polarized ion beams. We observe that the inclusion of radiation reaction into our simulations leads to a higher beam charge without affecting the polarization degree to a significant extent.


Introduction
In recent years, the acceleration of spin-polarized particle beams from laser-plasma interaction has gained a lot of interest.Applications i.a.include probing the nuclear structure of the proton using deep-inelastic scattering [1], or polarized fusion, as the cross section for fusion is increased for spin-polarized reactants [2].
Most of the recent publications utilize pre-polarized plasma targets, either consisting of HCl [3,4] or 3 He, which will be considered here.Such a Helium target has recently been experimentally realized and used in a first experiment at PHELIX (GSI Darmstadt) [5].The polarization process is described in detail in the work by Fedorets et al. [6].Pre-polarization is generally required, as no net polarization of the ions can be gained during laser-plasma interaction [7].For electron acceleration, an in-situ method of obtaining polarized witness beams in wakefields has been proposed by Nie et al. [8,9].A general overview of the state-of-the-art for laser-plasma based acceleration of spin-polarized particle beams is given by Büscher et al. [10].
The crucial problem with using pre-polarized targets is that the polarization has to be maintained during the acceleration process.However, the strong electromagnetic fields required for high-energy particle beams induce spin precession according to the T-BMT equation [11,12].This was also observed in particle-in-cell (PIC) simulations that study the acceleration of spin-polarized ion beams via Magnetic Vortex Acceleration [13].The oscillating laser fields as well as the field structure prevalent in the created plasma channel lead to the depolarization of the accelerated ion beam.Jin et al. showed that for a laser with a 0 = 25, the final beam polarization is in the range of 82%, while for a 0 = 100, only 56% could be obtained.Here, a 0 = eE 0 /(m e cω 0 ) denotes the normalized laser vector potential and e is the elementary charge, E 0 the peak laser electric field, m e the electron rest mass, c the vacuum speed of light and ω 0 the laser frequency.With the aim of utilizing polarized particle beams in combination with laser intensities in the range of 10 23 W/cm 2 [14], or near-future facilities even exceeding that [15,16], different acceleration schemes will be necessary.
In a proceeding publication [17], a dual-pulse Magnetic Vortex Acceleration setup was proposed that uses two co-propagating laser pulses with a carrier envelope phase difference of π.The presence of the two laser pulses creates an accelerating region for ions with decreased depolarization.The setup delivers proton beams with a polarization of 77% even when two pulses with a 0 = 100 are utilized.Still, this method will yield decreased polarization for higher intensities as the accelerated particles will still be subject to the strong laser fields.
An alternative mechanism avoiding this problem is that of Collisionless Shock Acceleration [18,19,20].One option to realize this mechanism is to place a solid foil in front of a gaseous target.When the laser pulse irradiates the foil, it heats up the electrons and displaces them with respect to the ions.In turn, an accelerating electric field is induced by which the ions in the gas are being reflected.Depending on the laser and target parameters, the laser pulse will not penetrate the foil.Thus, the oscillating laser fields do not reach the polarized ions in the gaseous target, preventing strong depolarization.This scheme of Collisionless Shock Acceleration (CSA) has been proposed for spin-polarized proton beams from an HCl target by Yan et al. [21,22].In a separate publication, they investigated shock acceleration using a micro-structured foil [23].Alternatively, CSA can also be realized without a foil, utilizing just a density-ramp [24].
In this paper, we extend the investigation of CSA to regimes of a 0 = 100 − 200.In our simulations, we utilize a 3 He target, whose density is based on current experimental capabilities [6].We include the effect of radiation damping which becomes important for the electron motion.In turn, the change in electron motion leads to differing electromagnetic fields which affect the ion motion.The setup of our PIC simulations as well as the included spin and radiation effects are explained in section 2. The results of the parameter scan concerning laser intensity, foil thickness and the inclusion of radiation reaction are presented in section 3. The relevance of other effects concerning spin polarization besides precession according to T-BMT is discussed in section 4.

PIC simulation setup
In the following, we present the simulation setup used for the parameter scans.For all of the 2D simulations, we have used the particle-in-cell code vlpl [25,26].The simulation box has a size of 120λ L × 30λ L with a grid size of h x = 0.012λ L , h y = 0.03λ L (with laser wavelength λ L = 800 nm).From here on out, x is the direction of laser propagation.As we utilize the rhombi-in-plane Maxwell solver [27], our time step is required to be ∆t = h x /c.For all of the simulations, the linearly polarized laser pulse has a duration of 10λ L /c ≈ 27 fs and a focal spot size of 5λ L = 4 µm (FWHM).At t = 0, the center of the laser pulse is located at x = −20λ L .Its normalized laser vector potential a 0 will be varied in the range 100-200.
The target the laser pulse interacts with consists of two components, a Carbon foil and a spin-polarized 3 He gas (cf.Fig. 1 for a schematic).The thickness of the Carbon foil will be varied throughout the parameter scans, but is in the range of multiple laser wavelengths.It is modelled as a rectangular slab with a Carbon ion density of 33.3n cr .For all simulations, the foil is simulated with 8 particles per cell and -for reasons of computational efficiency -fully ionized.The polarized Helium is modelled as a homogeneous gaseous target with 2λ L of linear up-ramp that is placed as pre-plasma in front of the foil itself.Afterwards, the Helium target has its peak density of 0.02n cr for 97λ L followed by a 20λ L down-ramp back to vacuum.Initially, the Helium is fully ionized, and spin-polarized in y-direction, i.e. s y = 1 for all ions.Two particles per cell are used for the Helium ions.
The spin precession of the macro-particles is calculated according to the T-BMT where Ω is the precession frequency given as The pre-factors in front of the terms containing electric field E, magnetic field B and particle velocity v are , Here, a is the particle's anomalous magnetic moment and γ its Lorentz factor.Since the electric field of the shock is pointed in x-direction, i.e. the direction of ion acceleration, the term v × E vanishes.Further, magnetic fields are mostly negligible for shock acceleration, meaning that the precession frequency Ω will be small.The relevance of other spin-related effects is discussed in section 4 as well as the publication by Thomas et al. [28].
For the investigation of the effects of radiation on the spin-polarization of the beam we use the following description: the radiation reaction (RR) force is calculated as where denotes the quantum parameter.The so-called Gaunt factor incorporates the fact that charges will emit less if χ → 1, i.e. when QED effects become important.As this integral is computationally expensive, we use the approximation which has been shown to be a reasonable simplification [29].For the macro-particles in the PIC simulations, the radiation-corrected momentum is then calculated as This effect will be switched on and off in the following simulations to observe the consequences for the final Helium beam.

Results
We will first discuss the results obtained without considering radiation reaction.When the laser pulse irradiates the foil target, it heats up the electrons.The displacement of the electrons with respect to the ions of the foil creates an electric field which reflects the Helium ions.The reflected ions propagate in the form of a shock wave (cf.Fig. 2 at t = 213 fs).If the foil remains opaque for the laser pulse, the Helium ions only experience the electric field of the shock.This field is much more homogeneous than the oscillating laser fields, thus maintaining a higher spin polarization than in the case of Magnetic Vortex Acceleration, where the laser pulse directly interacts with the target.Going to high laser intensities or very thin foils, the foil will be penetrated by the laser pulse.This, in turn, will expose the Helium to parts of the laser pulse itself.This will add further acceleration to ions, but also lower polarization as they now experience the oscillating laser fields.For a 0 = 200, we observe signs of a MVA-like plasma channel after t = 373 fs (cf.Fig. 2).

Laser intensity
The influence of the laser intensity on final Helium ion energies is in accordance with previous studies on the CSA process: higher intensity induces a stronger shock potential, i.e. higher energies are achieved.In the case of a foil with 2λ L thickness, a well-defined energy peak can be observed in all simulations for the range a 0 = 100 − 200 (cf.Fig. 3).For a 0 = 100, the peak energy is around E p ≈ 370 MeV, while for a 0 = 200 up to 943 MeV are obtained (see Table 1).The scaling found from our simulations is E p ∝ a 1.34  0 .This scaling is within the theoretical prediction of [18] and stronger than the typical scaling of Target Normal Sheath Acceleration.The relative width of the peak (FWHM) drops from 11.5% to 6.6% for higher intensities meaning that the condition for particle trapping into the shock wave becomes more restrictive for higher a 0 .The angular spectrum of the accelerated beam will depend on the laser parameters (duration, focal spot size) as well as the foil thickness.
The amount of charge accelerated within this peak is on the level of tens of pC.Higher intensities lead to lower beam charge as the shock is less able to trap a significant amount of particles the faster it becomes.For a 0 = 200, it drops to about 11.3 pC.This further indicates that our setup is limited by the low Helium density.While choosing a higher density in simulations should improve the amount of accelerated ions, experimentally the Helium-3 target is currently limited to densities of a few 10 19 cm −3 [6].Thus, beam charge can only be further controlled by our choice of laser and foil parameters.
The degree of polarization is calculated as P = P 2 x + P 2 y + P 2 z , where P j = i s i,j /N is the average over the spin components of all N particles in one direction j ∈ {x, y, z}.While the minimum degree of polarization per energy bin generally decreases with laser intensity (cf.Fig. 4), it stays on a 90% level throughout the a 0 scan.Considering only the ions in the FWHM around the energy peak, the degree of polarization is even higher (although the direction of the particle spins may have rotated uniformly compared to their initial direction).Once the laser pulse fully penetrates foil, the polarization abruptly drops due to the highly oscillating laser fields.The degree of polarization with the CSA setup is generally much higher than what was obtained in the MVA studies [13,17].This is in accordance with the scaling laws derived by Thomas et al. [28].   .Energy spectra for different foil thicknesses (a 0 ≡ 200, RR off).The maximum energy decreases for thicker foils due to the increased shielding of the laser fields and the weaker induced electric shock field.In the case of very thin foils (≤ 1λ L ), the laser pulse can easily penetrate the foil leading to significantly higher energies, but also highly reduced polarization of the ions.

Foil thickness
Keeping the laser intensity fixed at a 0 = 200, we also perform a scan of the foil thickness.For very thin foils ≤ 1λ L , the laser easily penetrates the foil, leading to significantly higher energies but also low beam polarization.In this case, CSA will no longer be the prominent acceleration mechanism at later stages, but rather processes like MVA (as indicated by the similar plasma channel at t = 373 fs, cf.Fig. 2) or, potentially, others.
For thicker foils in the range (2 − 5)λ L , we see the expected trend of decreasing ion energy: with increasing thickness of the foil, the rear part is heated less sufficiently by the laser pulse, which reduces the shock potential accordingly.At a 0 = 200, the original E p ≈ 1 GeV for 2λ L drops to 846 MeV for a foil with 5λ L thickness (cf.Fig. 5).In all cases, the peak feature is preserved with the relative energy spread staying at a 10% level.The beam charge is increased from the initial 11.3 pC at 2λ L to 28.7 pC which is due to the modified shock dynamics.Polarization within the FWHM around E p does not change significantly with foil thickness, while the minimum polarization in the whole simulation domain slightly increases due to improved shielding from the laser fields.

Radiation reaction
The simulations of the previous sub-sections have not included radiation reaction.In the regime of laser-plasma interaction considered here, however, electrons will radiate part of the energy, thus affecting their motion over the course of the interaction.In particular, the inclusion of RR leads to a slower propagation of the shock wave and therefore lower energy of the accelerated ions.The difference in shock velocity can be seen in Figure 6.
For a 0 = 100, the peak energy drops to approximately 353 MeV, while for a 0 = 200 it reaches 674 MeV (see Table 1 for an extensive list of the results).Here it has to be noted, that in the latter case, the shape of the energy spectrum changes significantly: while the main energy peak is moved to lower energies when considering RR, some high-energy ions in the region of the non-RR energy peak remain (see the second peak around 950 MeV in Fig. 7).These ions still experience a fast-moving accelerating field.By comparison, at lower intensities, the shape of the energy spectrum remains the same and only is dampened by RR, i.e. acceleration of the ions is performed solely by the shock wave.The relative width of the main peak stays rather constant at around 10% for all simulations.
A benefit of radiating damping and the slower moving shock wave is that more particles can be trapped by it: regardless of laser intensity, the beam charge in the FWHM around the energy peak stays consistently at around 25 pC.This leads to the conclusion that a point of saturation with respect to beam charge is reached here, and any further increase would necessitate a change of target parameters (which is, again, experimentally constrained).Polarization of the ion beams in the FWHM remains mostly unchanged when considering RR: while the reduced fields due to RR should generally improve the degree of polarization, it stands to reason that this effect is equilibrated by trapping a larger amount of particles in the shock wave.The minimum polarization throughout the whole simulation box stays at a high degree of approximately (85 − 90)% as well, the only exception being the simulation for a 0 = 200.Here, the minimum polarization drops to about 70% in the energy region between the two peaks at 674 MeV and 950 MeV, where significantly fewer particles are located.In general, however, we obtain a consistently high degree of polarization with the benefit of accelerating more charges when considering radiation reaction.

Discussion
Throughout the conducted simulations, the only spin-related effect considered has been precession according to the T-BMT equation.Other effects that would -depending on the parameter regime -need to be considered are those of the Stern-Gerlach force [30] and radiative polarization as in the Sokolov-Ternov effect [31].Taking the scaling laws derived in the publication by Thomas et al. [28], we can conclude that the Stern-Gerlach force on the Helium ions can be neglected as the differences in particle trajectories would Comparison of the Helium density after 213 fs for simulations where radiation reaction is turned off or on, respectively.In both cases, the laser had a normalized laser vector potential of a 0 = 200 and the foil has a thickness of 2λ L .In the case with RR, the shock propagates more slowly (cmp.the red/blue shock waves in the range x = (10−20) µm as well as the bow-like structures in the region x = (20−30) µm).The colorbar is clipped at ±0.1n cr for better visibility.
be miniscule.For electrons, it is well known from several publications that in the highintensity regime, their polarization will significantly be affected by radiative processes [32,33].Technically, the spin-dependence of the Gaunt factor would also need to be considered in that case [34].In the current publication, however, electron polarization is not of interest and their trajectories will not be significantly altered.The Helium ions themselves will not radiate to a large extent, as their Lorentz factor γ = E/E 0 + 1 is comparatively low even for a 0 = 200.Thus, polarization build-up due to radiation emission can be neglected for the Helium ions.
Another effect that can be neglected here is that of electron-positron pair creation as the laser vector potential never exceeds a 0 = 200 in our simulations.It should, however, be noted that in the regime where pair production can occur, the pre-plasma in front of the foil will likely play an important role for the acceleration process: Wang et al. showed that for ultra-intense laser pulses QED cascades in the pre-plasma can lead to the formation of an opaque particle layer in front of the foil [35].This, in turn, changes the interaction of the laser pulse with the foil and the subsequent particle acceleration.
The simulations conducted here are currently restricted to a 2D geometry for reasons of computational efficiency.As it is known from several publications, this means that the beam energy of particles will be overestimated compared to 3D-PIC simulations  [36,37].As the results already have indicated, which acceleration mechanism is prevalent during laser-solid ineraction (and in some circumstances, subsequent laser-Helium interaction) will strongly depend on choice of laser and target parameters.Thus, a future, separate study will further investigate the different acceleration mechanisms during high-intensity laser-solid interaction and their consequences for beam polarization.This will include a study of the influence of laser polarization on the acceleration mechanisms.It is e.g.known from the works by Tamburini et al. that for Radiation Pressure Acceleration, circular polarized laser pulses exhibit less spatial anisotropies and no significant RR effects compared to linear polarization [38].

Conclusion
We have studied the acceleration of spin-polarized Helium-3 via Collisionless Shock Acceleration for high laser intensities.A solid Carbon foil is placed in front of nearcritical, pre-polarized Helium to shield it from the oscillating laser fields.The obtained particle beams are of significantly higher polarization than what was obtained for previous Magnetic Vortex Acceleration studies.Depending on the laser and target parameters, the foil will be penetrated laser pulse, leading to additional acceleration mechanisms and reduced beam polarization.When radiation reaction is considered, the shock wave slows down due to the reduced shock potential, thus allowing a larger amount of Helium ions to be accelerated by it while maintaining a high degree of spin polarization.Future research could investigate the possibility of utilizing density ramps to induce similar shock waves in polarized gaseous targets without the need for the Carbon foil.Going to even higher intensities for future laser facilities, the inclusion of pair production and its effects on pre-plasma opaqueness would be of interest for CSA.

Figure 1 .
Figure 1.Schematic of the simulation setup.The laser pulse irradiates a target consisting of a solid-density Carbon foil and a gaseous, pre-polarized 3 He target.

Figure 2 .
Figure 2. Helium density at two different time steps, t = 213 fs and t = 373 fs, for an a 0 = 200 laser pulse and a 2λ L thick foil.At earlier times, the shock wave can clearly be identified (high-density region around x = 20 µm).Later on, a plasma channel similar to the one in the MVA process can be observed, indicating that multiple acceleration mechanisms occur depending on the specific laser and target parameters.The colorbar is clipped at 0.2n cr for better visibility.

Figure 3 .
Figure 3. Energy spectrum in dependence of the laser intensity without radiation reaction.For all of the simulations, a foil of 2λ L is used.Energy increases with higher laser intensity.A well-defined high-energy peak can is observed for all simulations.

Figure 4 .
Figure 4. Polarization per energy bin for different laser intensities using a 2λ L foil without RR.While higher intensity laser pulses will generally lead to lower polarization, the values obtained greatly exceed those observed in MVA simulations.

Figure 5
Figure 5. Energy spectra for different foil thicknesses (a 0 ≡ 200, RR off).The maximum energy decreases for thicker foils due to the increased shielding of the laser fields and the weaker induced electric shock field.In the case of very thin foils (≤ 1λ L ), the laser pulse can easily penetrate the foil leading to significantly higher energies, but also highly reduced polarization of the ions.
Figure 6.Comparison of the Helium density after 213 fs for simulations where radiation reaction is turned off or on, respectively.In both cases, the laser had a normalized laser vector potential of a 0 = 200 and the foil has a thickness of 2λ L .In the case with RR, the shock propagates more slowly (cmp.the red/blue shock waves in the range x = (10−20) µm as well as the bow-like structures in the region x = (20−30) µm).The colorbar is clipped at ±0.1n cr for better visibility.

Figure 7 .
Figure 7. Energy spectra obtained from simulations with a 2λ L foil at different laser intensities under consideration of radiation reaction.While the peak energy generally is reduced compared to the simulations without RR, in the case of a 0 = 200 highenergy ions exceeding the peak can be observed, indicating that other acceleration mechanisms besides CSA occur in this regime.

Table 1 .
Peak energy, peak width, charge and minimum polarization from simulations with different a 0 without and with radiation reaction.For all simulations, the foil was 2λ L thick foil.Note that the peak width and the charge correspond to the FWHM around the peak.