Hostname: page-component-cd9895bd7-q99xh Total loading time: 0 Render date: 2024-12-27T05:11:53.624Z Has data issue: false hasContentIssue false
  • English
  • Français

Electron acceleration at grazing incidence of a subpicosecond intense laser pulse onto a plane solid target

Published online by Cambridge University Press:  28 December 2015

N. E. Andreev ,
L. P. Pugachev ,
M. E. Povarnitsyn  and
P. R. Levashov
N. E. Andreev
Affiliation:
Joint Institute for High Temperatures RAS, 13 Izhorskaya str., Bldg. 2, Moscow 125412, Russia
L. P. Pugachev
Affiliation:
Joint Institute for High Temperatures RAS, 13 Izhorskaya str., Bldg. 2, Moscow 125412, Russia
M. E. Povarnitsyn*
Affiliation:
Joint Institute for High Temperatures RAS, 13 Izhorskaya str., Bldg. 2, Moscow 125412, Russia
P. R. Levashov
Affiliation:
Joint Institute for High Temperatures RAS, 13 Izhorskaya str., Bldg. 2, Moscow 125412, Russia
*
Address correspondence and reprint request to: M. E. Povarnitsyn, Joint Institute for High Temperatures RAS, 13 Izhorskaya str., Bldg. 2, Moscow 125412, Russia. E-mail: povar@ihed.ras.ru
Rights & Permissions [Opens in a new window]

Abstract

Generation of hot electrons at grazing incidence of a subpicosecond relativistic-intense laser pulse onto a plane solid target is analyzed for the parameters of petawatt class laser systems. We study preplasma formation on the surface of solid aluminum targets produced by laser prepulses with a different time structure. For modeling of the preplasma dynamics, we use a wide-range two-temperature hydrodynamic model. As a result of simulations, the preplasma expansion under the action of the laser prepulse and the plasma density profiles for different contrast ratios of the nanosecond pedestal are found. These density profiles are used as the initial density distributions in three-dimensional particle-in-cell simulations of electron acceleration by the main P-polarized laser pulse. Results of modeling demonstrate a substantial increase of the characteristic energy and number of accelerated electrons for the grazing incidence of a subpicosecond intense laser pulse in comparison with the ponderomotive scaling of laser–target interaction.

Keywords

Electron accelerationGrazing laser pulse incidenceLaser–matter interactionParticle-in-cell method
Type
Research Article
Information
Laser and Particle Beams , Volume 34 , Issue 1 , March 2016 , pp. 115 - 122
Copyright
Copyright © Cambridge University Press 2015 

1. INTRODUCTION

For relativistic intensities of a main subpicosecond pulse, even high-contrast beams can produce plasma on a target surface due to a long nanosecond and picosecond prepulses impact, which results in target surface modifications before the main pulse action, and influences the subpicosecond pulse interaction with the target extremely (Povarnitsyn et al., Reference Povarnitsyn, Andreev, Levashov, Khishchenko and Rosmej2012b, Reference Povarnitsyn, Andreev, Levashov, Khishchenko, Kim, Novikov and Rosmej2013; Esirkepov et al., Reference Esirkepov, Koga, Sunahara, Morita, Nishikino, Kageyama, Nagatomo, Nishihara, Sagisaka, Kotaki, Nakamura, Fukuda, Okada, Pirozhkov, Yogo, Nishiuchi, Kiriyama, Kondo, Kando and Bulanov2014). In this paper, the generation of hot electrons at grazing incidence of a subpicosecond intense laser pulse onto a plane solid target is analyzed for the parameters of the petawatt high-energy laser for heavy ion experiments (PHELIX) (Bagnoud et al., Reference Bagnoud, Aurand, Blazevic, Borneis, Bruske, Ecker, Eisenbarth, Fils, Frank, Gaul, Goette, Haefner, Hahn, Harres, Heuck, Hochhaus, Hoffmann, Javorková, Kluge, Kuehl, Kunzer, Kreutz, Merz-Mantwill, Neumayer, Onkels, Reemts, Rosmej, Roth, Stoehlker, Tauschwitz, Zielbauer, Zimmer and Witte2010; Wagner et al., Reference Wagner, Bedacht, Ortner, Roth, Tauschwitz, Zielbauer and Bagnoud2014a) using three-dimensional (3D) particle-in-cell (PIC) modeling (Pukhov, Reference Pukhov1999) and a wide-range hydro modeling (Povarnitsyn et al., Reference Povarnitsyn, Andreev, Apfelbaum, Itina, Khishchenko, Kostenko, Levashov and Veysman2012a) of the preplasma expansion under the action of the laser prepulse. Elaboration of wide-range models of the laser–matter interaction is necessary for planning and interpretation of experiments carried out in view of different applications aimed for the development of secondary sources of high energy particles and hard radiation (Morace et al., Reference Morace, Fedeli, Batani, Baton, Beg, Hulin, Jarrott, Margarit, Nakai, Nakatsutsumi, Nicolai, Piovella, Wei, Vaisseau, Volpe and Santos2014; Brabetz et al., Reference Brabetz, Busold, Cowan, Deppert, Jahn, Kester, Roth, Schumacher and Bagnoud2015; Rusby et al., Reference Rusby, Wilson, Gray, Dance, Butler, MacLellan, Scott, Bagnoud, Zielbauer, McKenna and Neely2015; Vaisseau et al., Reference Vaisseau, Debayle, Honrubia, Hulin, Morace, Nicola, Sawada, Vauzour, Batani, Beg, Davies, Fedosejevs, Gray, Kemp, Kerr, Li, Link, McKenna, McLean, Mo, Patel, Park, Peebles, Rhee, Sorokovikova, Tikhonchuk, Volpe, Wei and Santos2015). Surface acceleration and transport of energetic electrons in intense laser–matter interactions at the grazing incidence were investigated experimentally for femtosecond laser pulses of hundreds mJ energy in (Wang et al., Reference Wang, Chen, Mao, Huang, Ma, Zhao, Zhang, Yan, Li, Ma, Ly, Lu, Wei, Sheng and Zhang2013; Mao et al., Reference Mao, Chen, Huang, Ma, Zhao, Li, Yan, Ma, Aeschlimann, Wei and Zhang2015) and for hundred joule subpicosecond laser pulses of the laser system PHELIX in paper (Gray et al., Reference Gray, Yuan, Carroll, Brenner, Coury, Quinn, Tresca, Zielbauer, Aurand, Bagnoud, Fils, KГјhl, Lin, Li, Li, Roth, Neely and McKenna2011).

For different laser systems the contrast of a nanosecond prepulse of amplified spontaneous emission (ASE) can vary within wide limits from 10−6 to 10−12 and higher, when plasma mirror focusing is used. The normal ASE contrast (or pedestal) is between 10−6 and 10−7, if no special measures are taken. A typical time structure of the petawatt class laser pulse is shown in Figure 4 of (Wagner et al., Reference Wagner, João, Fils, Gottschall, Hein, Körner, Limpert, Roth, Stöhlker and Bagnoud2014b), where the pulse shape of the PHELIX laser at the first harmonic, λ0 = 1.053 µm, is shown for different parameters of the contrast boosting module. For simulation of the laser–matter interaction at non-relativistic prepulse intensities, we have elaborated and used the two-temperature single-fluid radiation hydrodynamic model (Povarnitsyn et al., Reference Povarnitsyn, Andreev, Apfelbaum, Itina, Khishchenko, Kostenko, Levashov and Veysman2012a, Reference Povarnitsyn, Andreev, Levashov, Khishchenko, Kim, Novikov and Rosmej2013). This wide-range model describes laser energy absorption, electron–ion coupling and two-temperature effects, radiation transport, thermodynamic properties of materials, and ionization from normal conditions at room temperature to weakly non-ideal high-temperature plasma. Using the model we have studied pump-probe experiments (Povarnitsyn et al., Reference Povarnitsyn, Andreev, Apfelbaum, Itina, Khishchenko, Kostenko, Levashov and Veysman2012a) as well as the action of a nanosecond prepulse on thin films (Povarnitsyn et al., Reference Povarnitsyn, Andreev, Levashov, Khishchenko and Rosmej2012b, Reference Povarnitsyn, Andreev, Levashov, Khishchenko, Kim, Novikov and Rosmej2013). As opposed to the Radiative Arbitrary Lagrangian–Eulerian Fluid dynamics in two dimensions code (Ortner et al., Reference Ortner, Faik, Schumacher, Basko, Blazevic, Busold, Bedacht, Cayzac, Frank, Kraus, Rienecker, Schaumann, Tauschwitz, Wagner and Roth2015), in our model, the laser absorption is described using the wave equation for the laser pulse electric field.

This paper is organized as follows. Preplasma formation under the action of different contrasts prepulses is modeled and described in Section 2. The results of 3D PIC modeling of the electron acceleration by the main subpicosecond relativistic-intense laser pulse (with parameters typical for the PHELIX laser system at the first harmonic, λ0 = 1.053 µm) in the plasma produced by the prepulse are presented in Section 3. Discussion and conclusion are summarized in Section 4.

2. PREPLASMA FORMATION UNDER THE ACTION OF DIFFERENT PREPULSE CONTRASTS

To describe plasma corona formation under the action of different intensity prepulses, the 1D version (Povarnitsyn et al., Reference Povarnitsyn, Andreev, Levashov, Khishchenko, Kim, Novikov and Rosmej2013) of the HELIO2D code (Andreev et al., Reference Andreev, Povarnitsyn, Veysman, Faenov, Levashov, Khishchenko, Pikuz, Magunov, Rosmej, Blazevic, Pelka, Schaumann, Schollmeier and Roth2015) was used, which takes into account 1D hydrodynamic motion of matter, laser energy absorption, two-temperature non-equilibrium states for electron and ion subsystems, electron thermal conductivity, and radiation transport in diffusion approximation. The evolution of material parameters is described using the conservation of mass, momentum, and energy of electron and ion subsystems in a single-fluid two-temperature Lagrangian form.

2.1. Picosecond laser prepulse

While a nanosecond ASE prepulse can be substantially reduced by different measures, the prepulse emission on the timescale of a few tens of picoseconds is practically always present [see Fig. 4 of (Wagner et al., Reference Wagner, João, Fils, Gottschall, Hein, Körner, Limpert, Roth, Stöhlker and Bagnoud2014b)]. The importance of the multi-picosecond target dynamics in relativistic laser–matter interactions was recently emphasized in a publication (Schollmeier et al., Reference Schollmeier, Sefkow, Geissel, Arefiev, Flippo, Gaillard, Johnson, Kimmel, Offermann, Rambo, Schwarz and Shimada2015), where it was shown that inclusion of realistic temporal laser pulse profiles in modeling is necessary for the correct description of the laser energy conversion to hot electrons.

In our modeling of the plasma formation, we analyzed firstly the bulk Al target heating by the PHELIX picosecond prepulse only during the time interval (−100, −5) ps, as shown in Figure 1 for high contrast level ~10−10 at time less than −100 ps (zero time moment corresponds to the maximum of the laser pulse intensity).

Fig. 1. Time structure of the PHELIX laser pulse at the first harmonic, λ0 = 1.053 µm, in a subnanosecond domain for high contrast level ~10−10 at time less than −100 ps. Solid (black) curve presents the cross-correlation measurement, dashed (red) curve is an approximation used in simulations. Laser intensity is normalized to the maximum of the main pulse intensity.

The electron density distribution at t = −5 ps is shown in Figure 2 for the P-polarized laser pulse with the angle of incidence 80° and maximum intensity 4 × 1019 W/cm2. For electron densities less than, or about the density of the reflection point y r, the density profile can be well approximated by the rarefaction wave exponent with a characteristic scale length L r = 1.8 µm:

(1)$${n_{\rm e}}(y) = {n_{{\rm e},{\rm r}}}\exp [(y - {y_{\rm r}})/{L_{\rm r}}],$$

${\rm where \,{\it n}_{{\rm e,r}}} \equiv {n_{\rm e}}({y_{\rm r}}) = {n_{{\rm cr}}}\mathop {\cos} \nolimits^2 {\rm \theta} \!\approx\! 0.03 \times {n_{{\rm cr}}} \!\approx\! 3 \times {10^{19}}{\rm c}{{\rm m}^{ - 3}}$ for the angle of incidence θ = 80° and the laser wavelength λ0 = 1.053 µm (${n_{{\rm cr}}} = {\rm \omega}_0^2 {m_{\rm e}}/4{\rm \pi} {e^2}$ is the critical density for the laser frequency ω0 = 2πc/λ0). Here and hereinafter the electron plasma density is in the centimetre-gram-second system of units, while the scale length is in micrometers. The sum of distribution Eq. (1) and two high gradient exponential density profiles with the scale lengths of 0.45 and 0.05 µm at the densities 1022 and 1023 cm−3, respectively, reproduces quite well the total density profile of plasma corona obtained in simulation (see empty green circles in Fig. 2).

Fig. 2. The electron density distribution at t = −5 ps (before the main pulse) by the action of the picosecond prepulse only. Thick (black) curve is the results of modeling; dash-dot (red) curve is the rarefaction wave exponent Eq. (1); empty (green) circles show a three-exponent approximation. Laser pulse parameters are indicated in the text.

This three-exponential electron density profile of the plasma produced by the prepulse was used for modeling of the electron acceleration by the main picosecond laser pulse with a high contrast nanosecond ASE pedestal (of the level ≤10−11 when its influence can be neglected and only the picosecond prepulse can produce plasma before the main pulse action).

2.2. Nanosecond ASE pedestals

The influence of the 1 ns ASE pedestal (which acts before −100 ps, see Fig. 1) on the preplasma electron density distribution was investigated for the pedestal intensities in the range of 109–4 × 1011 W/cm2 for the same main laser pulse parameters (P-polarization, angle of incidence 80°, and maximum intensity 4 × 1019 W/cm2) that corresponds to the pedestal contrast levels in the range 2.5 × 10−11–10−8. All over these pedestal intensities, for the electron densities less and about the reflection point y r [${n_{\rm e}}({y_{\rm r}}) = {n_{{\rm cr}}}\mathop {\cos} \nolimits^2 {\rm \theta} \approx 3 \times {10^{19}}\;{\rm c}{{\rm m}^{ - 3}}$] the density profile can be well approximated by the rarefaction wave exponent (1) with the characteristic scale length L r in the range of 2–40 µm (see Fig. 3). The reflection point position y r demonstrates the plasma propagation as a function of the pedestal intensity (see Fig. 3).

Fig. 3. The characteristic scale length of the prepulse produced plasma L r at the reflection point at t = −5 ps (before the main pulse) taking into account ASE pedestal of the length 1 ns, as a function of the laser pulse pedestal intensity – solid (blue) curve with empty circles. Circles are the result of approximation (1) of the plasma expansion, the line is the B-spline interpolation of the points. P-polarized laser pulse with the angle of incidence 80° is considered. The reflection point position corresponding to 0.03 × n cr – solid (red) curve with empty diamonds.

These modeling results evidently show that even at a rather high contrast of the laser pulse pedestal of the order of 10−10 and a nanosecond duration, the characteristic scale length of the produced preplasma can be increased two times (L r = 3.6 µm) before the main short laser pulse action. For a lower contrast of the laser pulse pedestal in the range 10−9–10−8 (and also for smaller angles of laser incidence), the characteristic scale length of the prepulse produced plasma can increase substantially and reach tens of microns as shown in Figure 3.

The exponential electron density profiles of the prepulse produced plasma with a characteristic scale length from L r ≈ 2÷4 to 30 µm near the reflection point were used for PIC modeling of the electron acceleration by the main subpicosecond laser pulse with different contrast levels of the nanosecond ASE pedestal.

3. LASER-PLASMA ELECTRON ACCELERATION

PIC simulations were carried out with the 3D PIC code Virtual Laser Plasma Laboratory  (Pukhov, Reference Pukhov1999). The geometry of interaction is shown in Figure 4.

Fig. 4. Geometry used in 3D PIC simulations: y = 0 corresponds to the initial position of the solid target surface. The laser pulse is incident on plasma, which is spread out to y > 0 with an exponential electron density profile produced by the laser prepulse (see Figs 2 and 3).

The density profiles obtained by the wide-range hydro modeling of the preplasma expansion under the action of the laser prepulse were attached at y ≥ 0 to the homogeneous plasma slab with the electron density of 100 × n cr and thickness of 5 µm over the OY-axis [y ∈ ( − 5, 0) μm, see Fig. 4]. The simulation box had the sizes of 400 µm along the OX-axis and 80 µm along the OZ-axis. The size along the OY-axis was varied from 90 to 230 µm for various preplasma density profiles to get the reflection point y r inside the simulation box far enough from the boundaries. Cells with sizes 0.05 × 0.5 × 0.5 µm3 and four particles per cell for the electrons were used together with immobile ions which produced neutralizing background.

The main P-polarized laser pulse with the Gaussian envelope in time and space, λ0 = 1 µm, the full width at half maximum (FWHM) pulse duration 400 fs and energy 120 J was focused in the FWHM spot size of 25 µm with the angle of incidence 80° and maximum intensity 4 × 1019 W/cm2. Initially, at t = 0, the center of the laser pulse was far enough from the target at x 0 = −200 µm, z 0 = 40 µm, and y 0 was varied so that the pulse reached the reflection point y r at the middle of the x-length of the simulation box (at x ≈ 200 µm).

Figure 5 shows the energy spectra of accelerated electrons for the initial density scale lengths L r = 1.8 and 3.6 µm, typical for the contrast of the laser pulse nanosecond pedestal from 10−11 to 10−10. The spectrum with a hot temperature T h = 2.3 MeV predicted by the ponderomotive electron energy at a peak laser intensity (Wilks et al., Reference Wilks, Kruer, Tabak and Langdon1992) is shown also by the dotted line in Figure 5 for the coefficient of transformation of the laser energy to the hot electrons energy η = 0.07. With the increase in L r within a few micrometers the number of accelerated electrons drops for the energies higher than 30 MeV.

Fig. 5. Electron energy spectra for the initial density scale length L r = 1.8 µm – solid (black) curve and L r = 3.6 µm – dashed (blue) curve. Dotted (red) curve shows the spectrum for a hot temperature T h = 2.3 MeV predicted by the ponderomotive electron energy.

The lower contrast nanosecond ASE pedestal of the laser pulse can create longer density scales of the preplasma with L r ≥ 10 µm (see Fig. 3). In this case, the accelerating mechanism changes and the maximum energy of accelerated electrons can reach hundreds of MeV. The dependences of spectra of accelerated electrons on the laser pulse polarization and the angle of incidence are illustrated in Figure 6 for the fixed initial density scale length L r = 10 µm.

Fig. 6. Electron energy spectra for the initial density scale length Lr = 10 µm at different laser pulse polarizations and angles of incidence. Solid (red) and dashed (green) curves are for the P-polarization, θ = 80 and 70°, respectively; dotted (blue) curve shows the spectrum for the S-polarization, θ = 80°.

With the increase of the characteristic density scale length, the maximum energy of accelerated electrons grows up to about 300 MeV for the density scale length L r = 20 µm (at the expense of some decrease of the number of particles for the energies less than 50 MeV in comparison with a longer density scale length 30 µm, see Fig. 7).

Fig. 7. Electron energy spectra for the P-polarized laser pulse with the angle of incidence θ = 80 and different initial density scale lengths: Solid (red) curve for L r = 20 µm and dashed (blue) curve for L r = 30 µm.

Figure 8 illustrates qualitatively different angular distributions of accelerated electrons for different scale lengths formed by different contrast laser prepulses. For the preplasma scale length L r = 20 µm, which can be produced by the nanosecond ASE pedestal of an intensity about 1011 W/cm2 (see Fig. 3), accelerated electrons propagate at a small angle to the target surface, and the width of the angular distribution decreases for higher electron energies (upper plots of Fig. 8). While for a lower level of the laser pulse nanosecond ASE pedestal (of intensity less than 109 W/cm2, when the preplasma scale length L r = 1.8 µm), there are two jets of accelerated electrons propagating outward and inside the target with the angle (relative to the target surface) about the angle of propagation of the reflected laser pulse (lower plots of Fig. 8).

Fig. 8. Angular distributions of accelerated electrons with energies > 50 MeV [left, (a) and (c)] and >100 MeV [right, (b) and (d)] for preplasma scale lengths L r = 20 [upper plots, (a) and (b)] and 1.8 µm [lower plots, (c) and (d)]. The laser pulse is P-polarized (along the y-axis) with the angle of incidence θ = 80°, ${\rm \varphi} = \arctan ({\,p_y}/{\,p_x})$, where p x and p y are the electron momentum components.

4. DISCUSSION AND CONCLUSIONS

Different energy spectra and angle distributions of hot electrons at ten times different preplasma scale lengths L r = 1.8 and 20 µm (see Figs 5, 7, and 8) indicate different mechanisms of the electron acceleration in the plasma corona. In the case of sharp (short scale length) plasma density gradients, produced quasi-static electric and magnetic fields are significant for generation and transport of high energy electrons along the surface of a solid target irradiated by high power laser pulses. The importance of quasi-static fields in the electron acceleration by a self-channeled relativistic-intense laser pulse was shown in paper (Pukhov et al., Reference Pukhov, Sheng and Meyer-ter Vehn1999). The mechanism and conditions for the formation of quasi-static surface fields at oblique incidence of a short intense laser pulse onto a solid target with a short scale length preplasma corona were analyzed in papers (Nakamura et al., Reference Nakamura, Kato, Nagatomo and Mima2004, Reference Nakamura, Mima, Sakagami and Johzaki2007). In full agreement with 2D PIC simulation results of (Nakamura et al., Reference Nakamura, Mima, Sakagami and Johzaki2007) (see Fig. 3 of this paper), the structure of quasi-static electric and magnetic fields obtained in our 3D PIC modeling (shown in Figure 9 for the preplasma scale length L r = 1.8 µm) indicates the formation of the focusing electric and magnetic fields reflecting electrons towards the vacuum. It should be noted that the enhanced transverse momentum of accelerated electrons (see Fig. 8 for the preplasma scale lengths L r = 1.8 µm), which exceeds substantially the normalized vector potential of the laser pulse, clearly indicates the resonant direct laser acceleration with assistance of quasi-static electric and magnetic fields (Pukhov et al., Reference Pukhov, Sheng and Meyer-ter Vehn1999).

Fig. 9. Averaged over laser period electric field <e y> ≡ <|e| E y/mcω0> – (a) and magnetic field <b z> ≡ <|e|B z/mcω0> – (b) on the plane OXY at z = 40 µm, ct = 200 µm (when the laser pulse center is located at x = 200 µm) for the preplasma scale length L r = 1.8 µm.

A relatively low contrast nanosecond pedestal of the laser pulse can produce longer preplasma corona extended to tens of microns with density scales of the order or in excess of the laser spot size. In this case the mechanism of the laser acceleration along the target surface can be attributed to stochastic electron acceleration in plasma waves driven by an intense laser pulse at grazing incidence [see (Bochkarev et al., Reference Bochkarev, Brantov, Bychenkov, Torshin, Kovalev, Baidin and Lykov2014) and references cited therein]. Plasma waves in subcritical plasma are excited by a subpicosecond laser pulse with the duration much longer than the plasma wave period in the process of the pulse self-modulation (Antonsen & Mora, Reference Antonsen and M. Mora1993; Krall et al., Reference Krall, Ting, Esarey and Sprangle1993; Andreev et al., Reference Andreev, Gorbunov, Kirsanov, Pogosova and Sakharov1996, Reference Andreev, Gorbunov, Kirsanov, Pogosova and Ramazashvili1992). The self-modulation instability develops from the front of the laser pulse and saturates due to the plasma wave breaking. Fourier spectra of the longitudinal electric field e x ≡ |e| E x/mcω0 in the regions x ∈ (0, 200) μm and x ∈ (200, 400) μm, where the rear and front parts of the laser pulse are situated at ct = 200 µm, are shown in Figure 10(a) and 10(b), respectively, for the preplasma scale length L r = 20 µm at different distances from the initial target surface.

Fig. 10. Absolute value of the longitudinal electric field e x ≡ |e| E x/mcω0 Fourier-component at time ct = 200 µm (when the laser pulse center is located at x = 200 µm) calculated for different intervals over OX: x ∈ (0, 200) μm – (a) and x ∈ (200, 400) μm – (b) for the different distances y from the initial target surface (y = 0): y = 12 µm – solid (red) curves, y = 20 µm – dashed (blue) curves, y = 50 µm – dotted (green) curves. Results presented for the preplasma scale length L r = 20 µm and the laser incidence angle 80°.

The enhanced Fourier-component on the right plot of Figure 10 in the spectral range $k/{k_0} \cong 0.1 \div 0.2$ (k 0 = ω0/c) reflects the excitation of plasma waves at plasma densities ${n_{\rm e}}/{n_{{\rm cr}}} \cong (1 \div 3) \times {10^{ - 2}}$ that is just about and below the unperturbed laser reflection point (y r ≈ 20 µm) for the angle of incidence 80°. The absence of these spectral lines in the left plot of Figure 10 (for the rear part of the laser pulse) indicates the saturation of self-modulation instability and transmission of the excited plasma waves energy to accelerated electrons.

Angular distributions of fast electrons emitted from intense laser–solid interactions were investigated in experiments at the PHELIX laser system at Gesellschaft fur Schwerionenforschung, Darmstadt for different angles of laser incidence (Gray et al., Reference Gray, Yuan, Carroll, Brenner, Coury, Quinn, Tresca, Zielbauer, Aurand, Bagnoud, Fils, KГјhl, Lin, Li, Li, Roth, Neely and McKenna2011). In these experiments the intensity contrast up to 50 ps prior to the peak of the pulse was about 10−6, that corresponds to the long scale plasma corona of tens of microns (see Fig. 3, L r). The measured angular distributions of energetic electrons for the laser incidence angle 80° is in good agreement with the modeling result shown in the upper plots of Figure 8.

In conclusion, preplasma formation under the action of different contrast prepulses is modeled and described using the wide-range two-temperature hydrodynamic model. The preplasma density profiles for different contrast ratios of the nanosecond pedestal are found and used as the initial density distributions in 3D PIC simulations of electron acceleration by the main P-polarized laser pulse.

The generation of hot electrons at the grazing incidence of a subpicosecond relativistic-intense laser pulse onto the plane solid target is analyzed for the parameters typical for petawatt class laser systems. Table 1 summarizes the parameters of electrons accelerated to high energies for two typical scale lengths of the preplasma density L r = 1.8 and 20 µm. For the grazing incidence of the laser pulse, the substantial increase of the characteristic energy, number, and collimation of electrons accelerated along the target surface is demonstrated (see Figs 5–7) in comparison with the ponderomotive scaling of laser–target interaction (Wilks et al., Reference Wilks, Kruer, Tabak and Langdon1992; Chen et al., Reference Chen, Wilks, Kruer, Patel and Shepherd2009).

Table 1. The parameters of electrons accelerated to high energies for two typical scale lengths of the preplasma density.

ACKNOWLEDGMENTS

This work was supported by the Russian Science Foundation grant 14-50-00124.

References

REFERENCES

Andreev, N.E., Gorbunov, L.M., Kirsanov, V., Pogosova, A.A. & Ramazashvili, R.R. (1992). Resonant excitation of wakefields by a laser pulse in a plasma. JETP Lett. 55, 6.Google Scholar
Andreev, N.E., Gorbunov, L.M., Kirsanov, V.I., Pogosova, A.A. & Sakharov, A.S. (1996). Theory of the resonance modulational instability of short laser pulses in a homogeneous plasma and plasma channels. Plasma Phys. Rep. 22, 379–389.Google Scholar
Andreev, N.E., Povarnitsyn, M.E., Veysman, M.E., Faenov, A.Y., Levashov, P.R., Khishchenko, K.V., Pikuz, T.A., Magunov, A.I., Rosmej, O.N., Blazevic, A., Pelka, A., Schaumann, G., Schollmeier, M. & Roth, M. (2015). Interaction of annular-focused laser beams with solid targets. Laser Part. Beams 33, 541–550.CrossRefGoogle Scholar
Antonsen, T. Jr & M. Mora, P. (1993). Self-focusing and Raman scattering of laser pulses in tenuous plasmas. Phys. Fluids B 5, 1440–1452.CrossRefGoogle Scholar
Bagnoud, V., Aurand, B., Blazevic, A., Borneis, S., Bruske, C., Ecker, B., Eisenbarth, U., Fils, J., Frank, A., Gaul, E., Goette, S., Haefner, C., Hahn, T., Harres, K., Heuck, H.-M., Hochhaus, D., Hoffmann, D.H.H., Javorková, D., Kluge, H.-J., Kuehl, T., Kunzer, S., Kreutz, M., Merz-Mantwill, T., Neumayer, P., Onkels, E., Reemts, D., Rosmej, O., Roth, M., Stoehlker, T., Tauschwitz, A., Zielbauer, B., Zimmer, D. & Witte, K. (2010). Commissioning and early experiments of the phelix facility. Appl. Phys. B 100, 137–150.CrossRefGoogle Scholar
Bochkarev, S.G., Brantov, A.V., Bychenkov, V.Y., Torshin, D.V., Kovalev, V.F., Baidin, G.V. & Lykov, V.A. (2014). Stochastic electron acceleration in plasma waves driven by a high-power subpicosecond laser pulse. Plasma Phys. Rep. 40, 202–214.CrossRefGoogle Scholar
Brabetz, C., Busold, S., Cowan, T., Deppert, O., Jahn, D., Kester, O., Roth, M., Schumacher, D. & Bagnoud, V. (2015). Laser-driven ion acceleration with hollow laser beams. Phys. Plasmas 22, 013105.CrossRefGoogle Scholar
Chen, H., Wilks, S.C., Kruer, W.L., Patel, P.K. & Shepherd, R. (2009). Hot electron energy distributions from ultraintense laser solid interactions. Phys. Plasmas 16, 020705.CrossRefGoogle Scholar
Esirkepov, T.Z., Koga, J.K., Sunahara, A., Morita, T., Nishikino, M., Kageyama, K., Nagatomo, H., Nishihara, K., Sagisaka, A., Kotaki, H., Nakamura, T., Fukuda, Y., Okada, H., Pirozhkov, A.S., Yogo, A., Nishiuchi, M., Kiriyama, H., Kondo, K., Kando, M. & Bulanov, S.V. (2014). Prepulse and amplified spontaneous emission effects on the interaction of a petawatt class laser with thin solid targets. Nucl. Instrum. Methods A 745, 150–163.CrossRefGoogle Scholar
Gray, R.J., Yuan, X.H., Carroll, D.C., Brenner, C.M., Coury, M., Quinn, M.N., Tresca, O., Zielbauer, B., Aurand, B., Bagnoud, V., Fils, J., KГјhl, T., Lin, X.X., Li, C., Li, Y.T., Roth, M., Neely, D. & McKenna, P. (2011). Surface transport of energetic electrons in intense picosecond laser-foil interactions. Appl. Phys. Lett. 99, 171502.CrossRefGoogle Scholar
Krall, J., Ting, A., Esarey, E. & Sprangle, P. (1993). Enhanced acceleration in a self-modulated-laser wake-field accelerator. Phys. Rev. E 48, 2157.CrossRefGoogle Scholar
Mao, J.Y., Chen, L.M., Huang, K., Ma, Y., Zhao, J.R., Li, D.Z., Yan, W.C., Ma, J.L., Aeschlimann, M., Wei, Z.Y. & Zhang, J. (2015). Highly collimated monoenergetic target-surface electron acceleration in near-critical-density plasmas. Appl. Phys. Lett. 106, 131105.CrossRefGoogle Scholar
Morace, A., Fedeli, L., Batani, D., Baton, S., Beg, F.N., Hulin, S., Jarrott, L.C., Margarit, A., Nakai, M., Nakatsutsumi, M., Nicolai, P., Piovella, N., Wei, M.S., Vaisseau, X., Volpe, L. & Santos, J.J. (2014). Development of X-ray radiography for high energy density physics. Phys. Plasmas 21, 102712.CrossRefGoogle Scholar
Nakamura, T., Kato, S., Nagatomo, H. & Mima, K. (2004). Surface-magnetic-field and fast-electron current-layer formation by ultraintense laser irradiation. Phys. Rev. Lett. 93, 265002.CrossRefGoogle ScholarPubMed
Nakamura, T., Mima, K., Sakagami, H. & Johzaki, T. (2007). Electron surface acceleration on a solid capillary target inner wall irradiated with ultraintense laser pulses. Phys. Plasmas 14, 053112.CrossRefGoogle Scholar
Ortner, A., Faik, S., Schumacher, D., Basko, M.M., Blazevic, A., Busold, S., Bedacht, S., Cayzac, W., Frank, A., Kraus, D., Rienecker, T., Schaumann, G., Tauschwitz, A., Wagner, F. & Roth, M. (2015). A novel double hohlraum target to create a moderately coupled plasma for ion stopping experiments. Nucl. Instrum. Methods B 343, 123–131.CrossRefGoogle Scholar
Povarnitsyn, M.E., Andreev, N.E., Apfelbaum, E.M., Itina, T.E., Khishchenko, K.V., Kostenko, O.F., Levashov, P.R. & Veysman, M.E. (2012 a). A wide-range model for simulation of pump-probe experiments with metals. Appl. Surf. Sci. 258, 9480–9483.CrossRefGoogle Scholar
Povarnitsyn, M.E., Andreev, N.E., Levashov, P.R., Khishchenko, K.V., Kim, D.A., Novikov, V.G. & Rosmej, O.N. (2013). Laser irradiation of thin films: Effect of energy transformation. Laser Part. Beams 31, 663–671.CrossRefGoogle Scholar
Povarnitsyn, M.E., Andreev, N.E., Levashov, P.R., Khishchenko, K.V. & Rosmej, O.N. (2012 b). Dynamics of thin metal foils irradiated by moderate-contrast high-intensity laser beams. Phys. Plasmas 19, 023110.CrossRefGoogle Scholar
Pukhov, A. (1999). Three-dimensional electromagnetic relativistic particle-in-cell code VLPL (Virtual Laser Plasma Lab). J. Plasma Phys. 61, 425–433.CrossRefGoogle Scholar
Pukhov, A., Sheng, Z.-M. & Meyer-ter Vehn, J. (1999). Particle acceleration in relativistic laser channels. Phys. Plasmas 6, 2847–2854.CrossRefGoogle Scholar
Rusby, D.R., Wilson, L.A., Gray, R.J., Dance, R.J., Butler, N.M.H., MacLellan, D.A., Scott, G.G., Bagnoud, V., Zielbauer, B., McKenna, P. & Neely, D. (2015). Measurement of the angle, temperature and flux of fast electrons emitted from intense laser–solid interactions. J. Plasma Phys. 81, 475810505.CrossRefGoogle Scholar
Schollmeier, M., Sefkow, A.B., Geissel, M., Arefiev, A.V., Flippo, K.A., Gaillard, S.A., Johnson, R.P., Kimmel, M.W., Offermann, D.T., Rambo, P.K., Schwarz, J. & Shimada, T. (2015). Laser-to-hot-electron conversion limitations in relativistic laser matter interactions due to multi-picosecond dynamics. Phys. Plasmas 22, 043116.CrossRefGoogle Scholar
Vaisseau, X., Debayle, A., Honrubia, J.J., Hulin, S., Morace, A., Nicola, P., Sawada, H., Vauzour, B., Batani, D., Beg, F.N., Davies, J.R., Fedosejevs, R., Gray, R.J., Kemp, G.E., Kerr, S., Li, K., Link, A., McKenna, P., McLean, H.S., Mo, M., Patel, P.K., Park, J., Peebles, J., Rhee, Y.J., Sorokovikova, A., Tikhonchuk, V.T., Volpe, L., Wei, M. & Santos, J.J. (2015). Enhanced relativistic-electron-beam energy loss in warm dense aluminum. Phys. Rev. Lett. 114, 095004.CrossRefGoogle ScholarPubMed
Wagner, F., Bedacht, S., Ortner, A., Roth, M., Tauschwitz, A., Zielbauer, B. & Bagnoud, V. (2014 a). Pre-plasma formation in experiments using petawatt lasers. Opt. Express 22, 29505–29514.CrossRefGoogle ScholarPubMed
Wagner, F., João, C., Fils, J., Gottschall, T., Hein, J., Körner, J., Limpert, J., Roth, M., Stöhlker, T. & Bagnoud, V. (2014 b). Temporal contrast control at the phelix petawatt laser facility by means of tunable sub-picosecond optical parametric amplification. Appl. Phys. B 116, 429–435.CrossRefGoogle Scholar
Wang, W.M., Chen, L.M., Mao, J.Y., Huang, K., Ma, Y., Zhao, J.R., Zhang, L., Yan, W.C., Li, D.Z., Ma, J.L., Ly, Y.T., Lu, X., Wei, Z.Y., Sheng, Z.M. & Zhang, J. (2013). Collimated quasi-monoenergetic electron beam generation from intense laser solid interaction. High Energy Density Phys. 9, 578–582.CrossRefGoogle Scholar
Wilks, S.C., Kruer, W.L., Tabak, M. & Langdon, A.B. (1992). Absorption of ultra-intense laser pulses. Phys. Rev. Lett. 69, 1383.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Time structure of the PHELIX laser pulse at the first harmonic, λ0 = 1.053 µm, in a subnanosecond domain for high contrast level ~10−10 at time less than −100 ps. Solid (black) curve presents the cross-correlation measurement, dashed (red) curve is an approximation used in simulations. Laser intensity is normalized to the maximum of the main pulse intensity.

Figure 1

Fig. 2. The electron density distribution at t = −5 ps (before the main pulse) by the action of the picosecond prepulse only. Thick (black) curve is the results of modeling; dash-dot (red) curve is the rarefaction wave exponent Eq. (1); empty (green) circles show a three-exponent approximation. Laser pulse parameters are indicated in the text.

Figure 2

Fig. 3. The characteristic scale length of the prepulse produced plasma Lr at the reflection point at t = −5 ps (before the main pulse) taking into account ASE pedestal of the length 1 ns, as a function of the laser pulse pedestal intensity – solid (blue) curve with empty circles. Circles are the result of approximation (1) of the plasma expansion, the line is the B-spline interpolation of the points. P-polarized laser pulse with the angle of incidence 80° is considered. The reflection point position corresponding to 0.03 × ncr – solid (red) curve with empty diamonds.

Figure 3

Fig. 4. Geometry used in 3D PIC simulations: y = 0 corresponds to the initial position of the solid target surface. The laser pulse is incident on plasma, which is spread out to y > 0 with an exponential electron density profile produced by the laser prepulse (see Figs 2 and 3).

Figure 4

Fig. 5. Electron energy spectra for the initial density scale length Lr = 1.8 µm – solid (black) curve and Lr = 3.6 µm – dashed (blue) curve. Dotted (red) curve shows the spectrum for a hot temperature Th = 2.3 MeV predicted by the ponderomotive electron energy.

Figure 5

Fig. 6. Electron energy spectra for the initial density scale length Lr = 10 µm at different laser pulse polarizations and angles of incidence. Solid (red) and dashed (green) curves are for the P-polarization, θ = 80 and 70°, respectively; dotted (blue) curve shows the spectrum for the S-polarization, θ = 80°.

Figure 6

Fig. 7. Electron energy spectra for the P-polarized laser pulse with the angle of incidence θ = 80 and different initial density scale lengths: Solid (red) curve for Lr = 20 µm and dashed (blue) curve for Lr = 30 µm.

Figure 7

Fig. 8. Angular distributions of accelerated electrons with energies > 50 MeV [left, (a) and (c)] and >100 MeV [right, (b) and (d)] for preplasma scale lengths Lr = 20 [upper plots, (a) and (b)] and 1.8 µm [lower plots, (c) and (d)]. The laser pulse is P-polarized (along the y-axis) with the angle of incidence θ = 80°, ${\rm \varphi} = \arctan ({\,p_y}/{\,p_x})$, where px and py are the electron momentum components.

Figure 8

Fig. 9. Averaged over laser period electric field <ey> ≡ <|e| Ey/mcω0> – (a) and magnetic field <bz> ≡ <|e|Bz/mcω0> – (b) on the plane OXY at z = 40 µm, ct = 200 µm (when the laser pulse center is located at x = 200 µm) for the preplasma scale length Lr = 1.8 µm.

Figure 9

Fig. 10. Absolute value of the longitudinal electric field ex ≡ |e| Ex/mcω0 Fourier-component at time ct = 200 µm (when the laser pulse center is located at x = 200 µm) calculated for different intervals over OX: x ∈ (0, 200) μm – (a) and x ∈ (200, 400) μm – (b) for the different distances y from the initial target surface (y = 0): y = 12 µm – solid (red) curves, y = 20 µm – dashed (blue) curves, y = 50 µm – dotted (green) curves. Results presented for the preplasma scale length Lr = 20 µm and the laser incidence angle 80°.

Figure 10

Table 1. The parameters of electrons accelerated to high energies for two typical scale lengths of the preplasma density.