Towards bright gamma-ray flash generation from tailored target irradiated by multi-petawatt laser

One of the remarkable phenomena in the laser-matter interaction is the extremely efficient energy transfer to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}γ-photons, that appears as a collimated \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}γ-ray beam. For interactions of realistic laser pulses with matter, existence of an amplified spontaneous emission pedestal plays a crucial role, since it hits the target prior to the main pulse arrival, leading to a cloud of preplasma and drilling a narrow channel inside the target. These effects significantly alter the process of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}γ-photon generation. Here, we study this process by importing the outcome of magnetohydrodynamic simulations of the pedestal-target interaction into particle-in-cell simulations for describing the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}γ-photon generation. It is seen that target tailoring prior the laser-target interaction plays an important positive role, enhancing the efficiency of laser pulse coupling with the target, and generating high energy electron-positron pairs. It is expected that such a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}γ-photon source will be actively used in various applications in nuclear photonics, material science and astrophysical processes modelling.

S ince the invention of the Chirped Pulse Amplification technique (1), development of multi-petawatt (multi-PW) laser systems is envisioned worldwide.The first laser to reach the 10 PW peak power level has been recently reported by ELI-NP (2), Romania, with a pulse duration of ∼ 25 fs.Another ∼ 10 PW laser is soon expected in ELI-Beamlines (3), Czech Republic, with approximately five times higher energy.In addition, ELI-ALPS, Hungary, aims at constructing an ultrashort ∼ 17 fs laser of ∼ 2 PW (4).Moreover, a laser combining pulses shorter than 20 fs in the 10 PW level is developed in Apollon facility (5), France.By focusing a multi-PW laser down to a micrometer-wide spot (6,7), intensities exceeding 10 27 W m −2 can be achieved, where this threshold has been recently surpassed by the ∼ 4 PW CoReLS laser (7), South Korea.
A typical high power laser consists of an ultrashort main pulse, preceded by a lower amplitude background field extending in the nanosecond scale.The laser contrast is defined as the ratio of the main pulse amplitude to the background field amplitude.Usually, in high power laser systems the contrast is increased through complex and/or expensive additions, such as Optical Parametric Chirped Pulse Amplification (8) and plasma mirrors (9,10).
The importance of a finite contrast for laser-matter interactions is highlighted both experimentally (11)(12)(13)(14) and theoretically (12,13,(15)(16)(17).By assuming an initially steep density gradient (flat-foil, or simply foil) target, in the aforementioned literature it is agreed that the background field modifies the initial density profile at an extent proportional to its amplitude and duration.A relatively thick (micrometer range) foil is curved in the vicinity of the laser focal spot, where a gradually increasing density profile appears in the target front region.On the other hand, if the target is thin enough then the background field drills the target resulting in no interaction when the main pulse arrives.As a result, the preplasma strongly affects the energy spectra and directionality of particles emitted due to the laser-target interaction.However, since computational studies of the interaction typically involve particle-in-cell (PIC) simulations which cannot be applied to model the nanosecond long duration required by the background field, usually a foil target is assumed, acknowledging only the effect of the main pulse on target.
A plethora of exotic target geometries has been already considered for laser-matter interaction experimnets.Among numerous examples, we mention proton-rich micro-dots (18), cylindircal micro-lenses (19,20), hollow micro-spheres (21), micro-coils (22) and wavelength-scale holed targets (23,24).All of the aforementioned target designs require explicit microfabrication techniques, while they add further complexity to a laser-target experiment since they require additional efforts on positioning and alignment of the target.However, it was noticed that the use of tailored targets is favorable for the laser-target interaction and their use is widely employed.
Apart of target tailoring methods outside the interaction chamber, it is also possible to manipulate the electron density profile by means of a secondary, long pulse duration, lower power laser (25).The secondary laser provides the required background field, where temporal control of the laser (26) can provide the required electron density distributions.This

Significance Statement
Gamma-ray flashes are phenomena of diachronic interest that occur not only in the microcosm but also in cosmic environments.In the laboratory, extreme gamma-flashes may be generated during laser-matter interaction.In finding optimal conditions of γ-photon generation, it is crucial to elucidate the role of laser pulse structure comprising of ultrashort main pulse, preceded by a relatively long background field that alters the target profile.Background field removal not only adds complexity to high power lasers but also comes at cost of reduced laser energy.In the present work, we find an approach for substantial enhancement of the laser energy conversion to γ-photon energy.

D R A F T
scheme is known as 'laser heater' and has been successfully implemented on improving the properties of laser driven particle beams (27,28).Notably, manipulating the target density with a laser heater resulted in record energy values of particles accelerated by optical means (29).
At the ultrahigh intensity limit, the theory foresees that the multi-PW lasers interacting with matter will provide not only energetic charged particles (electrons (e − ), protons (p + ), heavy ions (i + )) but also a plethora of high energy γ-photons and electron-positron (e − -e + ) pairs.Once an electron (or positron) collides with the incident field it is scattered, resulting in alteration of its momentum and a γ-photon is emitted, in a process known as multiphoton Compton scattering.In the present work, we focus on γ-photons produced via the multiphoton Compton scattering dominating at ultrahigh intensities (30)(31)(32)(33), produced within a time approximately equal to to the laser pulse duration.Notably, at lower intensities (< 10 27 W m −2 ) and significantly thick targets γ-photons can be also produced via Bremsstrahlung emission (34), where the present target geometry limits Bremsstrahlung contribution.The γ-photon production then makes possible the production of e − -e + pairs via the multiphoton Breit-Wheeler process (35), where a high energy γ-photon interacts with multiple laser photons.Maximizing the laser to γ-photon energy conversion efficiency, κγ, is required in photonuclear physics (36,37), for study of high energy density physics in materials science (38) and studies of astrophysical processes (39)(40)(41)(42).
The recently available multi-PW lasers in combination with the broad application range of the γ-photons, draw the interest of several research groups, suggesting various methods to maximize κγ.An early suggestion is based on the reflected part of the laser pulse incident on an overdense target (43)(44)(45).However, loss of significant laser energy towards the reflection region suggested the use of two counter-propagating pulses (46)(47)(48)(49), later extending the scheme to multiple colliding pulses (50,51).Apart from the all-optical approach, other groups suggested microfabricating sophisticated target schemes (52,53) or even combining a laser with either optically (54) or externally accelerated electrons (55).Moreover, it was shown that by adding a chosen preplasma on the target front, κγ can also be increased (32,56).

Hybrid Magnetohydrodynamic and Particle-in-Cell Simulations
In reference (57) it has been demonstrated through magnetohydrodynamic (MHD) simulations that a unique microfabrication of the foil target is achieved by a background field corresponding to multi-PW lasers.The radial symmetry of the results allows their expansion on a (x, ŷ, ẑ) grid, where the laser propagation axis is set along x.Therefore one can obtain the three-dimensional (3D) electron number density distribution.The results of the MHD simulations are then used as initial conditions in PIC simulations, in a hybrid modeling scheme that reveals the physical processes governing the complete interaction of the main laser pulse with the newly formed density distribution.
The background field sets a 60 ps duration region with intensity starting at 10 19 W m −2 , dropping exponentially to 10 14 W m −2 , and then continues with a constant intensity for 940 ps.The laser intensity spatial profile is represented by a Gaussian of ∼ 2.2 µm full-width-at-half-maximum.The afore-mentioned background field intensity is comparable with what was considered as cutting-edge laser main pulse intensity a few decades ago (25).From reference (57) we choose a set of elements with approximately equally increasing electron number density, ne, namely lithium (Li, ne ≈ 1.39 × 10 29 m −3 ), sodium (Na, ne ≈ 2.79 × 10 29 m −3 ), beryllium (Be, ne ≈ 4.95 × 10 29 m −3 ), carbon (C, ne ≈ 6.02 × 10 29 m −3 ), and aluminum (Al, ne ≈ 7.83 × 10 29 m −3 ).The effect of the background field on lithium is shown in Fig. 1, where lithium having the lowest density among all solids at room temperature it is affected the strongest.In agreement with previous works (16,32), a preplasma distribution is generated in the target front region, although with a large exponential coefficient for the electron number density gradient of ∼ 3.1 × 10 5 m −1 , shown by the red lineout in Fig. 1, along the laser propagation axis.Notably, the preplasma exceeding the critical density, ncr, extends for only ∼ 1 µm in the target front region.
The sizable target (12 µm diameter and 10 µm thickness) ensures neither target curvature on the target rear, nor complete drilling and destruction of the target (16,17).In addition to the preplasma formation, a conical-like cavity is generated within the target volume, with walls several times denser than the background.For lithium, the cavity has a depth of ∼ 4 µm and a base diameter of ∼ 6 µm.The Gaussian background field profile allows stronger target deformation in the center, where the cavity profile exhibits two off-center symmetric vertexes (where the derivative of the curvature is zero).This unique tailored target is achievable with finite contrast lasers currently available, revealing a preplasma regime which to the best of our knowledge has never been considered in PIC simulations.
The PIC code used is the relativistic quantum electrodynamic EPOCH code (58), compiled with the Higuera-Cary ( 59  In the present work, the default target referral is a tailored lithium target and the default laser power is ∼ 10 PW, except where it is stated otherwise.As the background field deposits most of its energy on target, a conical-like cavity is generated prior to the arrival of the main laser pulse.This scheme resembles laser-target interaction geometries on which a cone is purposely fabricated at the target, aiming at novel fast-ignition schemes (61,62), increasing laser induced γ-photon production (63), enhancing laser field intensity (64), and efficient proton acceleration (65,66).Since in our case the relativistic critical density is above lithium electron number density, one might expect the cavity formation to have no effect on focusing the laser field, due to dominance of relativistic self-focusing (67) of the laser field in the underdense solid.However, as seen by the spikes in the density lineout profiles in Fig. 1, a thin, overdense layer D R A F T is formed on the cavity walls.As a result, the laser field can penetrate only within the skin-depth of the walls, and reflected towards the cavity depth.

Cavity Propagation and Intensity Enhancement
The reflected fields then interfere, rapidly increasing the field intensity.Note that the laser focal spot is at the base of the cone, meaning that in vacuum the laser would be defocusing if the cavity did not exist.In addition, the cavity volume is filled with a low electron density, that also aids focusing due to weak relativistic self-focusing; although to a lesser extent than the cavity focusing, as its electron density increases exponentially.
Once the laser is intensified in the cavity, its intensity surpasses that expected in the focal spot.As a result, the thin overdense electron distribution becomes undercritical, and the laser breaks into the target volume (68), where the lithium density is not modified by the background field.In addition, field reflection by the cavity walls results in caustics, instantaneously increasing the field intensity.These effects are seen in Fig. 2(A), where the laser intensity is shown when reaching its maximum value of ∼ 2.6 × 10 28 W m −2 , an order of magnitude higher than the intensity expected at the focal spot.The electron number density is overplotted on the figure, visualizing the spatial location of the laser field with respect to the cavity.The highest intensity is recorded at ∼ 15 fs after the laser field reaches the focal spot.
Figs. 3(A) and 3(B) show line-outs (along laser propagation axis) of the transverse and longitudinal components of the electric field (Ey and Ex respectively) at successive times.When the laser field arrives at the cavity (at ∼ − 5 fs), Ex is practically absent.However, reflection of the laser field on the cavity walls allows not only interference of Ey but also generation of Ex, where the magnitude of the two field components is comparable since Ey ≈ 2Ex for times after ∼ 0 fs.Existence of Ex ceases only when both field components dissipate, near the end of the simulation.
Penetration of the laser field in the target at successive times is shown in Fig. 2(B), by overlaying three electron density contours.The color magnitude serves purely for visualization purposes, since appropriate contrast is applied to the figure to reveal the cavity propagation.Although the laser can penetrate the relativistically undercritical target, its propagation is further assisted by Ex, which further drills the target along its propagation axis.
As a result, the laser pulse can maintain high intensities until the end of the simulation, although dissipating due to energy transfer to electrons and γ-photons.This can be seen in Fig. 3(C ), where the left axis plots the laser peak intensity as a function of time.For comparison, the orange dashed line is the intensity in the focal spot in vacuum.At ∼ − 5 fs the laser field is in a region of extremely low electron density, also in the vicinity of the focal spot region.Therefore, the peak intensity coincides with the expected peak intensity in vacuum.Beyond that time, up to ∼ 15 fs, the intensity is increasing.At that time interval most of the field energy is transferred to the rest of particles (γ-photons, electrons, positrons, ions), shown by the dashed purple line in Fig. 3(C ).
Laser energy conversion to particle energy continues up to ∼ 35 fs, where the peak intensity drops since no significant laser energy remains within the propagating cavity.The laser field energy saturates at ∼ 15 % of its initial value, due to laser back-reflection.The κγ, shown by the solid purple line, exhibits a sigmoidal behavior with changing curvature at the time the laser intensity is maximized, saturating at ∼ 30 %. Remarkably, conversion of approximately one third of the laser energy to γ-photons with the currently available technology, is of interest for worldwide laser facilities.

Gamma-Ray Flash Scaling with Power and Target Material
The aforementioned laser-target interaction leads to a plethora of energetic electrons in the GeV-scale, which in turn result in ion acceleration with energies per nucleon of ∼ 300 MeV as seen in Fig. 4(A).We should clarify that optimizing ion acceleration is not among the goals of the present work, hence the thick target chosen.By considering that for optimal laser coupling to electrons (69) the condition a0 = π(ne/ncr)(l/λ) must be satisfied (where a0 is the dimensionless amplitude and l is the penetration depth per wavelength), it can be calculated that the intensity ranges shown in Fig. 3(C ) correspond to 1 µm < l < 3 µm.Since the temporal pulse length is ∼6 λ, it is estimated that a ∼ 10 µm thick target is required for the laser pulse to efficiently transfer its energy to the target.

D R A F T
scattering.However, the micrometer-thick, low atomic number target is unsuitable for Bremsstrahlung emission, although the emission continues for a significantly longer time than the pulse duration.In addition, the stochastic dynamics of electrons in the solid target results in a non-directional Bremsstrahlung emission, making γ-photons from Compton scattering dominate the γ-photon signal during ultraintense laser-matter experiments.The Compton γ-photon spectrum is comparable to the electron spectrum, as seen in Fig. 4(A), with a temperature of ∼ 55 MeV for γ-photons of > 200 MeV energy.For comparison, the Bremsstrahlung spectrum at the end of the simulation is also shown, which is several orders of magnitude lower than the Compton spectrum.However, its contribution can only increase the cumulative γ-photon energy, with a κγ of ∼ 30 % only from Compton scattering contribution.Hence, the default γ-photon referral is Compton scattering.
At high intensity laser-matter interactions, generation of e − -e + pairs is possible through the multiphoton Breit-Wheeler process.Although e − -e + pair generation is not expected to be observed for a ∼ 10 PW laser (corresponding to intensities of ∼2.8 10 27 W m −2 ), we demonstrate that the one order of magnitude intensity enhancement due to the tailored target makes multiphoton Breit-Wheeler e − -e + pair observation feasible.The positron spectrum is shown in Fig. 4(A), spanning up to the maximum γ-photon energy with a Maxwell-Juttner distribution of ∼ 80 MeV temperature, where ∼ 4 × 10 8 pairs are recorded.
The preformed cavity importance on κγ optimization and γ-photon maximum energy increase is seen in Fig. 4(B) by comparing the solid lines (tailored target) with the dashed lines (foil target) for various laser powers.In all cases, the cavity formation significantly amplifies the γ-photon spectrum, where the exponential behavior still persists but with decreasing temperature for decreased power.The quantitative results for κγ as a function of a0 are shown in Fig. 5(A), where the blue line corresponds to tailored targets and the green to foil targets.The figure reveals an almost linear dependency of κγ to a0 (or to the laser power) within the range of interest.
The right axis of Fig. 5(A) shows the ratio of κγ from tailored to that of foil lithium targets as a function of a0.Notably, for a ∼ 10 PW laser the ratio is ∼ 1.5, but rising to ∼ 10 for a ∼ 1 PW laser.This result is related to the relativistic transparency of lithium for the ∼ 10 PW case, where even in the case of no cavity formation the laser pulse is still efficiently penetrating the target.On the other hand, the ∼ 1 PW case relies on intensity enhancement within the cavity to reach the relativistic transparency threshold, whilst the foil target strongly reflects the laser pulse with little-to-no conversion to γ-photons.Therefore, even single-PW lasers can be used to efficiently create a γ-ray flash in the laboratory, when a proper background field is combined with a relatively low density target.
In connection to the high κγ for lithium targets, the results for denser materials are shown in the columns of 5(B), where the rows correspond to varying laser power.Study of the dependency of κγ on the target material is performed only for the ∼ 10 PW case, since in all lower powers it follows a similar trend.For the materials considered, it is found that κγ is inversely proportional to ne.By increasing the target density, two effects result in reduced κγ.Firstly, the cavity formation is less prominent in a denser material, resulting in lower intensity amplification.Secondly, the target electron number density shifts away from the relativistic critical density and the laser field can no longer be efficiently coupled to the target electrons.Notably, materials commonly used in experiments have an electron density of a few times higher than that of lithium (e.g.approximately five times for aluminum).As a result, if aluminum is used in a ∼ 10 PW γ-ray flash experiment, κγ will be significantly suppressed.The main drawback of lithium is its high chemical reactivity in air.Fortunately, a thin (sub-micron) polymer coating does not significantly alter the laser interaction with the target since it is drilled by the background field before the main pulse arrives.In addition, the laser-target interaction takes place in vacuum.

Gamma-Ray Flash Parameters
So far, we have quantitatively described the laser generated γ-ray flash at various laser powers and for various target materials, demonstrating that κγ can reach significantly high values if the laser background field is not suppressed.Here, the absolute value of κγ is obtained by integrating γ-photons emitted in a 4π solid angle.It has been previously demonstrated for a linearly polarized laser that either by employing specific targets (cylindrical channel (70,71) or tightly focused lasers (72), γ-photon emission is stronger in two directions along the polarization plane, forming a double-lobe pattern.We generalize these results, showing that the double-lobe pattern is a general characteristic of linearly polarized laser-matter interactions, as shown in Fig. 6(A).The figure shows the γ-photon radiant intensity, I Ω (radiant energy per unit time per solid angle), by approximating that the γ-photon emission duration equals to the laser pulse duration.

D R A F T
The double-lobe feature is observed either for tailored or for foil targets, and at all laser power levels.The main difference among the various cases is the γ-ray flash amplitude, as seen in Fig. 6(B) showing a line-out of the γ-photon radiant intensity, and by considering photons emitted within 10 o full-angle with respect to the polariztion plane.The peak of Fig. 6(A) is slightly higher than the peak of Fig. 6(B) due to the 10 o averaging of the radiant intensity.
Mapping the γ-photon radiant intensity is crucial in experiments aiming on γ-photon production through laser-matter interactions.If one needs to optimize γ-photon detection, the detection system must be aligned along the highest radiant intensity direction, which in all cases lies on the laser polarization plane.However, the polar angle depends on the laser-matter interaction parameters, where the peak can be seen in Fig. 7(A) for tailored targets (orange line) and foil targets (green line).The ranges on the plot define where atleast 90 % of the peak signal is detectable.At all laser power levels, the tailored targets compared to foil targets result in γ-photon peak emission closer to the laser propagation axis, reaching ∼ 37 o versus ∼ 47 o respectively for a ∼ 10 PW laser.Furthermore, in both target cases the peak emission angle moves closer to the laser propagation axis for increased laser power.The γ-ray flash is usually treated as a time-integrated quantity of a single γ-photon pulse.However, ballistic evolution of the γ-photons allows an angular and temporal discrimination, as shown in Fig. 7(B); the temporal axis has an arbitrary (detection distance related) offset, not relevant to revealing D R A F T the γ-flash dynamics.As seen in the figure, γ-photon emission is directly connected to the laser wavelength in the regions of strongest emission.Temporal discrimination of the γ-ray flash reveals that the localized radiant intensity can reach values as high as ∼ 6 PW sr −1 .In addition, the γ-photons emitted at symmetric azimuthal angles come with a time difference of half laser period.
The reduced γ-photon production on the laser propagation axis indicates that the γ-photons do not originate from electrons moving antiparallel to the laser propagation direction, where in that case a peak should had been observed rather than a dip.Instead, the γ-photons are mostly emitted from electrons moving at a non-zero angle along the laser propagation.The lower energy γ-photons that are emitted as an afterglow at later times have an antiparallel direction with respect to the laser propagation axis and are characterized by a larger period.If a circular detector array is centered at the laser focal spot, then the γ-ray flashes will be detected with a time delay, as indicated by the non-centered γ-photon fronts of Fig. 7(B), indicating that the major γ-photon emission happens at a shifted position, at a certain depth inside the target cavity.

Conclusions
We present the efficient generation of γ-ray flash through simulations of an ultrashort ultraintense laser interacting with a solid target, where both the background field and main pulse are considered.The effect of the background field is taken from previously published MHD simulation results.The MHD simulations exhibit a unique favorable tailoring of the initially flat target, drilling a conical-like cavity in its volume preceded by an exponentially increasing preplasma distribution.The resulting electron and ion number densities are used as initial conditions for PIC simulations, revealing the effect of the main laser pulse on the tailored target.
The conical-like cavity formation strongly alters the laser interaction compared to the interaction with a foil target.The laser pulse, being of similar temporal extent to the cavity depth is reflected by the cavity walls, increasing the laser field intensity by an order of magnitude.Furthermore, cavity reflection of the laser field allows appearance of a longitudinal electric field component, aiding in further propagation of the cavity.In addition, if a relatively low density solid is chosen as the target material, then the current multi-PW lasers are relativistically transparent to the target, penetrating the target volume more deeply.
The ultrahigh intensities reached are within the regime where nonlinear Compton scattering process dominates the laser-target interaction, being capable of generating a bright γ-ray flash.Although no e − -e + pairs were expected from currently available lasers, we estimate that a significant number of pairs can still be produced due to the multiphoton Breit-Wheeler process by the significant intensity enhancement.It is found that the emitted γ-photons can have an energy approximately equal to the one third of the initial laser energy, with γ-photon energies approaching the GeV-level, if a ∼ 10 µm thick lithium foil is used as a target and by introducing an appropriate background field.
The target material is of crucial importance, since a denser material shifts the interaction out of the relativistic transparency regime.In addition, a denser material results in a shallower cavity, that can not effectively intensify the main laser pulse.As a result, both dense and/or foil targets suppresses the γ-ray yield since the laser field is strongly reflected, which can also pose a higher risk for optical damage of the laser system.
It is found that the γ-photon radiant intensity forms a double-lobe pattern emitted in the polarization plane, which is a general feature of linearly polarized lasers interacting with matter.The γ-photon emission peaks at angles of 37 o − 55 o under the parameters examined.In all cases, the emission angle is smaller when the cavity is present in the target.
The temporal dynamics of the γ-ray flash is also revealed, exhibiting direct connection to the laser period.It is found that for azimuthally symmetrical angles the γ-photon emission is temporally shifted by half laser period and is suppressed on the laser propagation axis.This pattern indicates that γ-photon emission originates from electrons co-moving with the laser pulse at a certain angle.Those electrons moving exactly parallel to the pulse produce little-to-no γ-photons, where this behavior explains the double-lobe pattern observed.In addition, the temporal discrimination of the γ-ray flash reveals radiant intensities significantly higher than previously expected, opening the road to their application in studying astrophysical processes and to effects of photonuclear interations.

Fig. 1 .
Fig. 1.Electron number density as given by MHD simulations (data taken from reference (57)), following background laser field irradiation of a lithium foil.The yellow contour line is at the critical density and the white contour line is at the lithium solid electron density.The orange saturated contour is overcritical for laser intensities above 10 27 Wm −2 .
) (to more accurately resolve electron trajectories), Bremsstrahlung and Photons (60) prepocessor directives enabled.The simulations are performed in the 3D version of the code, where the laser pulse and target characteristics match D R A F T those of the MHD simulations, with main pulse focused intensity ranging from ∼2.8 × 10 26 W m −2 to ∼2.8 × 10 27 W m −2 .The corresponding energy is set from 20 J to 200 J, while a pulse duration of ∼ 17 fs (4, 7) corresponds to lasers of ∼ 1 PW to ∼ 10 PW power.The central laser wavelength is set to 815 nm, typical for titanium-sapphire lasers.The laser focal spot is set at the center of the front surface of a foil target.The MHD density array, after interpolation to a 10 × 40 × 40 nm cell size grid, is centered and expanded in a cubic PIC volume, extending from −15.36 µm to 15.36 µm in all three directions.The cell dimensions are chosen small enough to accurately resolve the relativistic skin depth, while 8 macro-electrons and 8 macro-ions are assigned to each cell.The simulation ran for 110 fs, enough for κγ to saturate and at the same time neither field energy nor energetic electrons escape the simulation box from the open boundaries.The laser is focusing at a simulation time of ∼ 65 fs, which is set as 0 fs for the laser-target interaction.

Fig. 2 .
Fig. 2. (A) Laser intensity (color surface plot) overlayed to the lithium electron number density (grayscale image) at ∼15 fs, when an intensity of ∼2.6×10 28 Wm −2 is reached.(B) Overlay of three successive electron number density distributions, with a time step of 20 fs.The first layer is at ∼ − 5 fs, when the main pulse is within the cavity.The overlay of the three layers reveals the temporal dynamics of the cavity formation within the target.

Fig. 3 .
Fig. 3. (A) Line-out of Ey, and (B) line-out of Ex along the laser propagation axis at various times, as noted on the legend.Note that the laser is linearly polarized, along ŷ. (C) The left axis shows the maximum field intensity as a function of simulation time; the black dashed baseline denotes the field intensity expected in the focal spot, without the target.The right axis shows with the dashed purple line the amount of laser energy transferred to all target particles, while the solid purple line shows the percentage of the remaining laser field energy, κγ , as a function of time.

Fig. 4 .
Fig. 4. (A) Energy spectra after a ∼10 PW laser interacts with a tailored lithium target.(B) Energy spectra of γ-photons after interaction of the lithium target with main laser pulses of varying power, as noted on the legend.The solid lines correspond to tailored targets, and the dashed lines to foil targets.

Fig. 5 .
Fig. 5. (A) The left axis shows κγ as a function of a0, for tailored (blue) and foil (green) lithium targets.The right axis shows the ratio of the two aforementioned cases.(B) κγ as a function of laser power, for various tailored target materials.The κγ value is proportional to the circle area.

Fig. 6 .
Fig. 6. (A) Radiant intensity of γ-photons due to the interaction of a ∼10 PW laser with a tailored lithium target.(B) Line-out of the radiant intensity of γ-photons, taken along the polarization plane and within a full-angle divergence of 10 o .The solid lines correspond to tailored lithium targets and the dashed lines to foil targets, for various laser powers as noted in the legend.

Fig. 7 .
Fig. 7. (A) Angle of the peak γ-photon radiant intensity as a function of a0.The tailored and foil target cases are shown by the orange and green lines respectively.The ranges for both cases show the region where the radiant intensity is above 90 % of its peak value.(B) Radiant intensity of γ-photons measured within a full-angle divergence of 10 o with respect to the polarization plane, as a function of time-of-flight.In the figure, for computational reasons, only γ-photons with energy larger than 5 MeV are considered.