1. INTRODUCTION
Laser propagation in plasma depends on the laser properties as well as the plasma density. In low-density plasmas (n 0 ≪ n c, where n c is the critical density), the laser pulse can propagate a considerable distance. In solid-density plasmas (n 0 ≫ n c), on the other hand, it can penetrate only about a skin length and is reflected. In the recent years, the interaction of short and intense laser pulses (Perry & Mourou, Reference Perry and Mourou1994; Umstadter, Reference Umstadter2003) with gaseous (Tajima & Dawson, Reference Tajima and Dawson1979; Wang et al., Reference Wang, Yu, Yu, Xu, Wang and Yuan2009; Reference Wang, Yu, Yu, Lei, Wang, Senecha, Wang, Murakami and Mima2010; Reference Wang, Sheng, Li, Chen, Dong, Lu, Ma and Zhang2011; Sun et al., Reference Sun, Ott, Lee and Guzdar1987; Borghesi et al., Reference Borghesi, Mackinnon, Barringer, Gaillard, Gizzi, Meyer, Willi, Pukhov and Meyer-Ter-Vehn1997), solid-density (Kruer & Estabrook, Reference Kruer and Estabrook1985; Luan et al., Reference Luan, Yu, Yu, Ma, Zhang, Sheng and Murakami2011; Reference Luan, Yu, Murakami, Zhuo, Yu, Ma and Mima2012a; Brunel et al., Reference Brunel1987), as well as near-critical density (Willingale et al., Reference Willingale, Nilson, Thomas, Bulanov, Maksimchuk, Nazarov, Sangster, Stoeckl and Krushelnick2011; Reference Willingale, Nagel, Thomas, Bellei, Clarke, Dangor, Heathcote, Kaluza, Kamperidis, Kneip, Krushelnick, Lopes, Mangles, Nazarov, Nilson and Najmudin1997; Kuznetsov et al., Reference Kuznetsov, Esirkepov, Kamenets and Bulanov2001; Mourou et al., Reference Mourou, Chang, Maksimchuk, Nees, Bulanov, Yu, Bychenkov, Esirkepov, Naumova, Pegoraro and Ruhl2002; Nakamura et al., Reference Nakamura and Mima2008) plasmas has been widely investigated. However, creating near-critical density plasma with steep density gradient and optically probing the laser-plasma interaction region is experimentally challenging (Willingale et al., Reference Willingale, Nilson, Thomas, Bulanov, Maksimchuk, Nazarov, Sangster, Stoeckl and Krushelnick2011; Reference Willingale, Nagel, Thomas, Bellei, Clarke, Dangor, Heathcote, Kaluza, Kamperidis, Kneip, Krushelnick, Lopes, Mangles, Nazarov, Nilson and Najmudin1997). Recently, low-density foam targets have been fabricated using supersonic gas jets and in situ polymerization (Willingale et al., Reference Willingale, Nagel, Thomas, Bellei, Clarke, Dangor, Heathcote, Kaluza, Kamperidis, Kneip, Krushelnick, Lopes, Mangles, Nazarov, Nilson and Najmudin1997), thus making available near-critical density plasmas suitable for experiments.
Electromagnetic (EM) solitary waves (Yu et al., Reference Yu, Shukla and Spatschek1978; Luan et al., Reference Luan, Zhang and Sheng2008; Reference Luan, Yu, Xu, Murakami, Zhuo, Wang, Wang and Wu2012b; Bulanov et al., Reference Bulanov, Esirkepov, Naumova, Pegoraro and Vshivkov1999; Esirkepov et al., Reference Esirkepov, Kamenets, Bulanov and Naumova1998; Reference Esirkepov, Nishihara, Bulanov and Pegoraro2002; Sentoku et al., Reference Sentoku, Esirkepov, Mima, Nishihara, Califano, Pegoraro, Sakagami, Kitagawa, Naumova and Bulanov1999; Weber et al., Reference Weber, Lontano, Passoni, Riconda and Tikhonchuk2005; Farina et al., Reference Farina, Lontano and Bulanov2000; Sanchez-Arriaga et al., Reference Sanchez-Arriaga and Lefebvre2011a; Sanchez-Arriaga et al., Reference Sanchez-Arriaga and Lefebvre2011b) can self-consistently appear when a laser pulse propagates in a near-critical density plasma (Yu et al., Reference Yu, Shukla and Spatschek1978; Luan et al., Reference Luan, Yu, Xu, Murakami, Zhuo, Wang, Wang and Wu2012b). On the ion time scale, they can appear as post-solitons (Yu et al., Reference Yu, Shukla and Spatschek1978; Luan et al., Reference Luan, Yu, Xu, Murakami, Zhuo, Wang, Wang and Wu2012b; Naumova et al., Reference Naumova, Bulanov, Esirkepov, Farina, Nishihara, Pegoraro, Ruhl and Sakharov2001; Bulanov et al., Reference Bulanov and Pegoraro2002; Borghesi et al., Reference Borghesi, Bulanov, Campbell, Clarke, Esirkepov, Galimberti, Gizzi, Mackinnon, Naumova, Pegoraro, Ruhl, Schiavi and Willi2002), in which the EM radiation sits inside a self-generated plasma cavity. Depending on the laser and plasma conditions, the latter can co-move with the laser pulse or be stationary (Yu et al., Reference Yu, Shukla and Spatschek1978; Luan et al., Reference Luan, Yu, Xu, Murakami, Zhuo, Wang, Wang and Wu2012b; Naumova et al., Reference Naumova, Bulanov, Esirkepov, Farina, Nishihara, Pegoraro, Ruhl and Sakharov2001; Bulanov et al., Reference Bulanov and Pegoraro2002; Borghesi et al., Reference Borghesi, Bulanov, Campbell, Clarke, Esirkepov, Galimberti, Gizzi, Mackinnon, Naumova, Pegoraro, Ruhl, Schiavi and Willi2002; Sarri et al., Reference Sarri, Singh, Davies, Fiuza, Lancaster, Clark, Hassan, Jiang, Kageiwa, Lopes, Rehman, Russo, Scott, Tanimoto, Najmudin, Tanaka, Tatarakis, Borghesi and Norreys2010; Zhu et al., Reference Zhu, Wu, Dong, Hong, Teng, Zhou, Cao and Gu2012). Experiments with near-critical density plasmas have demonstrated the existence of well-isolated solitons and post-solitons (Borghesi et al., Reference Borghesi, Bulanov, Campbell, Clarke, Esirkepov, Galimberti, Gizzi, Mackinnon, Naumova, Pegoraro, Ruhl, Schiavi and Willi2002; Sarri et al., Reference Sarri, Singh, Davies, Fiuza, Lancaster, Clark, Hassan, Jiang, Kageiwa, Lopes, Rehman, Russo, Scott, Tanimoto, Najmudin, Tanaka, Tatarakis, Borghesi and Norreys2010; Zhu et al., Reference Zhu, Wu, Dong, Hong, Teng, Zhou, Cao and Gu2012).
This paper investigates EM radiation trapping in self-generated as well as preformed cavities using two-dimensional (2D) particle-in-cell (PIC) simulation. In Section 2, simulations are performed to detect laser trapping in the self-generated cavity when the laser enters near-critical density plasma. The laser pulse is partially reflected and partially transmitted into the plasma. The ponderomotive force of the light pulse modifies the local plasma density and self-consistently creates in the cavity an oscillon (Stenflo & Yu, Reference Stenflo and Yu1989; Reference Stenflo and Yu1996), or a half-cycle EM standing wave. Section 3 considers laser-light trapping in a preformed cavity in near-critical density plasma. In this case, a multi-peaked EM wave is formed in the cavity. A discussion is given in Section 4.
2. LIGHT TRAPPING IN SELF-GENERATED CAVITY
We performed two-dimensional in space and three dimensional in velocity (2D3V) PIC simulations for laser light trapping in near-critical density plasma. A Gaussian laser pulse is incident along the z axis from the left vacuum region into a homogeneous plasma located at 10 < z/λ0 < 50. The computation box is 60λ0 × 30λ0 in size. The spatial mesh contains 1500 × 750 cells, and each cell contains 16 ions and 16 electrons. The ion-electron mass ratio is 1836, and the initial electron density is n = n(0) = Zn i (0) = 0.7, where n i is the ion density. The incident laser strength, pulse duration, and spot size are a L = 0.4, τ = 10T 0, and b = 5λ, respectively.
When the laser pulse is incident on the left vacuum-plasma boundary, it is partially reflected back to the left vacuum and partially transmitted into the near-critical density plasma. As the transmitted light propagates, the laser ponderomotive force, proportional to the gradient of EM energy density, pushes the plasma away from the region where the EM field is strong. Figure 1 is for the EM energy density E 2 + B 2, the electron and ion densities at t = 77T 0, which is well into the ion time scale. We can see in Figures 1a and 1b that the electrons pushed away by the laser light pile up to form an overcritical-density plasma wall around a cavity, which traps the light. The maximum EM energy density occurs at the cavity center, and the EM field vanishes a short distance (of the order of the effective skin depth) into the high-density plasma wall. Figure 1c is for the densities of the EM energy, electrons, and ions along the laser axis. As expected, the electrons and ions behave almost as a quasi-neutral plasma, with only very little charge separation. The density of plasma wall is more than twice the critical density, thereby preventing the EM radiation from leaving the cavity. Clearly, the near-critical density background plasma plays a crucial role in the existence of the localized structure: the laser light is effectively trapped and stopped by the overdense cavity wall induced by its ponderomotive pressure. Otherwise, a laser can propagate through underdense plasma and be reflected by overdense plasma.
The EM pulse trapped inside the self-generated cavity oscillates in time. Figure 2 is for the normalized electric field component E x at t = 102T 0 (1) and t = 103T 0 (2), about half an oscillation period (~2T 0) apart. The trapped half-cycle standing EM wave can be identified as an oscillon (Stenflo & Yu, Reference Stenflo and Yu1996), and the self-generated cavity a caviton (Wong et al., Reference Wong, Leung and Eggleston1977; Cheung et al., Reference Cheung, Wong, Darrow and Qian1982) with overdense plasma wall.
EM radiation can be trapped in an empty cavity surrounded by an overdense plasma wall. Furthermore, the trapped EM radiation should satisfy the vacuum wave equation
where time and space coordinates are normalized by ω0−1k 0−1, respectively. A simple cylindrically symmetric solution of the vacuum wave equation is
where E(r, t) is the radial electric field, the constant C = 0.2 is determined by best-fit of this solution to that of the simulation. Figure 3 shows E(r, t) at t = cos−1 (−1) and cos−1 (1), namely half an oscillation period apart. We see that the behavior of E(r, t) is consistent with that from the simulation, even though the boundary and initial conditions have not been applied.
3. LIGHT TRAPPING IN PREFORMED CAVITY
We now consider propagating a laser pulse in plasma containing a preformed cavity. As illustrated in Figure 4, a laser pulse enters near-critical-density plasma with a preformed vacuum cavity. The cavity radius is R = 2λ0, and its center is at (y, z) = (0, 16λ0).
Laser propagation in plasma depends on the plasma density, or the dielectric coefficient ɛ = 1 − n e/n c, such that ɛ ~ 0 in near-critical density plasma and ɛ ~ 1 inside the preformed vacuum cavity. Figure 5 shows the electric field component E x, and the electron density n e along the laser axis, at t = 21T 0 and 27T 0. The plasma density and the laser parameters are the same as that in the preceding case. As the laser pulse reaches the preformed cavity at t = 21T 0, its EM field is strongly modified because of the abrupt change in the plasma density: at the plasma-vacuum boundary a part of the laser light is back-scattered, enhancing its tail and exciting electron density oscillations with amplitude remaining less than n c. The rest of the laser pulse enters the cavity and interacts with the plasma at its front. A thin enhanced electron layer of maximum density n e ~ 1.2n c is formed in front of the cavity by the compressed electrons, as can be seen in Figure 5d for t = 27T 0. The laser light is reflected, forming a standing EM wave in the slightly modified cavity. However, a part of the laser energy can tunnel through the narrow above-critical-density electron layer, forming in its front an oscillon-caviton structure, similar to that in a plasma without preformed cavity, but much smaller.
Figure 6 shows the electric field component E x and the electron density n e at t = 81T 0, as well as the electron density along the laser axis direction. EM fields are trapped in the two cavities bounded by ponderomotive-force driven overdense plasma walls that prevent the light from leaving. As a result, the radiation can be trapped for rather long times. As aforementioned, the self-generated small cavity is attributed to radiation tunneling from the preformed cavity. The trapped radiation in this self-generated small cavity again undergoes single-peak oscillation. However, in the preformed cavity with radius larger than laser wavelength, EM radiation becomes a multi-peak structure, as shown in Figure 7a. Our simulation shows that at t = 81T 0 about 7% of the incident laser energy is still in the cavities.
Figure 7 shows the normalized laser electric field E x (1) and electron density (2) at t = 151T 0 for a larger preformed cavity of radius R = 4λ0. The other parameters are the same as that in Figure 6. We can see in Figure 8 for the evolution of the EM energy in the preformed cavity that about 4% of incident laser energy is still trapped there at such long times. This is because in the larger cavity EM energy depletion, which takes place only at the cavity boundary, occurs very slowly since there the EM wave intensity is very low. In fact, the laser light can be stably trapped in the preformed cavity for more than 300T 0.
4. DISCUSSION
In this paper, we have investigated light trapping in self-generated and pre-formed cavities during the propagation of a laser pulse in near-critical density plasmas using 2D PIC simulation. A considerable part of the laser light is trapped as an EM oscillon in the self-generated caviton bounded by an overdense plasma wall. In near-critical density plasma with a preformed cavity, however, the laser ponderomotive force only slightly increases the density of the bounding plasma, but sufficient to reflect and trap the light, forming a standing wave structure. However, a part of the laser light can tunnel through the thin overdense plasma layer at the boundary, and self-consistently form an oscillon-in-caviton structure in front of the preformed cavity. Clearly, the self-consistent overdense plasma boundary layer induced by the ponderomotive force plays an important role for the long-term trapping of the EM oscillon as well as the standing EM wave.
ACKNOWLEDGEMENTS
This work was supported by the National Natural Science Foundation of China (Nos. 11174303, 11204329, 11127901), and the National Basic Research Program of China (Nos. 2013CBA01500, 2011CB808100).