High-energy-density electron-positron pair plasma production and its dynamics in the relativistic transparency regime
W. Y. Liu, W. Luo, T. Yuan, J. Y. Yu, M. Chen, and Z. M. Sheng

TL;DR
This study uses multi-dimensional particle-in-cell simulations to explore how laser-driven electron-positron pair plasmas are produced and behave in the relativistic transparency regime, revealing optimal conditions for maximum pair production.
Contribution
It demonstrates that target transparency significantly enhances pair production via stable standing waves, identifying optimal foil densities and laser intensities for efficient plasma generation.
Findings
Optimal foil density for pair production is 200-280 n_c.
Laser energy conversion to pairs is four times higher at optimal density.
Stable standing waves trap and modulate created pairs.
Abstract
High-energy-density electron-positron pair plasma production and its dynamics in a thin foil illuminated by two counter-propagating laser pulses are investigated through multi-dimensional particle-in-cell simulations. We compare the production of electron-positron pairs and gamma-photons via quantum electrodynamics processes in the relativistic transparent and opaque regimes, and find that the target transparency can significantly enhance the electron-positron pair production due to the formation of stable standing wave (SW). An optimum foil density of 200 - 280 n_c (n_c is the laser critical density) is found for enhancing electron-positron pair production when laser intensity reaches a few 10e23 W/cm2. At such foil density, laser energy conversion to electron-positron pairs is approximately four times higher than at foil density of 710n_c, whereas laser energy conversion toโฆ
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsLaser-Plasma Interactions and Diagnostics ยท Laser-induced spectroscopy and plasma ยท Atomic and Molecular Physics
High-energy-density pair plasma production and its dynamics in the
relativistic transparency regime
W. Y. Liu
School of Nuclear Science and Technology, University of South China, Hengyang 421001, China
Key Laboratory for Laser Plasmas (Ministry of Education),School of Physics and Astronomy,
Shanghai Jiao Tong University, Shanghai 200240, China
IFSA Collaborative Innovation Center, Shanghai Jiao Tong University, Shanghai 200240, China
โโ
W. Luo
School of Nuclear Science and Technology, University of South China, Hengyang 421001, China
SUPA, Department of Physics, University of Strathclyde, Glasgow G40NG, United Kingdom
โโ
T. Yuan
Key Laboratory for Laser Plasmas (Ministry of Education),School of Physics and Astronomy,
Shanghai Jiao Tong University, Shanghai 200240, China
IFSA Collaborative Innovation Center, Shanghai Jiao Tong University, Shanghai 200240, China
โโ
J. Y. Yu
Key Laboratory for Laser Plasmas (Ministry of Education),School of Physics and Astronomy,
Shanghai Jiao Tong University, Shanghai 200240, China
IFSA Collaborative Innovation Center, Shanghai Jiao Tong University, Shanghai 200240, China
โโ
M. Chen
Key Laboratory for Laser Plasmas (Ministry of Education),School of Physics and Astronomy,
Shanghai Jiao Tong University, Shanghai 200240, China
IFSA Collaborative Innovation Center, Shanghai Jiao Tong University, Shanghai 200240, China
โโ
Z. M. Sheng
Key Laboratory for Laser Plasmas (Ministry of Education),School of Physics and Astronomy,
Shanghai Jiao Tong University, Shanghai 200240, China
IFSA Collaborative Innovation Center, Shanghai Jiao Tong University, Shanghai 200240, China
SUPA, Department of Physics, University of Strathclyde, Glasgow G40NG, United Kingdom
Abstract
High-energy-density electron-positron () pair plasma production and its dynamics in a thin foil illuminated by two counter-propagating laser pulses are investigated through multi-dimensional particle-in-cell simulations. We compare the production of pairs and -photons via quantum electrodynamics processes in the relativistic transparent and opaque regimes, and find that the target transparency can significantly enhance the pair production due to the formation of stable standing wave (SW). An optimum foil density of 200ย โย 280 ( is the laser critical density) is found for enhancing pair production when laser intensity reaches a few . At such foil density, laser energy conversion to pairs is approximately four times higher than at foil density of 710, whereas laser energy conversion to -photons keeps almost the same. Consequently high dense plasma with a maximum intensity above is produced. Modulation dynamics of created pair plasmas is further observed when target foil becomes transparent. It is shown that stable SWs formed directly by two counter-propagating lasers, not only trap the created pairs to their nodes, but also modulate periodically average energy and phase-space and angular distributions of trapped particles. However, similar trapping and modulation effects become obscure in the opaque regime due to the absence of stable SW field.
I Introduction
Relativistic electron-positron () pair plasma represents a unique state of matter, which is believed to exist in many extreme astrophysical environments, such as Active Galactic Nuclei (AGN), Pulsar Wind Nebulae (PWN), Gamma Ray Bursts (GRBs), and Black Holes (BHs) E. Weibel1959 ; B. Paczynski1986 ; M. J. Rees1992 ; E. Waxman1995 ; T. Piran1999 ; T. Piran2005 ; P. Mรฉszรกros2006 ; Hui Chen pop2015 . In particular, high dense pair plasmas are generally accepted to be invoked to explain energetic phenomena associated with quasars, GRBs and BHs J. Wardle1998 ; I. F. Mirabel1999 ; P. Meszaros2002 ; G. Weidenspointner2008 ; G. Sarri2015 . On earth, relativistic positron sources are of paramount importance in experimental physics due to their potential applications to a wide range of physical subjects, including nuclear physics, particle physics, and laboratory astrophysics wluo2013 ; nd2016 . Generation of such pair plasma is attracting an increasing amount of worldwide attention H. Chen2011 ; Hui Chen prl2015 ; Jian-Xun Liu2015 ; Jin-Jin Liu2016 ; Tongjun Xu2016 ; yuantao2017 , aiming to replicate the physics of small-scale astrophysical environments in the laboratory and to elaborate potentially board applications.
With everlasting and progressive development of laser technology, laser beams with intensities of are already realized V. Yanovsky2008 . The underway project Extreme Light Infrastructure (ELI) eli and proposed initiatives like XCELS xcels and iCAN G. Mourou2013 are expected to achieve intensities above within the next few years. At such intensities, the normalized amplitude of laser vector potential , where is the angular frequency of the laser, and laser-solid interaction enters the ultrarelativistic regime L.Willingale2009 , in which quantum electrodynamics (QED) processes may arise G. A. Mourou2006 ; E. N. Nerush2007 ; A. R. Bell2008 ; A. M. Fedotov2010 . One of the most representative QED processes is nonlinear Compton Scattering, A. DiPiazza2010 ; F. Mackenroth2011 ; C. P. Ridgers2012 . In this nonlinear process, foil electron accelerated (in the laser focus) spontaneously absorbs multiple laser photons and then radiate a high-energy photon . The radiated photon in strong electromagnetic (EM) fields can further decay into an pair via multi-photon Breit-Wheeler (BW) process, G. Breit1934 ; T. Erber1966 ; J. G. Kirk2009 . As QED effects become significant, prolific charged particles produced in the system can focus large current into the classical electron-ion plasma, which strongly modify the plasma dynamics. The EM fields, which are determined by the plasma dynamics, in turn affect the rates of the QED processes. As a result, the physics associated with the formation and dynamics of dense pair plasmas produced by QED-strong laser interacting with solid foil is of paramount importance that may prove essential for developing an astrophysical relevant research platform.
In this paper, we present enhanced pair plasma production and its dynamics in thin foils undergoing relativistic transparency. The case of a thin foil irradiated by two counter-propagating lasers (two-side irradiation scheme) is particularly interesting, because in contrast to single laser solid interaction, two-side irradiation scheme can significantly enhance the QED effects Wen Luo2015 ; H. X. Chang2015 . Here we revisit this irradiation scheme in the relativistic transparency regime and emphasize on how this affects the pair production and the following particle dynamics. When target foil undergoes transparency, stable standing-wave (SW) field can be produced directly by the overlap of two counter-propagating laser pulses. Such SW fields significantly enhance the pair production by the BW process and then produce high dense pair plasmas. Modulation dynamics of the positron energy and phase-space and angular distributions is further observed when transparency occurs.
The reminder of this paper consists of five sections. Section II describes the simulation setup. Section III demonstrates dense pair plasma and -ray burst production in both the relativistic transparency and opaque regimes. In Section IV, the modulation of stable SW fields on produced pair plasmas in the relativistic transparency regime is discussed. Section V presents a brief discussion and summary on this work.
II simulation setup
An open-source particle-in-cell (PIC) code EPOCH C. P. Ridgers2014 was used to perform multi-dimensional simulations. This code has been incorporated with a radiation-reaction module and a collective QED module, allowing self-consistent modeling of laser-plasmas interactions in the near-QED regime. Two linearly polarized (LP) laser pulses with an equal wavelength of and duration of focus on, respectively, the left and right boundaries of the simulation box at time . Here is the period of laser pulse. The lasers have transversely super-Gaussian spatial profiles with electric field as and focus to spots with radius of . The size of simulation box is . A thin foil is fully ionized carbons and hydrogen, which locates between to from the left boundary. Ratio of density of carbons and hydrogen is 1:1. The foil is discretized on a spatial grid with cell size of 10 nm and is represented by 500 macro electrons and 16 macro ions per cell. A set of initial foil densities () and laser intensities () were used in these simulations.
III enhanced \textnormal{\emph{missing}}{e^{-}e^{+}} pair production
When two counter-propagating laser pulses irradiate on a thin foil from opposite sides, laser hole-boring (HB) stage is initiated. At this stage, the foil electrons and ions are first accelerated in the laser field. Owing to the opposite and equal laser radiation pressures, these accelerated charged particles could hardly be escaped from the foil. The foil then becomes much denser with much higher reflectivity of the lasers, enabling formation of SW on each side of the foil. When HB stage ends, the foil starts thermal expansion stage, in which the sum of electrostatic and thermal pressures may exceed the laser radiation pressure on each side, leading to decompression and expansion of the foil H. X. Chang2015 . As the electron density of fully ionized foil becomes lower/higher than the relativistically corrected critical density, the foil plasma becomes transparent/opaque to the incident laser L.Willingale2009 . The results are shown in Figs.ย 1(a) and (e). It is seen that the foil plasma is underdense at such that the counter-propagating lasers can penetrate the thin foil [see Fig.ย 1(b)] and then form a stable SW. This SW has a period of 0.5. As foil plasma becomes denser, for example, at , it can be opaque to the incident laser [see Fig.ย 1(f)]. Then the incident and reflected lasers produce instable and relatively weak SWs at both sides of the irradiated foil, as shown in Fig.ย 1(e).
Figs.ย 1(c) and (g) shows the spatial distributions of hard photons with energies higher than 1.022 MeV. The number of radiated photons increases with the foil plasma density due to proportional increase of foil electron number. These photons have similar spatial patterns with those of foil electrons in the laser focus. At lower foil plasma density, we find that majority of photons are produced with relatively higher energies because the incident lasers can penetrate the thin foil and then enhances the amplitude of EM field. These photons in turn interact with the SW field sufficiently, producing more energetic pairs. The SW field can further bunch the produced pairs in space, leading to the formation of high dense plasma with a maximum intensity above , as shown in Fig.ย 1(d). Such bunching effect, which disappears in the opaque regime, is mainly caused by the SW formed directly by two incoming laser pulses. When the initial target density is 710, the interaction enters into the opaque regime. The majority of photons are produced at the high-density plasma region, where the incident lasers get reflected [see Fig.ย 1(f)]. In the opaque regime, energetic photons experience a relatively weak EM field when they escape from the interaction area, and then pair production via multi-photon BW process becomes deficient. Fig.ย 1(h) shows that the intensity of created positrons at is visibly lower than that at . Such intensity profile further displays a chaotic pattern due to the lack of stable SW.
In order to identify the relativistic transparency and opaque regimes, the characteristic behaviors of the right- and left-moving field components are shown in Fig.ย 2, which is obtained in simulations with laser amplitude . These components are defined by C. Baumann2016
[TABLE]
One can clearly see that in the case of , the onset of transparency begins at , which leads to a transient SW [see Figs.ย 2(a) and (c)]. The velocity of HB front is estimated as and the HB stage should be terminated at , which is in accordance with the simulation result. However, in the case of a visible interruption of laser propagation occurs until [see Figs.ย 2(b) and (d)], which indicates that the plasma remains opaque for this interaction stage. In the following stage, the injection of laser energy is almost finished and the reminder of the SWs formed directly by two incident lasers can hardly affect the production of pairs. This is the also the reason that the number of pairs obtained at higher foil density is visibly less than at lower foil density.
We investigate laser energy conversion to pairs () and -ray photons () in the relativistic transparency and opaque regimes. Fig.ย 3 shows such energy transfer as a function of foil plasma density at , which corresponds to the moment that the laser fields have vanished. It is seen that the increases first and then decreases rapidly with the foil density, whereas the increases first and then approaches a saturation value. For both laser amplitudes, and 660, the has a maximum value at foil plasma density of . Such density is slightly lower than the separatrix density, approximately given by , for which the relativistic transparency occurs at , that is, more than half of the laser pulses have reached the target. The separatrix density observed in the simulations is lower than the theoretical , due to target compression and QED plasma effects. The QED plasma effect can deplete the laser pulse energy, and thus reduces the laser amplitude in laser foil interactions. For laser amplitude , the target compression and laser depletion due to QED effect W. M.Wang increase significantly such that the separatrix density does not shift to higher value, compared to that for laser amplitude .
As displayed in Fig.ย 3, the obtained in the relativistic transparency regime is obviously higher than that in the opaque regime. For laser amplitude , when reducing foil plasma density from 710 to 280 the laser energy conversion increases from to , which is enhanced to about four times, so does the number of created pairs , from at to at . Meanwhile, the laser energy conversion keeps almost the same value, . The enhancement of pair production is mainly due to the target transparency. As discussed above, the transparency results in the formation of stable SW from a direct overlap of two counter-propagating laser pulses. Such pulse overlapping enhances the laser field strength, which can accelerate charged particles to higher Lorentz factors and then radiate more high-energy photons per each charged particle. In the following stage, the propagation of these energetic photons through the stable SW field increases the dynamical quantum parameter Wen Luo2015 ; H. X. Chang2015 , which controls pair production via multi-photon BW process. As a consequence, more QED pairs can be produced in foils undergoing relativistic transparency. We see in Fig.ย 3 that an optimum foil density region, namely, is found for realizing maximum conversion efficiency from laser photons to pairs.
IV modulation dynamics of \textnormal{\emph{missing}}{e^{-}e^{+}} pair plasma in the relativistic transparency regime
IV.1 Energy modulation
After transparency has occured, the stable SW can modulate the average pair energy. Fig.ย 4 shows temporal evolutions of ensemble-averaged energy for the created pairs at laser intensities of and , respectively. At the initial stage, the HB and thermal expansion effects are dominant over the formation of stable SW. The energy modulation is ambiguous. At , the plasma foil becomes transparent and stable SW fiels are formed consequently [see Figs.ย 2(a) and (c)]. During this stage, the average positron energy is modulated periodically [see Fig.ย 4]. The modulation period is 0.5, as a result of temporal evolution of the SW. On the contrary, the oscillation of average positron energy becomes insignificant when laser foil interaction occurs in the opaque regime, for example, in the case of for . This is due to merely the existence of instable SW fields formed by incident and reflected laser pulses. Furthermore, the foil does not undergo transparency at lower intensity of . Only a few pairs are produced and the modulation dynamics are invisible as well.
We present detailed modulation processes as follows. When SW field becomes strong at with being the integer, charged particles experience such field and radiate more energetic photons. As a result, average energy per particle decreases and these particles can be readily trapped to the electric nodes due to stronger radiative trapping (RT) L. L. Jiprl2014 and tight longitudinal confinement, as discussed later. On the contrary, when the overlap of the SWs diminishes at , charged particles radiate less; instead they can absorb some energy from the varying laser fields, which finally enhance their average energy. Such energy oscillation can maintain a few laser periods. After , it becomes weak and even disappears along with the fading of SW field. It is seen that the final average energy for BW-electrons approaches a median value, approximately 200 MeV. Meanwhile the positron average energy has considerable enhancement on acceleration of sheath potential generated by fast electrons as they leave the target. The existence of sheath potential also gives rise to the average energy difference (20 MeV) between the positron and BW-electronโs in the oscillation stage, as displayed in Figs.ย 4 (a) and (b).
The energy oscillation of pairs can be explained with conservation of canonical momentum of positrons (or BW-electrons) after they are spontaneously accelerated to relativistic velocities. The conservation of canonical momentum reads J. Meyer-Ter-Vehn192
[TABLE]
Here is the transverse momentum of positron, is the elementary charge, is the vector potential of EM field perpendicular to particle momentum, and the can be approximated as the median value mentioned above. During the interaction, the peak amplitudes of the SW field varied periodically, leading to that the trapped pairs have visible oscillation in the average transverse momentum. Meanwhile the longitudinal momentum keeps almost unchanged. The simulations have verified this behavior. A detailed discussion will be given in the following subsection.
At higher laser intensity of , regular modulation comes earlier than at lower laser intensity, as displayed in Fig.ย 4. The reason is that the relatively strong laser radiation pressure can shorten the stages of both laser HB and target thermal expansion, hence forming stable SW sooner. Furthermore, the overall charged particle energy decreases due to accumulation of considerable radiation loss, which cannot be compensated by the energy gain from the laser fields.
IV.2 Phase-space modulation
Besides particle energy oscillation, a visible phase-space modulation can be observed when thin foils undergo relativistic transparency. It is seen in Fig.ย 5 that the strength of SW field is always zero at nodes of , and maximum at antinodes of . The field strengths at antinodes alternate between strong and weak overlaps during each half of laser period. The created pairs are trapped to the nodes of the SW field when the field strength becomes strong at , which further helps forming stripe distributions along the direction of laser polarization [see an example shown in Fig.ย 5(a)]. The pair bunches have regular space interval of . After one-fourth of laser period, the bunched pattern evolves into a chaotic lobe when the field strength diminishes [see an example in Fig.ย 5(b)].
The positron density peaked at electric nodes can be clearly seen in Fig.ย 5(c), which plots the longitudinal distribution of positron number for a fixed moment in time . At this moment, charged particles experience an electric field almost doubles that they experience at [see Fig.ย 5]. Consequently these particles quiver violently in the intense field and loose substantial amount of energy by radiating high-energy photons. They in turn experience strong radiation-reaction (RR) force, leading to visible RT. A direct evidence for RT dynamics is the decrease of trapped particle transverse momentum L. L. Jiprl2014 , which is visibly shown in Fig.ย 6. Instead of being focused only in the electric nodes, more positrons tend to appear in the subspace between two electric nodes at [see Fig.ย 5(d)].
Figs.ย 6(a) and (b) show positron phase-space distributions at time . The distribution has transverse stripe structure periodically, which agrees with its spatial pattern shown in Fig.ย 5(a). Instead of being expelled transversely by the ponderomotive force, positrons are mainly trapped in the high-field region of SW. Consequently the distribution in Fig.ย 6(b) displays in a bright elliptical pattern in the center. However, the particle dynamics is different at time , which corresponds to that the field strength of SW at their antinodes. At this moment both the confinement of charged particles and RT effect loose. These positrons with tend to move forward, and those with move backward. As a result, the transverse stripe structure becomes oblique with respect to the longitudinal direction [see Fig.ย 6(c)]. The weakened SW field further disperses the positrons in phase-space of , since the RR effect becomes insignificant such that it is difficult to compensate their transverse dispersion.
IV.3 Angular distribution
We continue to investigate angular distributions of foil electrons, -photons and pairs. Since it is necessary to consider additional plasma effects, such as thermal expansion along -direction, more realistic three-dimensional (3D) simulations were performed. We define along the direction of the laser propagation, and along the direction of -polarization. Two-side irradiation scheme enables irradiating and compressing a thin foil more symmetrically, thus foil electron, -photon and positron have symmetric -intensity patterns with respect to the azimuthal plane [see Fig.ย 7(a)]. Since the laser intensity is much higher than the relativistic intensity of the order of , the effect of the magnetic and electric fields on the motion of charged particles should become comparable. Relativistic pairs produced in such field are predicted to quiver nonlinearly, moving in figure-of-eight trajectory, rather than in straight lines.
The modulation of polar distribution of positron intensity is seen from Fig.ย 7(a). In the presence of SW field, the polar patterns of positrons at time and 10.25 are โpinchedโ along the direction of laser polarization due to the confinement effect of the SW field. The polar pattern at time becomes more slim relative to that at . This is due to a tighter confinement effect longitudinally. In order to confirm such effect, we display in Fig.ย 7(b) that the โpinchedโ effect disappears when the SW field has already faded. It is also shown in Fig.ย 7(b) that different from BW-electron, the positron intensity have an isotropic distribution, due to the influence of sheath acceleration. For -polarized lasers, foil electrons experience transverse electric fields. They and hard photons they emit have tendency to move along the direction of laser polarization, ., .
Figs.ย 7(c) and (d) show azimuthal distributions of foil electron, -photon and positron intensities for and -polarized lasers, respectively. In both laser polarization cases, foil electrons have angular patterns with maximal intensity in their respective laser polarization directions. However, such pattern is a little โfatโ as a result of target thermal expansion after it symmetrically compressed by laser radiation pressures. Azimuthal distribution of -photon intensity is peaked on-axis of the laser polarization, suggesting high-order harmonic generation based on Compton scattering in highly nonlinear regime. The positron intensity displays similar patterns as the -photon. Also shown in Figs.ย 7(c) and (d) is temporal evolutions of SW field can hardly affect positron azimuthal pattern, which is different from modulation dynamics on positron polar structure.
V summary and conclusion
Thin foil undergoing relativistic transparency can enhance the pair production. The produced pair plasma has a density as high as . This implies the beam proper density of while considering the bulk Lorentz factor of . The beam of the produced pair plasma has transverse size of the order of . The relativistic corrected collision-less skin depth of the beam is thus , with being the relativistic plasma frequency. This value is smaller than the beam transverse size, indicating that collective (that is, plasma-like) behaviors are likely to occur in the beam. Since the created plasmas are dense enough to permit collective and kinetic behaviors to play roles, they are further similar to the condition of astrophysical events such as jets of long -ray bursts J. Wardle1998 ; I. F. Mirabel1999 .
In summary, the generation of overdense pair plasmas and following modulation dynamics have been investigated by simulations of two QED-strong lasers irradiating a thin foil from opposite sides. In the relativistic transparency regime the laser energy conversion to pairs can increase four times compared to that in the opaque regime, thus demonstrating enhanced pair production. At laser intensity of , the produced pair plasmas have a high energy density exceeding . Although the beam transverse size is only a few micros, the created pair plasmas are dense enough to allow collective and kinetic behaviors to play roles. The modulation dynamics of pairs is further demonstrated after transparency has occurred, which shows that the positron average energy, phase-space and angular distributions can be modulated periodically by the stable SW field formed directly by the counter-propagating laser pulses.
VI acknowledgements
This work was supported by the National Natural Science Foundation of China (Grant Nos. 11405083, 11347028 and 11675075), and the Graduate Student Innovation Project Foundation of Hunan Province (Grant No. CX2016B450). W.L. appreciates the support from China Scholarship Council and Young Talent Project of the University of South China. M.C. appreciates the support from National 1000 Youth Talent Project of China.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) E. Weibel, Phys. Rev. Lett. ๐ 2 \bf{2} , 83 (1959).
- 2(2) B. Paczynski, Astrophys. J. ๐๐๐ 308 \bf{308} , L 43 (1986).
- 3(3) M. J. Rees and P. Mรฉszรกros, Mon. Not. R. Astron. Soc. ๐๐๐ 258 \bf{258} , 41P (1992).
- 4(4) E. Waxman, Phys. Rev. Lett. ๐๐ 75 \bf{75} , 386 (1995).
- 5(5) T. Piran, Phys. Rep. ๐๐๐ 314 \bf{314} , 575 (1999).
- 6(6) T. Piran, Rev. Mod. Phys. ๐๐ 76 \bf{76} , 1143 (2005).
- 7(7) P. Mรฉszรกros, Rep. Prog. Phys. ๐๐ 69 \bf{69} , 2259 (2006).
- 8(8) H. Chen e โ t โ a โ l ๐ ๐ก ๐ ๐ et~{}al Phys. Plasmas ๐๐ 22 \bf{22} , 056705 (2015).
