Ultra-bright γ-ray emission and dense positron production from two laser-driven colliding foils.

Matter can be transferred into energy and the opposite transformation is also possible by use of high-power lasers. A laser pulse in plasma can convert its energy into γ-rays and then e - e + pairs via the multi-photon Breit-Wheeler process. Production of dense positrons at GeV energies is very challenging since extremely high laser intensity ~1024 Wcm-2 is required. Here we propose an all-optical scheme for ultra-bright γ-ray emission and dense positron production with lasers at intensity of 1022-23 Wcm-2. By irradiating two colliding elliptically-polarized lasers onto two diamondlike carbon foils, electrons in the focal region of one foil are rapidly accelerated by the laser radiation pressure and interact with the other intense laser pulse which penetrates through the second foil due to relativistically induced foil transparency. This symmetric configuration enables efficient Compton back-scattering and results in ultra-bright γ-photon emission with brightness of ~1025 photons/s/mm2/mrad2/0.1%BW at 15 MeV and intensity of 5 × 1023 Wcm-2. Our first three-dimensional simulation with quantum-electrodynamics incorporated shows that a GeV positron beam with density of 2.5 × 1022 cm-3 and flux of 1.6 × 1010/shot is achieved. Collective effects of the pair plasma may be also triggered, offering a window on investigating laboratory astrophysics at PW laser facilities.

Introduction

The rapid development of laser technologies promises substantial growth of peak laser intensities. As of today, laser intensity of 1022 Wcm−2 has been demonstrated[1] and several laser facilities such ELI[2] and APOLLON[3], towards focused intensities above 1023 Wcm−2 are under construction. The electron dynamics in such intense laser fields approaches the radiation-reaction dominated quantum-electrodynamics (QED) regime[4-6]. Bright γ-rays emission, e − e + pairs production, QED-cascade as well as energetic particles acceleration are highly coupled, forming an e − e +-γ plasma[7,8]. This is a totally unexplored research area, which opens up new avenues in high energy density physics, particle and nuclear physics, and high energy astrophysics such as γ-ray bursts, pulsars, and active galactic nuclei in laboratories[9-12].

Breit-Wheeler (BW) pair production is a physical process in which an e − e + pair is created in the collision of two photons[13]. At extremely high laser intensities, an important mechanism for e − e + pair production is the multi-photon BW process[7], which enables the laser energy to convert into copious e − e + pairs via photon-photon annihilation (). Here, ω is the laser photon frequency. Generally, the low energy photons are laser photons or thermal emission[14], while high energy γ-photons result from accompanying radiation by relativistic electrons propagating in medium (bremsstrahlung[14-17]), colliding with laser pulses (Compton back-scattering[6,18-21]), or oscillating in intense electric fields (synchrotron/betatron radiation and skin-depth emission[22-27]). At laser intensities of ~1023–24 Wcm−2, several schemes are recently proposed to generate bright γ-rays and dense e − e + pairs in the laser-driven QED regime[6,7,20-33]. Broadly, these attempts can be classified into three main categories: a super-intense laser or lasers interaction with (i) a solid Al target[7,23,24,28], (ii) near-critical-density (NCD) plasmas[20,29], or (iii) an energetic electron beam (or single electron[31]) driven by laser wakefield[6,21,32] in gas plasmas or by a conventional accelerator[19,33]. Among them, significant attention has been paid to the first category for overdense e − e + pair production, which may trigger collective effects or ‘medium-like behavior’, e.g., Debye Shielding, required for modeling astrophysical phenomena in laboratories. However, these studies are limited to either one-dimensional (1D) theory or 2D particle-in-cell (PIC) simulations, and the inherent deficiency, i.e, high laser reflection by the immobile target, results in relatively low laser energy conversion. A full understanding of QED physics during the laser-foil interaction and improving the laser energy conversion in real 3D configuration become pivotal and pressing.

At laser intensities of ~10[22,23] Wcm−2, much less attention has been given to e − e + pair production due to the small flux of the high energy γ-photon emission. Here, we propose a novel scheme for ultra-bight γ-ray emission and dense positron production at achievable laser intensities in the foreseeable future. As shown in Fig. 1, two elliptically-polarized (EP) lasers are incident onto two diamondlike carbon (DLC) foils. The electrons in the center area of foils (i.e. the region of the laser focal spot) are rapidly accelerated to high velocities by the laser pressure and form overdense relativistic electron layers (RELs). Due to strong electron heating and foil expansion, relativistic transparency of both foils occurs under certain conditions, such that the laser penetrates through one foil and collide with the counter-propagating compressed layer of electrons accelerated from the other foil. This symmetric configuration enables efficient Compton back-scattering and increases equivalent quantum invariants, resulting in ultra-bright γ-ray emission with an unprecedented peak brightness and dense GeV positron beam production. The laser intensity required is within the capabilities of future multi-PW laser facilities, paving the way to potential applications in nuclear and particle physics for fundamental research, laboratory study of astrophysics, medical imaging and material science[4,5,9-11].

Figure 1

Schematic diagram of ultra-bright γ-ray emission and dense e − e + pair production by counter-propagating lasers irradiating two diamondlike carbon (DLC) foils.

Results

Radiation pressure acceleration of both foils by EP laser pulses

There has been particular interest in radiation pressure acceleration (RPA) of ultra-thin foils using circularly-polarized (CP) laser pulses, as they are capable of suppressing fast electron generation and resulting in quasi-monoenergetic ion beams[34-42]. In experiments, however, EP laser pulses are far more general because it is usually very difficult to get pure CP pulses at extreme intensities. When such an EP pulse irradiates a thin foil, both the characteristics of a CP and linearly polarized (LP) pulse are present. On the one hand, like a CP pulse, the laser radiation pressure pushes the thin foil forward as a whole; on the other hand, using EP laser pulses can suppress the transverse instabilities, which is beneficial to the foil acceleration[43]. In our scheme, two colliding EP laser pulses with a polarization ratio as predicted in ref.[43] are employed. Meanwhile, DLC foils are used instead of common metals like Au and Al employed previously[15-17,20,22-24]. It is shown that the DLC foils can offer exceptionally high tensile strength, hardness, and heat resistance, making them ideally suited materials for self-supporting targets in experiments[36,44]; unlike immobile metal targets, the foil dynamics in our proposed scheme is coupled with the γ-ray emission, enabling efficient laser energy conversion to electrons, and then γ-photons and positrons. By using such low-Z thin foils, e − e + pair production is dominated by the multi-photon BW process, while alternative one, i.e., Bethe-Heitler (BH) process[45,46] related with high-Z atoms and thick targets can be ignored (see Methods).

Due to the surge of abundant secondary particles during the laser-target interaction, previous studies were usually limited to 1D or 2D simulations as 3D simulations are extremely computational expensive. Here we perform the first 3D simulation using QED-PIC code EPOCH[47,48] (see Methods) to demonstrate ultra-bright γ-ray emission and dense positron production from laser-driven DLC foils. Fig. 2 presents the main simulation results. During the foil acceleration at the initial stage, carbon ions and protons are separated from each other with protons at the leading edge of the carbon ion front, due to their higher ion charge-to-mass ratio[39,49,50]. Note that , which is not strong enough to blow all electrons out of the foil. Thus a depletion region with a length of x for electrons forms in the foil. This generates strong charge separation fields, which in turn pull ions forward. At and (see Methods) on both sides, the electron density has a maximal value in the x-axis, which is high enough to prevent the laser pulse from penetrating through the foil. During this initial stage, both foils almost keep the structure intact and the degree of electron heating is not significant in such a short time duration. Due to the transversely Gaussian distribution of the laser intensity, both foils are distorted and a cone structure forms at 12T 0, as shown in Fig. 2(a) and (d). Some electrons within the skin-depth of the foil’s inner wall are pulled out by the traveling laser pulses and interact with the reflected waves[22]. However, both the electron energy and number are relatively small and the reflected waves are much weaker due to Doppler red shifting, so that only a few γ-photons are emitted and no positrons are produced at this stage (see Fig. 2(g)). This is markedly different from the cases in refs[22-24], where ~1024 Wcm−2 laser pulses are adopted and e − e + pairs can be produced during the hole-boring RPA.

Figure 2

Distributions of the transverse electric field (a–c), electron density (d–f), and photon energy density (g–i) at t = 12T 0, 16 T 0 and 20 T 0, respectively. Here, Vm−1. The white dashed circle in (c) refers to the interaction zone (IZ) and the red arrow in (e) points to the relativistic electron layers (RELs).

With the foil deformation and thermal expansion, the electron density along the laser axis decreases, while its energy increases. At the same time, Rayleigh-Taylor-like instability (RTI) grows quickly[51-54]. The physics associated with the foil acceleration is complex in full dimensions here, however, we may get deeper insights from the the 1D light-sail (LS) mode[37-40]. Taking into account only the center foil with a size of the laser focal spot where is of interest for the multi-photon BW process, the foil dynamics is well described by the motion equation , where is the laser intensity, , E is the laser electric field in the y-axis, is the reflectivity in the rest frame of the foil with , , β is the foil velocity, , and is the foil areal mass density with , n and Z being the ion density and mass number, respectively. Since we are interested in only the carbon ion front, we may assume perfect reflection of the incident laser pulse by the foil. Simplifying the motion equation above, we thus get

Introducing a new parameter and solving Eq. (1) analytically, we obtain the γ evolution of the foil in the 1D LS mode aswhere , and with (SI). Fig. 3(a–e) present the carbon ion energy distribution in space and Fig. 3(f) shows the evolution of the maximal carbon ion energy per nuclei. Here we only consider the ion acceleration in a smaller size of the same order of the laser spot. We see the center part of both foils already collide with each other at 16T 0, before which the simulation results agree well with the predictions of Eq. (2), demonstrating the dominance of RPA ion acceleration in the region of the laser focal spot.

Figure 3

Carbon ion energy distributions (a–e), energy evolution of the carbon ion front (f), electron energy spectrum (g), and evolution of and the maximum electron density n (h). The solid curves in (h) serve as guides to the eye.

For the laser intensity under consideration, the radiation friction force does not strongly affect the ion acceleration[55,56]. By use of DLC foils, the electron acceleration becomes inefficient because of the heavier carbon ions. The electrons average energy is lower, leading to a much smaller ration of the friction force to the Lorenz force, . Here, , r is the classical electron radius, and is the electron velocity. It is shown that the friction force only equals to the Lorenz force when the threshold a 0 ~ 350 is met (assuming in the center zone of the simulation box). Meanwhile, the radiation reaction for electron propagation in the direction of a propagating wave is much smaller than that produced due to counterpropagation[57]. This further mitigates the radiation reaction effect on the electron dynamics. The agreement of the results from the 1D RPA model with the simulation results in Fig. 3(f) justifies the accuracy of the 1D model without consideration of radiation reaction for the ion acceleration in this stage. Although the radiation reaction has limited impact on the ion acceleration, the electron dynamics is affected, which becomes more important at later stages and results in the foil expanding in both longitudinal and transverse directions.

Relativistically induced transparency of both foils

As the foil undergoes significant deformation, the laser pulses are further focused in the inner bent surfaces and the laser intensity along the axis increases, so that a large number of electrons escape the foils at t = 14T 0. The accompanying electron heating becomes important and the RTI develops quickly, leading to the fast expansion of foils along both the longitudinal and transverse directions. However, the later development of the RTI is favorable here, which enables the foil slightly transparent for the incident pulse at t = 13T 0 and completely transparent before t = 15T 0. Finally, most of the laser waves penetrate through the foil, propagate along the laser axis, and collide with the relativistic electrons from the other side, as is seen in Fig. 2(b) and (e).

Figure 3(g) shows the evolution of electron energy spectrum. At t = 14T 0, most high energy electrons with originate from the formed RELs as pointed out by the red arrow in Fig. 2(e). The cutoff energy is up to 5 GeV, but it decreases to 3 GeV at t = 18T 0. Meanwhile, the corresponding average electron energy also reduces from 92 MeV to 83 MeV. In order to explore the underlying physics, we first consider the time evolution of average electron energy () and the maximum density () in the center area as shown in Fig. 3(h). We see the electron density decreases with the time, gets closer and closer to the curve and finally intersects the curve at about 15T 0. This intersection point refers to the time instant when the laser pulse penetrates through the remainder of the foil in the center area, since the relativistically corrected plasma frequency becomes smaller than the laser frequency[58-60]. The electron temperature is related to the RTI and also the radiation reaction effect, whose increase may result in the foil expanding in both longitudinal and transverse directions. We thus consider two limit mechanisms, that of a transversely expanding thin foil and of a heated plasma that expands longitudinally. Firstly, we can interpret the transparency occurrence of a thin transversely expanding foil by Bulanov’s condition[58,59]. For a laser pulse with , a foil is transparent to relativistically strong laser radiation if is satisfied. Here and is the shell transverse expansion factor. We can rewrite the above condition as

From our simulations, it is shown that the surface stretching factor scales as in two laser periods, which is compatible with the traverse expansion velocities of the order of c/3. Thus we get a reduction of the opacity by a factor more than 2 between 13T 0 and 15T 0. It is sufficient to make the foil locally transparent. As time progresses, the whole foil becomes transparent to the incident laser pulse. On the other hand, considering the significant electron heating and foil longitudinally expanding in space, we may assume the foil to be a bulk plasma and hole-boring RPA may also occur at the later stage[61]. This can be roughly modeled as a semi infinite plasma, since the expanded foil’s width is larger than the electron skin depth during the acceleration. For an initial overdense plasma foil with a sharp boundary, the pondermotive force pushes electrons into the plasma, which creates strong peaking of the plasma electron density and results in significant enhancement of the laser threshold of penetration. Here, the laser threshold intensity a at the electron depletion boundary and on both sides satisfies the relation[62,63] . When the incident laser amplitude approaches its maximum value, the corresponding boundary amplitude a has exactly the same value at which the electron density vanishes. Considering here, we obtain the electron density threshold for the foil transparency by , which approximates as . In our scheme, due to the multi-dimensional effects such as the foil expanding, laser focusing, electron heating, and the RTI, partial foil transparency is first induced at locations with the electron density smaller than the threshold . This agrees well with simulation results in Fig. 2(e) and (f). In spite of two different models, both predict the occurrence of the foil transparency at about t = 15T 0.

Enhanced γ-photon emission and dense e−e+ pair production

In the QED regime, the probability rate for photon emission is characterized by the quantum invariant[4] , where is the normalized QED critical field[64], Vm−1, is the electron energy and is the electron momentum, E and B are the electromagnetic field experienced by electrons. For an electron co-propagating with the laser pulse, this results in and there are almost no high energy γ-photons emitted. On the contrary, when the electron counter-propagates with the pulse, one gets and efficient γ-photon emission is achieved if [28]. As the high energy γ-photons collide with the low energy ones, it has a probability for e − e + pair production via the multi-photon BW process, which is determined by another quantum parameter[4], . Here , , and ω the photon frequency. One sees that the pair production severely hinges on the local photon energy and the electric field E in the interaction zone. As χ approximates to the unity, QED-cascade induced avalanche-like pair production occurs[31].

Figure 4 illustrates the numerically calculated two key quantum parameters. Here the maximal value of at t = 16T 0 and it gets much smaller on the lateral side of the deformed foils. The large results from the symmetric configuration in our scheme, which increases the equivalent quantum invariants of both and χ . Finally, QED effects are triggered and a large number of γ-photons are emitted with the maximal photon density up to 250n at t = 18T 0 as seen in Fig. 5(a). Figure 5(b) presents the energy spectrum evolution of the scattered photons. The spectrum has a wide distribution ranging from MeV to several GeV. The average can be approximated theoretically by , where is the amplitude of mean electric fields. Thus the characteristic photon energy is written as[47] MeV at 18T 0, in excellent agreement with the simulation results in Fig. 5(c). Figure 5(d) presents the laser energy conversion efficiency to the electrons and secondary particles. We see the electron’s energy-share increases almost linearly until 13T 0 as the foil transparency occurs, and then starts to decrease. This can be attributed to the abundant γ-photon emission on both sides. As time progresses, the photon number increases but the average photon energy remains after t = 18T 0, since the high energy photons are continually consumed by the multi-photon BW process. The total laser-to-photon energy conversion efficiency is about 6%, which is comparable with that reported in a single foil[22], though our laser intensity is only 1/4 of that. Note that this is a full 3D QED-PIC simulation and only γ-photons with energy larger than 1 MeV are counted.

Figure 4

Numerically calculated quantum parameters (a) and (c) at t = 16T 0. (b) and (d) present the corresponding transverse distributions of both parameters at x = 10 μm.

Figure 5

Density distribution of γ-photons at t = 18T 0 (a) and energy spectrum of γ-photons (b). Evolution of the total number of γ-photons and their mean energy (c). Laser energy conversion efficiency to electrons, positrons, and γ-photons (d). In (c) and (d) the points refer to the simulation results and the curves serve as guides to the eye.

The significant increase of the high energy γ-photons by comparing Fig. 2(g–h) can be attributed to Compton back-scattering of the transmitted laser photons from one side off the relativistic electrons from the other side. In this regime, nonlinear interactions in which an electron scatters simultaneously with two or more laser photons occurs. This is a highly-nonlinear process and prevails at ultra-high laser intensities. For a head-on collision, the maximum energy of scattered photons in a single-photon process[20] is peaked at , with average photon energy of . Apparently, it is insufficient for an even 5 GeV electron to back-scatter the laser photons to such high photon energies in a single-photon process. Meanwhile, multiple inverse Compton scattering of scattered photons colliding with the energetic electrons is in principle possible, but the Compton cross section decreases by a factor 106  [33], so that it fails to account for the huge number of γ-photon emission on both sides. This further proves that the high-order nonlinear Compton scattering dominates the photon emission in our scheme. Overall, the total γ photon number above 1 MeV is 7.23 × 1013 at t = 20T 0 in such a μm × μm size with >6% laser energy conversion to the γ-photons. Conservatively, we take 2 × 1013 photons, a 1/e 2 source radius of 2λ 0 × 2λ 0, a divergence of 0.1 × 0.1 rad2, a bandwidth of 100%, and a pulse duration of about 5T 0. This gives an unprecedented peak brightness of about 1025 photons/s/mm2/mrad2/0.1%BW at 15 MeV, which is four orders of magnitude higher than the current experiment record in laser-based Compton sources[18] and 20-fold increase over that in the laser wakefield scheme with similar laser parameters[32]. The corresponding γ-rays intensity is 5 × 1023 Wcm−2 and the total emission power is 33 PW, almost 103 times larger than reported in ref.[65].

The numerically calculated parameter χ determining the multi-photon BW process is given in Fig. 4(b). Using the above , we rewrite the average value of χ as . Although it is far from the threshold for the QED-cascade, it is sufficient to trigger the multi-photon BW process. Figure 6(a) and (b) show the positron density distribution at t = 16T 0 and 20T 0. The mean positron density achieved is about 4n with the maximal density up to 25 n along the laser axis at 20T 0. This may provide us with an unique opportunity to investigate astrophysical phenomena in laboratories, which otherwise is impossible in the recent BH experiments[15,46]. Figure 6(c) presents the evolution of the positron energy spectrum. The positrons are continually accelerated by the laser pulse and their numbers increase significantly. However, these positrons also oscillate like electrons in the laser fields and emit γ-photons so that the number of high energy positrons decrease from 18T 0 onwards. Meanwhile, the beam divergence decreases from 18T 0 (see Fig. 6(d)) and most of the positrons are directed along the laser axis, indicating a good prospect for potential applications in future.

Figure 6

Distribution of the positron density at t = 16T 0 (a) and 20T 0 (b). The density above 1n is marked as red and the maximum density is up to 25n . Energy spectrum of the positrons (c) and their divergence angle distribution at given times (d).

Our full 3D simulations show that the final positron number is 1.6 × 1010 at 20T 0, which is much larger than that in the laser wakefield scheme[6,21,32]. The average positron energy is as high as 230 MeV. Assuming the volume of the focal area where e − e + pairs are created to be ~(λ/2)3, we may estimate the maximal positron flux by[8,29].where is the fine-structure constant. Taking the peak electric field E ~ 285E 0, we get , which is on the same level as the simulation results above. The corresponding laser energy conversion efficiency from the lasers to positrons is about 0.02%, as shown in Fig. 5(d). Taking and , the relativistically corrected skin depth of the produced e − e + pair plasma is δ  ~ 2λ 0. Thus the ratio of the plasma transverse size to the skin depth, , so that the collective effects of the pair plasma may be also triggered, e.g., Debye Shielding and plasma oscillation, which in turn affect the electron and positron dynamics. Due to the small resolution of the pair plasma (), it is impossible to resolve the collective dynamics of the pair plasma. This is beyond the scope of this paper, but it is of significance for investigation of laboratory astrophysical phenomena in laboratories[9-11].

The positrons obtained are well confined in a small size of the same order of the laser focal spot, which can be clearly seen from their phase space distributions exhibited in Fig. 7. This is credited to the standing waves formed directly by the colliding laser pulses, which help to trap the pairs at the nodes (E = 0) such as x/λ 0 = 9.2,9.7,10.0,10.7, and prevent them from dispersing in space, so that the positrons can remain high energy and density until the laser pulses pass away. The wave structure with half laser wavelength spacing in Fig. 7(a) and (c) show the clear evidence of the positrons dynamics modulated by the standing waves.

Figure 7

Phase space distribution (x,p ) of positrons at t = 16T 0 (a) and 18T 0 (c). The positron density distribution along the laser axis and the corresponding electric field E at t = 16T 0 (b) and 18T 0 (d).

Discussion

There are several key factors to be considered for the experimental test and development of this concept. Firstly, the gap between the two foils plays an important role, because it determines the electron energy, density and flux. Since the carbon ion front moves with the proton layer, we can estimate the foil transparency time τ implicitly as

In our simulations, the foils become completely transparent at about 15T 0. Thus the critical foil gap is G > 2 × 6.94 − 2 × 0.32 = 12.88λ 0 for the laser and foil parameters under discussion. This is in excellent agreement with the simulation results. Choosing a proper foil gap thus becomes a key issue in experiments. For example, a 10 μm gap is insufficient because the electrons have not yet been accelerated to high enough energies before the collision. In contrast, a 16 μm gap is so large that the transverse RTI develops significantly and the foil collapses more quickly before they collide. The resulting electron density and total number become lower, which is definitely undesirable for the γ-photon emission. Meanwhile, compared with a single foil[22-24], which lacks flying foil’s dynamics at the initial stage and thus has smaller electron energies, we achieved more photons and positrons production due to the enhanced quantum parameters and χ . However, the spacing between the foils will also depend on the temporal profile of the laser pulses and an ultra-high laser contrast is required to facilitate the efficient RPA of DLC foils.

Secondly, we also find that the foil thickness and laser intensity play an important role in the foil dynamics and QED process. The laser intensity should be such that the foil can undergo a sufficient degree of acceleration and the transmitted light is still intense enough to penetrate the foil at the later stage. Considering the fact that κ depends on the foil thickness and laser intensity significantly, we can reduce the laser intensity and foil thickness simultaneously to keep κ unchanged. This leads to almost the same foil acceleration as above. For example, taking two 8 × 1022 Wcm−2 EP lasers with the same and DLC foils with thickness of , similar γ-photon emission and positron production can be achieved. The resulting brightness of the γ-rays and the maximum density of positrons are on the same order of magnitude. We have carried out a serials of 2D simulations to identify the optimal condition. It is finally shown that the optimal laser intensity should obey to enable highly efficient γ-photon emission and dense positron production.

Thirdly, in order to check the influence of the laser temporal profile on the γ-photon emission and pair production, we considered a Gaussian time profile. In this case, the full-width-at-half-maximum (FWHM) of the laser pulse is τ = 6T 0, to keep the total laser energy unchanged. The details of the laser parameters and simulation results are shown in Table 1 and Fig. 8. One can see that the magnitudes of the average energy, total number and energy spectrum of γ-photons in both cases are comparable. This indicates a weak dependence of γ-photon emission on the laser temporal profile in our scheme. Since the key parameter χ is determined by the average laser intensity, the positron production will be more affected by use of a Gaussian time profile. The final positron number obtained is up to 0.6 × 1010, which is about 40% of that in the case with a trapezoidal profile. However, the average energies of positrons obtained are almost the same in both cases.

Table 1

Comparison of γ-photon emission and positron production in two cases with different laser temporal profiles at 20T 0.

	Parameters	Trapezoidal profile case	Gaussian profile case
Laser	Time duration, τ 0	(1-8-1)T 0	6 T 0
	Focal spot, σ 0	4λ 0	4λ 0
	Peak intensity, I 0	1.1 × 1023 Wcm−2	1.1 × 1023 Wcm−2
	Total energy, ε L	1580 J	1538.6 J
γ-photons	Average energy, ε_γ¯	14.6 MeV	13.3 MeV
	Total number, Nγ	7.23 × 1013	4.80 × 1013
Positrons	Average energy, ε_e^+¯	230 MeV	225 MeV
	Total number, N_e^+	1.61 × 1010	0.6 × 1010

Figure 8

Energy spectrum of γ-photons and positrons in both cases at t = 20T 0.

In our configuration, bremsstrahlung is ignored since its cross section is highly dependent on the target atomic number[66,67] and is proportional to . Thus the cross section is significantly reduced by use of a low-Z target in our scheme. Additional simulations with the bremsstrahlung considered show that both the photon energy and numbers resulting from the bremsstrahlung is 105 times smaller than those from the nonlinear Compton scattering. In our simulations, the trident process is switched off and the BH process is also ignored because the yield of the BH is expected to scale as Z 4  [46]. It is interesting to compare our results with those observed in recent BH experiments. We see the obtained positron density is six orders of magnitude higher than reported in BH experiments (~1016 cm−3)[15,46] and 80 times higher than predicted in the laser wakefield scheme[32]. The laser-to-positron energy conversion efficiency is at the same level as that in the BH experiments[46] and much larger than in the laser wakefild scheme[32].

Since the laser-driven DLC for ion acceleration has been demonstrated in recent experiments[36]. So the main issue for our scheme in future experiments would be the co-aligning of the two laser beams as well as their synchronization. Techniques have been developed for the co-aligning of two laser beams[68]. In addition, there is active work to synchronize high power laser beams to better than 10 fs (approaching femtosecond level). The two beams of 0.5 PW of Astra-Gemini laser can be synchronized to better than ±50 fs at present and this value will be improved in the near future[69]. The femtosecond synchronization of the two fs laser pulses can be obtained because both pulses are split from one fs pulse (after the laser oscillator), travel nearly identical optical paths (in the laser amplifier chains) and the small temporal differences are compensated at the end, before recombination. In additional simulations, we assume the left pulse irradiates the target normally but the right pulse is incident with 2 μm off the laser axis of the left pulse. All other parameters are the same as above. Finally, we find that the photon number obtained decreases from 5.96 × 1013 to 5.85 × 1013 (i.e., 2%) and the positron number also decreases from 1.45 × 1010 to 1.27 × 1010 (i.e., 14%) at t = 18T 0. We therefore see a small reduction in the signal if the alignment is off by about the size of a laser focal radius. This further demonstrates the robustness of the scheme.

In conclusion, we propose an all-optical scheme for ultra-bright γ-ray emission and dense e − e + pair production by irradiating 1022–23 Wcm−2 colliding lasers onto two DLC foils which is investigated using full 3D PIC simulations. The electrons in the focal region of one foil are rapidly accelerated by the laser radiation pressure, forming flying RELs, and collide with the other laser pulse which penetrates through the second foil from the other side due to relativistically induced transparency. By incorporating the moving foil’s dynamics in the scheme, this dual interaction arrangement enables efficient Compton back-scattering and results in ultra-bright γ-photon emission with an unprecedented brightness of ~1025 photons/s/mm2/mrad2/0.1%BW at 15 MeV and intensity of 5 × 1023 Wcm−2. Finally, a GeV positron beam with a maximum density of 2.5 × 1022 cm−3 and flux of 1.6 × 1010/shot is achieved and collective effects of the pair plasma may be triggered, offering a window on investigating laboratory astrophysics at PW laser facilities.

Methods

Numerical modeling

The open source QED-PIC code EPOCH is employed, which has incorporated the binary collision, γ-ray emission, radiation reaction, and pair production by both the trident process and multi-photon BW process. A Monte Carlo algorithm with quantum correction is implemented in the code for calculating the photon emission and pair production by assigning each electron, positron, and γ-photon an optical depth, at which the emission occurs when the optical depth reaches a derived value from a pseudo-random number between 0 and 1. Bremsstrahlung is not included in EPOCH and is thus ignored in our scheme. Meanwhile, the trident process is switched off and the BH is also ignored because of the use of low-Z ultra-thin DLC targets. More details about the physical modules and algorithms can be found in the refs[47,48].

The simulation box size is X × Y × Z = 20λ 0 × 14λ 0 × 14λ 0, sampled by 2000 × 560 × 560 cells with 8 particles in each cell. Here λ 0 = cT 0 = 1 μm is the laser wavelength and T 0 = 3.3 fs is the laser cycle. The foil electron density is  = 200 n , mixed with 20% protons in number density. Both foils have a thickness of L = 0.32 μm and radius of R = 5 μm, and are placed symmetrically at a distance of D = 3 μm away from the left and right boundaries, as shown in Fig. 1. The foil gap is G = 13.36 μm. Two identical EP Gaussian laser pulses with intensity component of a  = 237 and a  = 154 and focal size of σ 0 = 4 μm irradiate the two foils simultaneously. Here, the ellipiticity of the laser pulse is chosen according to the simulation result[43], which is favorable for the ion acceleration by suppressing the transverse instabilities. For comparison to the 1D theory as discussed in ref.[8], each pulse has a trapezoidal profile in time with a duration of 10T 0 (1T 0 − 8T 0 − 1T 0). Thus the total laser energy of each laser pulse is about 1600 J. Here a typical 3D simulation presented requires two-week runtime using 600 CPU-cores on TianHe High Performance Computers at NUDT.

Data availability

Data associated with research published in this paper can be accessed at http://dx.doi.org/10.15129/5a84ecba-17bd-4750-971c-82593cf45d77.