Electromagnetic Burst Generation during Annihilation of Magnetic Field in Relativistic Laser-Plasma Interaction.

We present the results of theoretical studies of formation and evolution of the current sheet in a colliosionless plasma during magnetic reconnection in relativistic limit. Relativistic magnetic reconnection is driven by parallel laser pulses interacting with underdense plasma target. Annihilation of laser created magnetic field of opposite polarity generates strong non-stationary electric field formed in between the region with opposite polarity magnetic field accelerating charged particles within the current sheet. This laser-plasma target configuration is discussed in regard with the laboratory modeling of charged particle acceleration and gamma flash generation in astrophysics. We present the results of 3-dimensional kinetic simulations and theoretical studies on the formation and evolution of the current sheet in a collisionless plasma during magnetic field annihilation in the ultra-relativistic limit. Annihilation of oppositively directed magnetic fields driven by two laser pulses interacting with underdense plasma target is accompanied by an electromagnetic burst generation. The induced strong non-stationary longitudinal electric field accelerates charged particles within the current sheet. Properties of the laser-plasma target configuration are discussed in the context of the laboratory modeling for charged particle acceleration and gamma flash generation in astrophysics.

Introduction

Magnetic reconnection (MR) is a fundamental process in astrophysical and laboratory plasmas which provides a mechanism to convert the magnetic field energy to the energy of plasmas and charged particles accompanying with the topology variation of the magnetic field lines[1-3]. Various impulsive phenomena in space plasmas can be attributed to MR such as solar flares[4-6], coronal mass ejections[7,8], pulsar winds[9-11], open and closed planetary magnetospheres[12-14] and gamma-ray bursts[15-18]. In laboratory plasmas, MR is related to the nuclear fusion plasma dynamics[19-22] and the weakly ionized plasmas[23].

The observations of high-energy gamma rays from the Crab Nebula reported by Agile and Fermi-LAT in 2010 is tightly related to MR[24-27]. According to refs. [28], the underlying mechanism of gamma-flares in the Crab Nebulas is based on the electron acceleration during magnetic reconnection in the relativistic limit. The features of the time dependence of gamma-flares indicate that the electron acceleration occurs in non-stationary electric field. Modeling of such the phenomena under the conditions of terrestrial laboratories is an intriguing problem.

The dynamics of laser produced plasmas has been shown in refs. [29-32] to provide a framework where MR can be investigated in the regimes of interest for astrophysical environment. One of the pioneering experiments on MR studying in laser-plasmas was proposed by Nilson et al.[33] with two nanosecond laser beams and a planar solid target. Li et al.[34] irradiated a thin plastic (CH) foil with two or four 1-ns long OMEGA laser beams. Both experiments observed the plasma jets with keV energy in the reconnection region. Recent experimental works with intense and high-energy laser pulses have shown plasma outflows with keV electrons and plasmoid generation in current sheets formed during reconnection on the time scale of nanoseconds[35]. Fan-like plasma outflows with MeV electrons and a plasmoid were obtained in a similar setup with Al foils by Dong et al.[36]. Lezhnin et al. presented the results of the magnetohydrodynamics simulations of driven magnetic reconnection on colliding magnetized laser-produced plasmas[37]. In the presence of a very powerful (petawatt) short laser pulse MR transits into the collisionless relativistic regime. Several numerical studies have been performed in the relativistic regime. Ping et al.[38,39] reported a fast reconnection driven by two ultra-intense laser pulses using 3D kinetic simulations and the corresponding change of the topology structure was observed. MR driven by nonthermal and relativistic electron beams have been discussed recently in refs. [40-43]. Magnetic reconnections under the extreme condition of QED critical field are proposed in refs. [44,45].

In ultra-relativistic regime, the MR acquires novel features. Due to the relativistic constraint on the electron velocity, the variation of the magnetic field cannot be sustained by the upper limit of the electron current. The magnetic annihilation is principally dominated by the displacement current and the corresponding inductive electric field. This so-called dynamic dissipation of the magnetic field was first proposed by S. I. Syrovatskii[46]. Recent numerical studies with particle-in-cell (PIC) simulations presented clear signatures of the particle acceleration in the magnetic annihilation regime[47-50].

In this paper, we present the results of kinetic simulations on the collisionless relativistic MR regime in the 3D configurations. The purpose is to investigate the electron acceleration via the MR generated electric field. We consider the configuration in which two sub-petawatt short laser pulses interact with the hydrogen plasma target comprising two density steps. The magnetic fields with opposite polarities generated by the laser driven electron current is expected to annihilate in the low density region due to the transverse expansion of the magnetic field. In such the non-adiabatic region where the magnetic field freezing-in breaks, the plasma cannot react collectively anymore. The fast annihilation creates strong electric field accelerating electrons up to high energy.

The paper is organized as follows. First we describe the simulation setup. Section of Electron Current Filaments and Strong Magnetic Field Generation presents the results of the kinetic simulations showing two electron current filaments generated by two laser pulses interacting with underdense plasma target and the associated magnetic field configuration produced by these electron current filaments. The next section is devoted to the description of opposite magnetic polarities merging resulting in the magnetic X-line formation and its evolution to thin current sheet. Then we discuss the excitation of the electromagnetic burst and the corresponding charged particle acceleration. The dynamics of the electrons and the typical trajectories are presented. Reconnection of the magnetic field lines in collisionless plasmas is closely related to the Hall effect resulting in the transverse electric field and current excitation. The Hall effect and the quadruple magnetic field formed in the vicinity of the X-line are considered. The tearing mode-like instability leading to the current sheet break-up into magnetic islands is described. The last part contains discussions and conclusions.

Results

Simulation setup

Two linearly polarized Gasussian pulses with the peak intensity of W cm propagating in parallel along the -axis are focused on the plane , where is the laser wavelength. The optical axes of the two pulses are transversely separated by a distance equal to . The normalized amplitude is , where and are the laser electric field strength and frequency, and are the electron charge and mass, respectively; and is speed of light in vacuum. The pulse duration is fs and the spot size (FWHM) is of . The laser pulses radiate on a hydrogen plasma target with two density plateau. The high density region is , which then linearly decreases to the low density of . Here is the plasma critical density, which is approximately cm for m wavelength laser radiation. By employing a density downramp region, the magnetic field is forced to expand in the lateral direction quickly as discussed below in the section of X-line formation in the process of opposite magnetic polarities merging. The second density plateau suppresses the strength of the longitudinal electric field arrising due to the electric charge separation effect so that the inductive electric field effect can be clearly distinguished. The details of the simulations and the code information are presented in the Methods.

Electron current filaments and strong magnetic field generation

The laser pulses generate plasma channels in the underdense hydrogen target since the power of the pulses is higher than the relativistic self-focusing threshold, , where . Balancing the electron energy gain from the charge separation field and the laser field, one obtains the radius of the plasma channel as: , where is the amplitude of the laser pulse vector potential in the channel[51]. Due to the self-focusing effect of the laser field, becomes higher than the initial dimensionless amplitude . The corresponding channel radius in this case is , which is well consistent with the simulation results as shown in Fig. 1(a). The electron density distribution (at , is the laser period) on the bottom plane of Fig. 1(a) presents a double channel structure. In each channel, there is an electron beam trapped and accelerated by the wakefield. The strength of the wakefield is given by Sprangle et al.[52] as: . The maximum amplitude reaches so that the trapped electrons can be accelerated to high energy in a short distance. The energy spectra of the total electrons at , and are plotted in the right plane. The corresponding energy spectra of the high energy electrons, which are trapped and accelerated by the wakefield, are presented in the back plane. The peak energy increases from to within about . Each accelerated beam contains a large charge of the relativistic electrons ().

Figure 1

(a) The electron density distribution (for the plane cut at z = 0) at is plotted in the bottom plane. In the right plane, the curves represent the total energy spectrum at (green), (blue) and (red). The energy spectra of the wakfield accelerated electrons are shown in the back plane with the corresponding time. (b) The transparent cloud represents the electron density distribution in 3D. The current density (along z = 0 plane) is projected in the bottom. The distribution of the (along z = 0 plane) is projected on the back plane. The laser intensity distribution (in the center of the pulse) is projected on the left plane.

Figure 1(b) presents the electron density distribution in the 3D space. The cloud represents the walls structure of the plasma channels. The two accelerated high charge electron beams (the bunches inside the channel) generate the parallel currents, which are shown in the bottom plane. The laser intensity distribution is projected on the left plane. Due to the self-focusing effect, the intensity becomes as large as two times of the initial peak intensity . The currents produce magnetic fields according to Ampere-Maxwell law. The z-component of the azimuthal field () is projected in the back plane, which displays two dipole structures (the positive and negative polarities generated by one electron beam). The strength of the magnetic field can be estimated from Ampere-Maxwell law, By assuming a quasistatic condition of , one obtains that , which reaches about .

X-line formation in the process of opposite magnetic polarities merging

The radius of the magnetic dipole structure depends on the radius of the plasma channel . In Fig. 2(a), the distribution of the longitudinal components of the current density () is shown (zooming in the region of the upper half of the simulation box). The magnetic field generated by the forward accelerated electrons () is shielded by the return electron currents () which move along the channel wall. The return currents provide a magnetic field with opposite direction in the region out of the channel. As mentioned in the section of Electron current filaments and strong magnetic field generation, the radius of the channel is a function of the local plasma density, and . Therefore, the channel expands with the downramp of the density according to . The expansion velocity depends on the gradient of the density as , where is the forward propagating velocity of the channel structure, i.e. the group velocity of the laser pulse propagating in the plasma.

Figure 2

(a) The current density distribution (z = 0 plane) at is shown. The dashed arrows indicate the direction of the electron motion. The solid arrows show the radius of the channel structure. and represent the forward and return currents, respectively. (b) The schematic of the magnetic field expansion. (a–e) are the surfaces of the magnetic field in the plane.

With the expansion of the channels, the size of the magnetic dipoles also increase in the transverse direction. As illustrated in Fig. 2(b), the return currents shift from the red ones () to the green ones (). The black curves represent the amplitude of the magnetic field generated by the accelerated electron beam (), i.e. the amplitude of the z-component of the field (). The dashed parts are shielded before the expansion of the channel. Due to the outward shift of the return currents, the dashed parts also manifest. The contour lines for the magnetic field on the (y,z) plane also shows the expansion from the solid to the dashed region.

The distribution in the plane of at and are shown in Fig. 2(c–e). It’s clear that the boundaries of the magnetic field extend from to , while the strength of the magnetic field decreases. It is due to the outward shift of the return currents reducing their contribution to the magnetic field. The curves describe the profiles of , in which the gradient around reflects the variation of the magnetic field (arising from the expansion and collision of opposite magnetic polarities). In Fig. 2(c) at , the magnetic fields with opposite polarity do not overlap each other in the center and the curve is smooth. However, with the propagation in the density downramp, the opposite polarities start to interacting and the curve become sharp as seen in Fig. 2(e).

When the two regions with opposite polarity magnetic field merge, the so-called X-line structure in MR is formed. Fig. 3(a) depicts the longitudinal currents and the corresponding magnetic field lines in the plane at at . The vectors indicate the direction of the magnetic field lines. Around the two centers of the accelerated electron beams, the magnetic field loops independently belong to the left and right currents respectively. The separatrix surface forms in the transverse cross-section an eight-like curve with a null line locating in the middle of the two currents. Zooming into the region of the X-line in Fig. 3(b), the magnetic field lines have the hyperbolic structure in its vicinity as indicated by the green dashed line. This field can be expressed in terms of the local vector potential as: . Here is the coefficient which describes the different expansion distance in and direction. The magnetic field expand freely in the direction. However the expansion in direction is limited by the distance between the two currents. The X-line is the place where redistribution of magnetic fluxes occurs, which changes the connectivity of field lines. Along the separatrix, a strong electron current is formed. This thin and wide current sheet is one of the signatures of MR due to the particle acceleration via the variation of the magnetic field, which will be discussed in the next section. From the transverse current distribution shown in Fig. 3(c,d), one can find the two channels are still expanding to merge. This expansion pushes more magnetic field lines to annihilate and reconnect in the current sheet region.

Figure 3

(a) The curves are the contour magnetic field strength and the vectors represent the magnetic field lines with directions. The background is the current density distributions in the plane of at . (b) Zooming to the zeroth-point region in (a). The dashed green lines represents the hyperbolic structure. (c,d) are the distributions of and .

Electromagnetic burst and particle acceleration

As mentioned in the previous section, the current sheet is formed due to the particle acceleration in the region where the magnetic fields with opposite polarities annihilate. The particle acceleration in the X-line of MR has been interpreted by several mechanisms including fast shock wave[53], second type Fermi acceleration by turbulence[54], the reflection in magnetic islands[55] and the direct acceleration by reconnection electric field[56]. Here we proposed the regime of dynamic dissipation of the magnetic field and electromagnetic burst in the ultra-relativistic regime. Recalling the Ampere-Maxwell law 1, we see that the variation of the magnetic field is sustained by the conduction current with the current density equal to and the displacement current with the density . Due to the relativistic constraint on the particle velocity, the conduction current density has its limit as , where is the electron density of the current sheet. The local plasma density is relatively low in the downramp region of the target so that the variation of the magnetic field cannot be balanced by only the conduction current. In this case, the displacement current come to play an important role, which is actually the growth rate of the electric field (as noted above, ). The magnetic field distribution on the surface is presented in Fig. 4(a) for time at . The contributions of different terms in Ampere-Maxwell law along the current sheet () are plotted. In the region of , the conduction current (, green) is almost negligible and the variation of the magnetic field (, black) is identical with the displacement current (, red). With the increase of the displacement current, a strong longitudinal electric field (blue) is induced. The longitudinal electric field grows as a result of the onset of the displacement current. The peak of the displacement current propagates along the x-direction. Therefore the longitudinal electric field has its peak value behind the summit of the displacement current. This inductive electric field accelerates the electrons in the backward direction. The magnetic field energy is then transferred to the kinetic energy therefore it is called as dynamic dissipation by S.I. Syrovatskii.

Figure 4

(a) The surface represents the distribution at in the plane. The terms in the Ampere-Maxwell law along are plotted as the curves. (black), (green) and (red) are normalized to . The longitudinal electric field (blue) is normalized to . All the values are transversely averaged inside the current sheet (). (b) The real trajectories of the electrons ejected from the vicinity of X-line (from to ). The distribution of in the plane of at is on the bottom. The curves in the bottom plane represent the contour of magnetic field strength. (c) The energy spectra of the electrons initially locating in the current sheet region. The inset shows the quasi-mono-energetic electrons spectrum.

The dynamics of the particles in the vicinity of the X-line has been discussed in ref. [57]. Here we employ the main conclusions suitable for our case. The electric and magnetic field around the reconnection region can be approximately described as: Then the dynamic equations of the particle acceleration are:

In the ultra-relativistic case, the scale-length (Larmor radius) characterizing the electron trajectory is and the characteristic time of electron passing through the X-line region is . With R and T being the size and time scale related to the non-adiabatic region in the vicinity of the X-line, where the charged particles are not magnetized. Within this region the electron trajectory is given by the solutions of Eqs. 4–6) which read as Here and are the ordinary and the modified Bessel function of zeroth order. From Eqs. 7–9), one can see the electron trajectories oscillate in the -direction and exponentially expand in the -direction. The typical ejected electrons accelerated by the inductive electric field through the vicinity of X-line are selected and the corresponding real trajectories are plotted in Fig. 4(b). The trajectories have clear one period oscillation in the -direction and quick expansion in the -direction, which is well consistent with the theoretical description. The bottom surface of Fig. 4(b) shows the longitudinal electric field distribution in the plane of at . The curves are the contour of the magnetic field strength. The bright spot in the center shows the strong inductive electric field induced by MR, where the magnetic field is almost zero. The strength of the inductive electric field reaches about . The inductive electric field region moves in the forward direction with the propagation of the opposite magnetic polarities. In this case, the backward accelerated electrons only experience a short range in the acceleration phase and then are ejected away from the field. The relative small region and the strong field strength decide that the electrons experience almost the same interaction time. Therefore the energy of the ejected electrons can be estimated as , here is the distance in which the electron experiences the field and can be approximately equal to the size of the X-line region as . Then the characteristic energy is about . Here the energy spectrum of the electrons which are initially localized in the current sheet is provided in Fig. 4(c). An energy peak appears at around which is expected by our estimation. The mono-energetic radiation in the astrophysics is one of the difficulties in explaining by other acceleration mechanisms. Here we found by the MR induced particle acceleration, it is natural to obtain the mono-energetic beam.

Pattern of the longitudinal magnetic field induced by the hall effect

The Hall effect manifests the basic properties of the magnetic reconnection in collisionless plasmas (e.g. see refs. [58-60]). The Hall effect is closely related to the in the transverse electric field and current excitation, seen in the quadruple magnetic field formed in the vicinity of the X-line. With the magnetic field annihilation and the inductive field increasing, it generates in-plane currents due to the decoupling between electrons and ions. The plasma transversely drifts according to the distribution of the magnetic fields and the inductive electric field with the drift velocity as shown in Fig. 5(a). In the vicinity of the magnetic null line, the currents generate the longitudinal magnetic field with the feature of a characteristic quadrupole. It is also such the quadruple patterns are distinctly seen in Fig. 5(b).

Figure 5

(a) The schematic of the field and drift distribution in the vicinity of the magnetic null point. (b) The longitudinal magnetic field () distribution in the plane of at .

Tearing-like mode instability in current sheet and magnetic islands

The current sheet is unstable and it may break-up into filaments due to the development of the tearing mode instability. The corresponding schematics are shown in Fig. 6(a–c). In Fig. 6(d–f), we present the evolution of the current density in the plane with the magnetic field lines in the 3D simulation results. The snapshots correspond to the time at , and . The current sheets clearly extend in the -direction accompanying with filamentation. The electrons accelerated by the inductive electric field via magnetic field annihilation generate a current in the X-line region as shown in Fig. 6(a,d). The trajectories of the accelerated electrons experience an expansion in the -direction as discussed in the section of Electromagnetic burst and particle acceleration. The corresponding magnetic field produced by this current changes the local field topology. The X-line now splits into two symmetric zeroth lines in the -direction as in Fig. 6(b). This transverse expansion continuous and finally form a thin but wide current sheet. Accompany with the current sheet formation, the tearing mode instability leads to the current filamentation and pinching. As shown in Fig. 6(c,f), the current breaks into separated pieces and secondary magnetic islands are formed in the current sheet. Magnetic islands formation is often related with the maxima and minima of the plasma density and the magnetic field separatrix and X-lines[61]. However, in incompressible plasmas where the density does not change but the magnetic islands can be formed, as in the tokamak plasmas where the plasma incompressibility is provided by strong toroidal magnetic field[20] or in the EMHD magnetic reconnection where the electron density remains unchanged due to electric charge separation electric field effects[59]. Since in the case under consideration the charge separation electric field also contributes to the electron density evolution the relationship between the electron density distribution and the magnetic field configuration is quite complicated. The filamentation and the breaking consistent with the electron trajectories shown in Fig. 4(b).

Figure 6

(a–c) are the constant vector potential surfaces corresponding to the magnetic field in the vicinity of the current sheet presenting the current sheet expansion and the bifurcation of the X-line. (d–f) are the results from the 3D simulations. The corresponding times are , and . The curves are the contour magnetic field strength and the vectors represent the magnetic field lines with directions. The background is the current density distributions in the plane of .

Discussions and Conclusions

In conclusion, we investigate the magnetic reconnection driven by laser-plasma interaction by using the 3D kinetic simulations. It presents the formation and evolution of the current sheet in a collisionless plasma during magnetic field annihilation in ultra-relativistic limit. The accelerated electron beams generated in the gas jet plasma create strong magnetic fields. The annihilation of opposite magnetic polarities is accompanied by an electromagnetic burst generation whose strong non-stationary electric field accelerates the charged particles within the current sheet. It is found that the displacement current plays important role in the ultra-relativistic MR to induce the significant growth of the longitudinal electric field. In the vicinity of the magnetic null line, charged particle acceleration is observed. Since the inductive electric field moves in the forward direction with respect to the propagating of the laser field, the electrons, which are accelerated in the backward direction, experience only an instantaneous kick. Therefore the corresponding electron bunch has a relative small energy spread. The dynamics of the particles accelerated by this field in the region in the vicinity of the X-line has been studied. Narrow energy spectrum electron beam is obtained which will be useful in explaining the radiation spectrum obtained in the astrophysics. One of the intriguing problems standing in astrophysics for a number of years is explanation of the detected gamma-ray spectrum pointing towards a very narrow particle spectrum, which is one of the arguments against the shock-acceleration regime. Our results of the mono-energetic bunch generation provide a clear signature to support the particle acceleration via MR regime. Due to the development of the tearing-like mode instability, the current sheet breaks into separated pieces. It leads to formation and evolution of the secondary magnetic islands in the current sheet. The requirement of the laser energy can be expected to be fulfilled by the upcoming facilities like ELI-Beamlines[62]. The regime proposed can be used for formulating the program of forthcoming experiments, including the research in laboratory astrophysics[30,31,60]. The diagnostics for the backward accelerated electron beam will be one of the signatures for the relativistic magnetic reconnection. Furthermore, the growing of the displacement current generates a low frequency wave which has a characteristic wavelength about . On the extreme low density side where the plasma becomes transparent, such a wave serves as another signature in the potential experiment.

Methods

The simulations are performed with the relativistic electromagnetic code EPOCH[63,64] in 3D cases. Two linear-polarized laser pulses with the peak intensities of incident from the left boundary and propagates along the -axis. The optical axes of the two pulses locate at . The simulation box has the size of and . The mesh size in the simulations is . All the quasiparticles (8 per cell) are initially at rest with a total number of . The real mass ratio of electron and proton () is used in the simulations. Open boundary conditions are applied for both fields and particles. The near critical density hydrogen plasma target has a thickness of in the x-direction and remains uniform in the transverse direction in the region of . The density linearly increases from 0 to in the region and then remains constant for . From to , the density decreases to .