INTRODUCTION
The channeling or boring of an intense laser pulse in plasma and ion acceleration have been investigated extensively (Wilks et al., Reference Wilks, Kruer, Tabak and Langdon1992; Borisov et al., Reference Borisov, Borovskiy, Korobkin, Prokhorov, Shiryaev, Shi, Luk, McPherson, Solem, Boyer and Rhodes1992; Tabak et al., Reference Tabak, Hammer, Glinsky, Kruer, Wilks, Woodworth, Campbell and Perry1994; Zepf et al., Reference Zepf, Castro-Colin, Chambers, Preston, Wark, Zhang, Danson, Neely, Norreys, Dangor, Dyson, Lee, Fews, Gibbon, Moustaizis and Key1996; Pukhov & Meyer-ter-Vehn, Reference Pukhov and Meyer-ter-Vehn1997; Borghesi et al., Reference Borghesi, MacKinnon, Barringer, Gaillard, Gizzi, Meyer, Willi, Pukhov and Meyer-ter-Vehn1997, Reference Borghesi, Campbell, Schiavi, Willi, Galimberti, Gizzi, Mackinnon, Snavely, Patel, Hatchett, Key and Nazarov2002, Reference Borghesi, Kar, Romagnani, Toncian, Antici, Audebert, Brambrink, Ceccherini, Cecchetti, Fuchs, Galimberti, Gizzi, Grismayer, Lyseikina, Jung, Macchi, Mora, Osterholz, Schiavi and Willi2007; Fuchs et al., Reference Fuchs, Adam, Amiranoff, Baton, Gallant, Gremillet, Héron, Kieffer, Laval, Malka, Miquel, Mora, Pépin and Rousseaux1998; Vshivkov et al., Reference Vshivkov, Naumova, Pegoraro and Bulanov1998; Tanaka et al., Reference Tanaka, Allen, Pukhov, Kodama, Fujita, Kato, Kawasaki, Kitagawa, Mima, Morio, Shiraga, Iwata, Miyakoshi and Yamanaka2000; Willi et al., Reference Willi, Campbell, Schiavi, Borghesi, Galimberti, Gizzi, Nazarov, Mackinnon, Pukhov and Meyer-ter-vehn2001; Kodama et al., Reference Kodama, Mima, Tanaka, Kitagawa, Fujita, Takahashi, Sunahara, Fujita, Habara, Jitsuno, Sentoku, Matsushita, Miyakoshi, Miyanaga, Norimatsu, Setoguchi, Sonomoto, Tanpo, Toyama and Yamanaka2001; Najmudin et al., Reference Najmudin, Krushelnick, Tatarakis, Clark, Danson, Malka, Neely, Santala and Dangor2003; Hoffmann et al. Reference Hoffmann, Blazevic, Ni, Rosmej, Roth, Tahir, Tauschwitz, Udrea, Varentsov, Weyrich and Maron2005; Flippo et al. Reference Flippo, Hegelich, Albright, Yin, Gautier, Letzring, Schollmeier, Schreiber, Schulze and Fernández2007; Laska et al. Reference Laska, Badziak, Gammino, Jungwirth, Kasperczuk, Krása, Krouský, Kubeš, Parys, Pfeifer, Pisarczyk, Rohlena, Rosinski, Ryć, Skála, Torrisi, Ullschmied, Velyhan and Wolowski2007; Lei et al., Reference Lei, Pukhov, Kodama, Yabuuchi, Adumi, Endo, Freeman, Habara, Kitagawa, Kondo, Kumar, Matsuoka, Mima, Nagatomo, Norimatsu, Shorokhov, Snavely, Yang, Zheng and Tanaka2007). In the relativistic-intensity regime, the relativistic ponderomotive force, which is proportional to the gradient of the laser intensity, dominates the laser-plasma interaction (). The laser pulse enters the plasma like a piston, pushing and compressing the plasma in front of it. However, the high-intensity laser-plasma interaction also leads to other nonlinear processes such as laser filamentation and breakup (Schifano et al., Reference Schifano, Baton, Biancalana, Giulietti, Giulietti, Labaune and Renard1994; Tanaka et al., Reference Tanaka, Allen, Pukhov, Kodama, Fujita, Kato, Kawasaki, Kitagawa, Mima, Morio, Shiraga, Iwata, Miyakoshi and Yamanaka2000; Osman et al., Reference Osman, Yu, Hora, Cao, Liu, He, Badziak, Parys, Wolowski, Woryna, Jungwirth, Králikova, Krása, Láska, Pfeifer, Rohlena, Skála and Ullschmied2004), scattering and propagation instabilities (Najmudin et al., Reference Najmudin, Krushelnick, Tatarakis, Clark, Danson, Malka, Neely, Santala and Dangor2003; Bret & Deutsch, Reference Bret and Deutsch2006; Laska et al., Reference Laska, Badziak, Gammino, Jungwirth, Kasperczuk, Krása, Krouský, Kubeš, Parys, Pfeifer, Pisarczyk, Rohlena, Rosinski, Ryć, Skála, Torrisi, Ullschmied, Velyhan and Wolowski2007; Gupta et al., Reference Gupta, Sharma and Mahmoud2007), etc., which can reduce the quality and length of the channel. Efficient channeling of a relativistic laser in higher density plasmas has been found to be rather difficult.
In this paper, instead of a single pulse, we consider the channeling of two or more laser pulses into homogeneous as well as inhomogeneous plasmas. It is found that with suitably chosen duration and interval of the pulses, a train of short laser pulses can channel into the plasma deeper and more stably than a single longer pulse of similar peak intensity and total energy. The improvement can be attributed to the repeatedly applied ponderomotive force in the forward direction, the suppression of laser-driven plasma instabilities due to the intermittent laser-energy cut-offs, as well as the continuous between-pulse channel lengthening by the still evacuating ions.
THE SIMULATION PARAMETERS
We shall investigate the interaction of a train of short, intense laser pulses with dense plasma using two-dimensional (2D) relativistic particle-in-cell (PIC) simulation (Xu et al., Reference Xu, Chang, Zhuo, Chang and Yue2002; Yu et al., Reference Yu, Yu, Xu, Tian, Chen and Wong2007). Circularly polarized laser pulses with Gaussian envelop in both longitudinal (x) and lateral (y) directions are incident from the left vacuum region along the x axis into the plasma layer, where a L = 4 is the laser strength, τ is the pulse duration, and w = 5 µm is the spot radius, respectively. The corresponding laser intensity is thus about 1019 W/cm2, now available in many laboratories. The laser wavelength is λ = 1.06 µm. The simulation box is 60λ along the x-axis and 30λ along the y-axis. The plasma layer is bounded by 10λ wide vacuum regions on both its left and right sides. The spatial mesh contains 1024 × 512 cells with 12.8 million each of electrons and ions. The initial velocity distributions of the plasma electrons and ions are taken to be Maxwellian, with temperatures of 4 keV for the electrons and 1 keV for the ions. The time step of the simulation is 0.033T, where T = 3.5 fs is the laser period. The spatial and time coordinates are normalized by the laser wavelength and period, respectively, and the ion density is normalized by n c, where n c is the plasma critical density. The electromagnetic energy density is E 2 + B 2, where E and B are the electric and magnetic fields normalized by mω0c/e, where e, m, c, and, ω0 denote the electron charge and rest mass, the speed of light in vacuum, and the laser frequency.
SINGLE PULSE PENETRATION
For the purpose of comparison, we shall first consider channeling by a single laser pulse in a 40λ long uniform plasma layer of density n e = 2n c. Figures 1a–1c shows the distributions of the plasma ion density for the pulses with τ = 12.5T, 25T, and 50T, respectively. All three snapshots were taken when the corresponding channel extends deepest in the plasma, i.e., after the laser light has dissipated or reflected. As is well known, the laser ponderomotive force pushes away the electrons as it propagates in plasma, and the resulting charge-separation field then expels the massive ions, creating a channel. Here we are mainly interested in the effect of the forward-streaming ions on channel lengthening. We note that even after the laser action vanishes due to reflection and absorption, the still evacuating ions continue to lengthen the channel, until the forward momentum of the ions is totally spent, and the final channel length is reached. Figure 1 also shows that some of the expelled ions pile up in the plasma near the channel boundary, forming shock-like compressed layers.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160701183731-71061-mediumThumb-S0263034609000160_fig1g.jpg?pub-status=live)
Fig. 1. (Color online) Channeling of a laser pulse into a plasma layer: distributions of the ion density for pulse durations τ = 12.5 T (a), 25 T (b), and 50 T (c), respectively, at the maximum channel lengths. The case (c) shows that a 50 T pulse leads to considerable deterioration of the channel quality without lengthening it, as compared with that of the 25 T pulse.
For the same peak intensity, a pulse with longer duration has more electromagnetic energy. Figures 1a and 1b show that the τ = 25T pulse can channel deeper into the plasma than the τ = 12.5T pulse. However, further increase of the pulse length or laser energy does not increase the final channel length. Instead, instability-driven deterioration of the channel quality occurs, as can be seen in Figure 1c for the τ = 50T pulse. Here the additional laser energy is diverted to deflecting the light energy, creating a branch channel below the original one without lengthening the latter. For the laser and plasma parameters under consideration, we found that the τ = 25T pulse is optimal for channeling. Here most of the laser energy is spent on the latter, without too much excitation of unfavorable instabilities.
TWO-PULSE SYSTEM
We next examine plasma channeling by two short laser pulses. Instead of a single pulse of τ = 50T, we now employ two τ = 25T pulses. The time delay between the peaks of the two pulses is 70T. The other parameters are the same as those for Figure 1. The total energy in the two-25T pulses is about the same as that in the 50T pulse (Fig. 1c). Figure 2 shows the electromagnetic energy and ion densities at t = 69.62T, when the first laser pulse is channeling into the plasma and the second pulse is still far away. We see that at this time, the laser has already started to deflect, causing bending of the plasma channel. Deterioration of the pulse and channel quality indicates that unfavorable instabilities have set in.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160701183731-96204-mediumThumb-S0263034609000160_fig2g.jpg?pub-status=live)
Fig. 2. (Color online) Channeling by two τ = 25 T laser pulses at t = 69.62 T. (a) The electromagnetic energy density and (b) the ion density. Directional instability and channel bending is clearly occurring.
Figure 3 shows the distributions of the electromagnetic energy and ion densities at t = 79.56T, when the action of the first pulse has almost vanished and the second pulse is still outside the plasma. The shock-compressed layers are visible at the sides of the channel but not at its very front, which is still extending because of the inertially evacuating ions. Furthermore, there is no further development of the instabilities because of the cut-off of the laser-energy supply. Figure 4 shows the electromagnetic energy and ion densities at t = 139.23T, when the second laser pulse is acting on the plasma at the front of the channel created by the first pulse. The plasma channel is further extended in the laser direction. Finger-like laser and channel branching appears, but on smaller space scales. Figure 5a shows the ion density at t = 248.63T, when the action of the second laser pulse is also over. Again, the inertially evacuating ions continue to lengthen the channel, until the forward momentum of the ions is totally spent and the plasma channel attains its maximum length. The finger-like branching structure remains small because of the laser-energy cut-off. For comparison, Figure 5b is for maximum channeling by a single pulse with τ = 50T, i.e., the same as Figure 1c. The significant difference in the length and quality of the channels made by a single pulse and two shorter pulses indicates that the multi-pulse scheme is more efficient, since the total laser energy in the two cases is similar. Nevertheless, to justify the multi-pulse scheme, we need to show that the plasma channel can further lengthen if more laser pulses are applied.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160701183731-94366-mediumThumb-S0263034609000160_fig3g.jpg?pub-status=live)
Fig. 3. (Color online) Channeling by two τ = 25 T pulses at t = 79.56 T. (a) The electromagnetic energy density and (b) the ion density. Here the second pulse is still outside the plasma while the light energy of the first pulse has almost vanished. The energy cutoff stops further bending of the channel. Irregular high-density shock-like layers due to ion pile up can be seen at the sides of the channel.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160701183731-75531-mediumThumb-S0263034609000160_fig4g.jpg?pub-status=live)
Fig. 4. (Color online) Channeling by two τ = 25 T pulses at t = 139.23 T. (a) The electromagnetic energy density and (b) the ion density. Instabilities and channel branching still occur, but at much smaller scales. The difference in the penetration distance of the light and the plasma channel is due to the much slower speed of the evacuating ions.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160701183731-97419-mediumThumb-S0263034609000160_fig5g.jpg?pub-status=live)
Fig. 5. (Color online) The ion density in channeling (a) by two τ = 25 T pulses at t = 248.63 T and (b) by a single τ = 50 T pulse at t = 139.23 T, i.e., same as Figure 1c, with the observation times corresponding to that at maximum channel length. The two-pulse scheme can clearly channel deeper into the plasma. The single longer (50 T) pulse leads to considerable deterioration of the channel quality without lengthening it.
MULTIPLE-PULSE SYSTEM
In view of the above results, we now consider the effect of a third τ = 25T laser pulse. The other parameters are the same as in the two-pulse case, except that the interval between the second and third pulses is (corresponding to the optimum case) 90T.
Figure 6a shows the electromagnetic energy density at t = 228.74T, when the third pulse is inside the plasma channel. Its action is similar to that of the second pulse and it further lengthens the channel. Figure 6b shows the distribution of the light energy later (t = 248.63T). We see that the nearly dissipated light pulse has propagated through the plasma slab and exited at its rear surface. Figure 6c shows the ion density at a still later (in order to allow for the inertial evacuation of the ions) time, namely at t = 288.41T. One can see that the 40λ plasma layer has been completely bored through by the laser pulse train. For comparison, Figure 6d shows the ion density at t = 129.29T for a single pulse with τ = 75T. The total energy in the single 75T pulse is about the same as that in the three-25T pulses. In Figure 6d, one sees bending of the plasma channel. The length of the channel in the laser direction made by the 75T single pulse is much shorter than that made by the three-25T pulses.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160701183731-73875-mediumThumb-S0263034609000160_fig6g.jpg?pub-status=live)
Fig. 6. (Color online) Channeling by three τ = 25 T laser pulses. The electromagnetic energy density at (a) t = 228.74 T and (b) t = 248.63 T. The panel (c) shows the ion density at a still later time (t = 288.41 T), when the 40 λ plasma layer has been bored through. For comparison, the panel (d) shows the ion density at t = 129.29 T for a single pulse with τ = 75 T.
In order to see the effect of density inhomogeneity on multi-pulse channeling, we next consider a plasma slab whose density increases linearly from 0 to 7n c between x = 10λ and x = 50λ. The laser parameters are the same as that for Figure 2, and the time interval between the two peaks is 70T. Figure 7a shows the ion density at t = 198.90T (corresponding to maximum channel length in the inhomogeneous plasma) for channeling by two τ = 25T pulses. For comparison, Figure 7b shows the ion density at t = 169.07T for maximum channeling by a single pulse, as in Figure 5b. Although the total energy in the single 50T pulse is similar to that in the two 25T pulses, the length and quality of the channel made by the two-pulse laser train is again superior. We note that the single 50T pulse can only penetrate to a density of n e < 3n c, and the resulting channel is severely bent. On the other hand, the two-25T pulses can efficiently penetrate up to a density of 4.5n c without causing deterioration of the channel quality. In order to see if the channel in the inhomogeneous plasma can be further lengthened, we consider the effect of a third 25T laser pulse. The time interval between the peaks of the second and third pulses is again 90T. Figure 8 shows the ion density at (a) t = 258.57T as the third pulse starts to act on the plasma, and (b) t = 278.47T for maximum channeling. One can see the third pulse can indeed extend the channel to a still higher density region, namely to n e = 6n c, without causing serious deterioration of the channel quality.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160701183731-18951-mediumThumb-S0263034609000160_fig7g.jpg?pub-status=live)
Fig. 7. (Color online) Laser channeling in inhomogeneous plasma with linear density n i = 0 → 7n c across the slab x = 10 → 50λ. (a) The ion density at t = 198.90 T, when the channel longest, for two consecutive pulses with the same parameters as in Figure 2. For comparison (b) shows channeling into the same plasma by a single τ = 50 T pulse whose parameters are the same as in Figure 5b.
CONCLUSION
In conclusion, we have investigated multi-pulsed laser channeling in both homogeneous and inhomogeneous plasmas using relativistic 2D PIC simulations. It is found that two shorter laser pulses can penetrate deeper into the plasma than a single long pulse of similar peak intensity and total energy. This is because, in the two-pulse case, the ponderomotive force acts on the plasma twice and the laser-induced plasma instabilities that can severely deteriorate the channel quality and length are greatly suppressed because of the cutoff of their source energy from the laser. A third laser pulse can further lengthen the plasma channel. One can expect that laser channeling into dense plasma can be improved by employing a suitably designed train of laser pulses. The scheme is similar to precision hole boring of solid materials by multiple short-pulse lasers at the 1014 w cm−2 level (Bogaerts et al., Reference Bogaerts, Chen, Gijbels and Vertes2003; Zeng et al., Reference Zeng, Mao, Mao, Wen, Greif and Russo2006; Wolowski et al., Reference Wolowski, Badziak, Czarnecka, Parys, Pisarek, Rosinski, Turan and Yerci2007), although the physics is somewhat different. There the repeated action of the short pulses allows the material to be photo-ionized and removed layer by layer at the atomic level, without overheating the hole walls and causing uncontrolled particle emission.
ACKNOWLEDGMENTS
This work was supported by the National High-Tech ICF Committee of China, the Natural Science Foundation of China under Grant Nos. 10576007, 10576035, 10675024, 10676010, 10734130, 10775165, 10835003, and 10875158, the Science and Technology Commission of Shanghai Municipality under Grant No. 08PJ14102, National Basic Research Program of China (973 Program) under Grant Nos. 2007CB815101 and 2006CB806004, Japan-Korea–China cooperative project on High Energy Density Sciences for Laser Fusion Energy, and the JSPS Japan-China Core University Program.